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SPECIFIC  PROGRAM  OBJECTIVES 

The  specific  objectives  of  the  project  were  to  design  and  develop  novel  multifunctional  thermal 
interface  materials  (TIMs)  comprised  of  carbon  nanotubes  (CNTs)  and  epoxy  with  a  particular  emphasis 
on  understanding  and  improving  thermal  energy  transport  in  load-bearing  bonded  aerospace  structures. 


To  accomplish  the  goals  of  the  proposal  the  following  specific  tasks  were  defined  and  pursued: 

Task  1:  Fabricate  thin  film  VACNT/epoxy  composites  as  thermal  interface  materials 

Task  2:  Characterize  thermal  conductivity  in  2-D  thin-film  VACNT-epoxy  composites.  Various 
CNT  diameters,  volume  fraction  and  interfacial  characteristics  will  be  investigated. 

Task  3:  Understand  interfacial  thermal  contact  resistance  between  1-D  CNTs  and  2-D  metal 
capping  layers,  e.g.  tin  thin  film;  Characterize  thermal  conductivity  of  individual  metal  capping 
layers;  Characterize  thermal  conductivity  of  epoxy  matrix  material. 

Task  4:  Develop  experimental  procedures  to  characterize  thermal  conductivity  of  individual 
CNTs.  Various  CNT  diameters  and  annealing  conditions  (graphitic  condition  of  the  CNTs)  will  be 
investigated. 

Task  5:  In  support  of  tasks  1  i  4,  pursue  analytical  and  numerical  modeling  to  better  understand 
thermal  transport  in  nanowires  and  tubes  and  the  VAMWCNT  epoxy  composites. 
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1.  INTRODUCTION 


Thermal  management  at  structural  interfaces  operating  in  high  thermal  environments,  such 
as,  supersonic  and  hypersonic  air  and  space  vehicles,  have  initiated  a  demand  for  high  performance 
load-bearing  thermal  interface  materials  (TIMs).  When  two  nominally  flat  surfaces  come  together  at 
at  a  typical  material  interface,  the  solid- solid  contact  area  is  limited  to  1-2%  of  the  apparent  contact 
area  due  to  surface  asperities  [1].  As  a  consequence,  heat  conduction  through  the  contact  junctions 
and  the  medium  surrounding  the  non-contact  area  provide  parallel  paths  for  heat  flow.  Ideally,  TIMs 
are  designed  to  fill  the  space  in  place  of  surrounding  medium,  requiring  them  to  be  highly  thermally 
conducting  and  also  to  effectively  rWetO  the  bounding  surfaces.  However,  real  TIMs  have  a  finite 
thickness  (called  the  bond-line  thickness,  BLT),  and  invariably  involve  dissimilar  material  surfaces 
with  the  bounding  surfaces.  Moreover,  they  do  not  completely  wet  the  surfaces,  resulting  in  gaps 
with  contact  resistances  at  the  two  bounding  surfaces. 

For  a  typical  TIM  with  a  bond- line  thickness  BLT  and  an  effective  thermal  conductivity  kjiM, 
the  total  thermal  resistance  across  an  interface  can  be  written  in  terms  of  the  thermal  interface  contact 
resistances  on  either  side  of  the  TIM,  Rci  and  RC2,  as 

BLT 

RTotal  =  Rc\  +  T - +  Rc2  ■ 

^ TIM 

In  order  to  make  the  TIM  effective,  the  broader  goal  of  thermal  management  is  to  minimize  the  total 
thermal  interface  resistance  across  the  TIM  by  reducing  BLT  and  increasing  kjiM,  and  minimizing  the 
interfacial  contact  resistances  Rci  and  RC2,  at  the  bounding  surfaces  to  the  TIM. 
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Particle-laden  polymers  are  one  of  the  most  prominent  TIMs  used  in  the  electronics 
cooling  industry  today.  As  described  in  Equation  (1),  the  thermal  contact  resistance  of  such 
material  systems  is  not  only  dependent  on  the  thermal  conductivity  of  the  TIM,  but  also  on  the 
bond  line  thickness  (BLT).  Prasher  et  al.  (2006)  [2]  introduced  a  rheology  based  semi-empirical 
model  for  the  prediction  of  the  BLT  of  these  TIMs.  The  BLT  was  modeled  to  depend  on  the  yield 
stress  of  the  particle-laden  polymer  and  the  applied  pressure.  The  BLT  model  was  then  combined 
with  a  thermal  conductivity  model  to  estimate  the  total  effective  thermal  resistance  of  the 
particle-laden  polymer  TIMs  that  includes  factors,  such  as,  base  polymer  (matrix)  viscosity, 
particle  volume  faction  and  shape,  and  particle  boundary  resistance,  etc.  The  analysis  showed 
that  there  exists  an  optimal  filler  (particle)  volume  fraction  at  which  the  total  thermal  resistance 
of  the  particle  laden  TIM  becomes  a  minimum. 

Recent  studies  indicating  relatively  high  intrinsic  thermal  conductivities  in  single-walled  and 
multi-walled  carbon  nanotubes  as  well  as  other  graphite  materials  [3-11],  suggests  that  these  nano- 
structured  materials  and  their  combinations  are  promising  candidate  materials  for  the  development  of 
high-performance  TIMs.  The  earliest  attempts[12]  for  the  use  of  randomly  orientated  SWCNTs  as 
filler  materials  in  oil  suspensions  reportedly  increased  the  effective  thermal  conductivity  of  nanotube- 
in-oil  suspensions  by  2.5  times  over  that  of  the  base  fluid  (matrix)  with  only  1%  volume  fraction  of 
carbon  nanotubes.  Another  study[13]  involving  randomly  orientated  SWCNT  in  an  epoxy  matrix 
reported  enhancements  in  thermal  conductivity  of  nearly  125%,  while  the  vapor-grown  carbon-fiber 
(VGCF)-epoxy  composites  exhibited  an  enhancement  of  45%  in  thermal  conductivity  when 
compared  to  the  neat  epoxy  for  the  same  loading  (1  wt%).  However,  the  maximum  thermal 
conductivity  achieved  by  randomly  orientated  carbon  nanotube  composites  was  still  less  than  1  W/m- 
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K.  This  surprisingly  low  value  has  been  attributed  to  the  substantially  weak  thermal  cross-linking 
between  contacting  nanotubes  [14]. 

In  recent  years,  vertically  aligned  mutiwalled  nanotube  arrays  have  generated  much  interest 
since  by  placing  CNTs  perpendicular  to  and  spanning  the  system  components,  there  is  a  lesser 
number  of  CNT/epoxy  interfaces  in  the  through-thickness  direction,  thus  minimizing  the 
effective  thermal  interfacial  resistance  between  the  CNTs  and  the  epoxy  matrix.  Assuming  the 
epoxy  matrix  does  not  affect  the  heat  transport  through  individual  MWCNTs,  it  is  expected  that 
the  overall  thermal  conductivity  of  the  composite  may  be  enhanced  by  the  intrinsically  higher 
thermal  conductivity  of  the  MWCNTs. 

Thermal  transport  through  vertically  aligned  multi-walled  carbon  nanotube  (VA-MWCNT) 
arrays  [15]  has  been  observed  to  be  highly  anisotropic.  The  longitudinal  diffusivity  was  reported  to 
be  72  times  that  of  the  in-plane  value,  signifying  their  worth  as  1-D  heat  pipes  between  two 
contacting  surfaces.  A  similar  study  over  a  wider  temperature  range  180K  -  300K  [16,  17]  have 
reported  thermal  diffusivity  measurements  along  the  alignment  direction  in  a  VA-MWCNT  array 
(20-50  nm  diameter,  1.64  mm  long)  to  be  25  times  higher  than  that  across  the  transverse  direction. 
Choi  et  al.  (2003)  [14]  reported  enhancement  in  thermal  conductivity  of  SWCNT-epoxy  composites 
by  nearly  300%  at  3  wt%  loading  of  randomly  dispersed  SWCNTs,  and  an  additional  10% 
enhancement  in  thermal  conductivity  after  applying  a  magnetic  field  during  processing  to  improve 
nanotube  alignment.  The  maximum  room  temperature  thermal  conductivity  obtained  after  magnetic 
alignment  during  processing  was  about  6.5  W/m-K,  which  was  much  higher  those  reported  in 


5 


previous  studies.  These  results  indicate  that  vertically  aligned  CNTs  which  are  in  isolation  from  each 
other  have  a  potential  to  serve  as  promising  candidates  as  TIMs. 

With  regards  to  VAMWCNTs,  Ivanov  et  al.  (2006)  [15]  reported  a  thermal  conductivity  of 
5.5  ±  0.7  W/m-K  for  an  epoxy-infiltrated  VA-MWCNT  array  (8+1  vol  %,  2mm  long)  and  6.4  +  0.8 
W/m-K  of  the  same  array  in  air.  Borca-Tascuic  et  al.  (2007)  [17]  reported  a  maximum  thermal 
conductivity  of  3.8  W/m-K  in  a  MWCNT  array-polymer  composite  along  the  alignment  direction  at 
room  temperature  with  2%  volume  fraction  of  aligned  MWCNTs.  Xuejiao  et  al.  (2004)  [18]  and 
Tong  et  al.  (2004)  [19]  proposed  the  use  of  CNTs  along  with  other  conventional  TIM  filler  materials 
in  a  polymer  matrix,  and  vertically  oriented  growth  of  CNTQ  on  both  contacting  surfaces  to  form  a 
CNT-CNT  interface  in  between  them.  Hu  et  al.  (2007)  [20]  characterized  the  thermal  interface 
resistance  of  CNT-CNT  interfaces  using  diffraction  limited  infrared  microscopy.  They  reported  a 
thermal  interface  resistance  of  the  CNT-based  interface  structure  to  be  much  below  expectations  at 
3.8  xlO-4  K-m2/W.  Xu  and  Fisher  (2005)  [21]  experimentally  evaluated  a  minimum  thermal 
resistance  of  2  x  10'5  Km2/W  for  a  dry  joint  of  VA-CNT  array  in  between  Cu  and  Si,  which  is  very 
close  to  that  of  TIM  comprising  indium  and  CNT.  However,  incorporating  a  phase  change  material 
(PCM)  into  a  VA-MWCNT  array  yielded  a  lower  thermal  resistance  of  5.2  x  1 0"6  Km2/W.  More 
recently,  Cola  et  al.  (2007)  [22]  measured  the  thermal  resistance  across  (a)  one-sided  (Si-MWCNTs- 
Ag),  wherein  carbon  nanotubes  are  grown  on  one  of  the  bounding  surface;  and  (b)  two-sided  (Si- 
CNT-CNT-Cu)  MWCNT  arrays.  For  the  one-sided  interface  the  total  interfacial  thermal  resistance 
was  found  to  be  dominated  by  the  thermal  resistance  of  the  CNT-Ag  interface  and  was  measured  to 
be  1.4  x  10'5  Km2/W,  while  the  thermal  resistance  of  the  two-sided  configuration  was  dominated  by 
the  interface  resistance  between  the  free  tips  of  the  mating  CNT  arrays,  and  was  measured  to  be  an 
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order  of  magnitude  lower  at  2  x  10'  Km  /W.  A  potentially  promising  variation  of  the  two-sided 
CNT  interface  was  the  introduction  of  a  thin  deformable  foil  in  between  the  nanotube  ends  of  the 
CNT-CNT  arrays  such  that  the  foil  could  adjust  to  the  deformation  of  CNTs  with  far  field  pressure.  In 
this  way,  using  a  10  pm  thick  Cu  foil  (with  MWCNT  on  both  sides)  the  thermal  interfacial  resistance 
of  a  rough  Cu-Ag  interface  was  observed  to  be  greatly  reduced  [23]. 

Recent  works  by  Shaikh  (2007)  [24]  and  Cola  (2007)  [22]  have  explored  through-thickness 
thermal  conductance  in  systems  comprising  CNTs  on  glass,  copper  and  silicon  substrates.  In  these 
systems  it  was  hoped  that  the  weld-contacts  made  via  the  catalyst  during  synthesis  of  vertical  aligned 
MWCNT  arrays  would  be  better  than  simple  Van  der  Waal  type  bonds.  However,  it  was  found  that 
even  in  the  best  case  scenario,  the  CNT-substrate  contact  resistance  was  still  quite  high.  Tong  et  al. 
(2007)  [25]  and  Sihn  et  al.  (2008)  [26]  used  a  combination  of  MWCNT  array  with  thin  gold  and 
indium  capping  layers  to  improve  the  thermal  conductance  between  the  MWCNT  array  and  the 
bounding  surfaces;  the  thermal  conductivity  of  the  joint  (device)  was  observed  to  increase  by  more 
than  1-2  orders  of  magnitude  compared  to  the  absence  of  the  TIM. 

Encouraged  by  these  findings,  in  this  paper  we  report  the  development  of  vertically- 
aligned  multi-wall  carbon  nanotube  (VA-MWCNT)  array  composites  for  thermal  management  in 
load-bearing  structural  applications.  Unlike  previous  studies  on  the  use  of  VA-MWCNT-based 
thermal  interface  materials  for  essentially  non-load  bearing  thermal  management  applications, 
the  TIMs  of  interest  here  involve  the  use  of  VA-MWCNTs  in  an  epoxy  matrix.  The  epoxy  matrix 
imparts  mechanical  strength  to  these  systems  while  the  VACNTs  provide  avenues  for  high  thru¬ 
thickness  thermal  conductivity  across  a  typical  material  interface.  It  is  assumed  that  the  thermal 
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properties  of  the  aligned  MWCNTs  in  the  epoxy  matrix  are  similar  to  those  of  individual  free 
standing  MWCNTs.  Furthermore,  we  introduce  a  transition  zone  (TZ)  comprising  of  a  tin  thin 
film  at  the  interface  between  the  MWCNTs  and  the  surrounding  material  (Si02)  to  minimize  the 
thermal  impedance  mismatch  and  thermal  resistance  due  to  interfacial  roughness.  Figure  1 
shows  a  schematic  of  our  proposed  Sn-YA-MWCNT/epoxy  TIM  system.  The  overall  thermal 
resistance  of  the  TIM  is  expected  to  be  governed  by  the  interfacial  thermal  contact  resistance 
between  the  bounding  solids  and  the  mating  surfaces  of  the  Sn  YA-MWCNT/epoxy  layer.  These 
interfaces  are  generally  neither  fully  conforming  nor  smooth  and  thus,  may  lead  to  a  significant 
increment  in  the  total  thermal  contact  resistance. 


Low  melting 
point  metal 
thin  film 


Solid  1 


VA  MWCNTs  in 
jn  epox/  matrix 


Solid  2 


Figure  1:  Schematic  of  VA-MWCNT  TIM.  The  TIM  facilitates  heat  transport  between  bounding  solids 
(solid  1  and  solid  2).  A  low  melting  point  (Sn)  thin  film  (green)  fills  the  transition  zones  between  the 
CNT-polymer/solid  interfaces. 


In  particular,  we  obtain  the  thermal  conductivity  of  the  Sn-coated  VA-MWCNT  epoxy  device 
over  a  temperature  range  of  240  K-  300  K  using  the  three  omega  method.  In  order  to  estimate  the 
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thermal  conductivity  of  the  VA  MWCNT  epoxy  composite  as  well  as  the  thermal  interfacial 
resistance  at  the  various  layer  interfaces,  we  characterize  thermal  conductivity  in  individual 
constituents  of  the  VA-MWCNT  composite,  namely  (1)  individual  free  standing  nanotubes  selected 
from  the  same  VA-MWCNT  array  as  that  used  in  the  fabrication  of  the  VA-MWCNT-epoxy 
composite;  (2)  the  epoxy  matrix;  and  (3)  Sn  thin  films  (capping  layer)  of  500  nm  thickness.  These 
measurements  will  be  used  in  conjunction  with  a  1-D  multilayer  thermal  model  based  on  Feldman© 
notation  [27]  to  estimate  the  thermal  conductivity  of  the  CNT-epoxy  composite  as  well  as  the  thermal 
interfacial  resistance  between  the  layers. 

2.  MATERIALS  AND  METHODS 

2.1.  Sample  Processing 

The  VA-MWCNT  arrays  used  (Figure  2a  and  Figure  2b)  in  the  present  investigation  were 
procured  from  the  laboratories  of  Prof  Dai  at  Case  Western  Reserve  University  and  Prof  Shanov 
at  the  University  of  Cincinnati.  In  Prof  Dai©  laboratory,  VA-MWCNT  arrays  (LD  Samples)  are 
grown  via  pyrolysis  of  iron  (II)  phthalocyanine  (FePc)  in  Ar/H2  at  900  °C  on  a  silicon  wafer.  [28, 
29]  These  MWCNTs  are  ~  25-80  nm  in  diameter  and  have  lengths  varying  from  2.25  mm  in  the 
middle  of  the  turf  to  2.9  mm  at  the  edges.  VA-MWCNT  arrays  procured  from  Prof.  Shanov© 
laboratory  (UC  samples),  were  synthesized  using  water-assisted  chemical  vapor  deposition 
(CVD)  [30-32],  The  nanotube  diameter  of  the  water-assisted  CVD  arrays  varied  from  15-45  nm. 
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Figure  2:  SEM  micrographs  of  VA-MWCNT  arrays  prepared  by  thermal  CVD  taken  at  two  different 
magnifications  (a)  740X  and  (b)  6200X  showing  the  region  near  CNT  tips. 

The  epoxy-based  carbon  nanotube  composites  were  fabricated  in  our  laboratory  by 
immersing  the  MWCNT  array  into  a  solution  of  Epon  862  epoxy,  EPICURE  curing  agent  W,  and 
an  acetone  solvent.  Prior  to  immersion  of  the  YA-MWCNT  array,  the  Epon  862  epoxy  is 
ultrasonicated  for  approximately  8  minutes.  The  solution  is  poured  onto  the  YA-MWCNT  array 
and  then  spin-coated  to  allow  the  epoxy  to  infiltrate  the  array.  The  volume  fraction  of  nanotubes 
in  the  composite  is  estimated  to  be  8%.  Degassing  is  performed  under  high  vacuum  (30  in  of  Hg) 
to  remove  the  bubbles  generated  during  mixing.  After  casting  the  MWCNT/epoxy  composite, 
both  surfaces  are  cut  and  lapped.  The  surfaces  are  then  sequentially  polished  at  100-150  rpm 
with  15,  6,  1  and  0.1  £m  diamond  abrasives  while  applying  a  5  N  constant  force. 

Atomic  force  microscopy  (AFM)  topology  measurements  (Figure  3d,  3e,  and  3f)  indicate 
that  the  polishing  process  reduces  the  roughness  of  the  samples  from  nearly  1pm  to  100  nm  or 
less.  Following  polishing,  the  sample  is  ethanol  washed  and  air-dried.  The  MWCNT  tips  are 
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exposed  from  the  epoxy  by  reactive  ion  etching  (RIE)  using  O2  plasma  with  13.5  MHz  125  W 
RF  power  for  6-8  min.  Once  the  tips  are  exposed,  a  thin  layer  of  tin  of  thickness  ~  500  ±  50  nm 
is  deposited  using  RF  sputtering  at  2  x  10'7  Torr. 
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Figure  3:  (a)  AFM  Phase  imaging  of  MWCNT-Epoxy  composite  before  plasma-etching;  (b)  Tapping 
mode  AFM  reveals  CNT  tips  exposed  after  plasma-etching  on  composite;  (c)  3D  view  of  roughness 
characteristic  on  top  surface  of  epoxy  before  polishing  using  AFM;  (d)  Surface  topology  of  epoxy 
composite  before  polishing  indicating  roughness  600nm-lpm;  (e)  Surface  topology  traced  by  AFM  tip  on 
composite  sample  after  polishing  indicating  roughness  of  the  order  ~  lOOnm. 

2.2.  The  ‘T’  Type  Wollaston  Probe  Method  for  Individual  Nanotubes 

The  thermal  conductivity  of  an  individual  MWCNT  is  measured  using  a  three-omega- 
based  Wollaston  T-Type  probe  inside  a  scanning  electron  microscope  (SEM).  Details  regarding 
these  measurements  are  presented  in  APPENDIX  1  of  this  report.  Figure  4  depicts  the  Wollaston 
probe  wire®  temperature  profile  (a)  before  the  specimen  is  placed  in  contact,  and  (b)  following 
contact  with  the  specimen.  The  drop  in  spatially  averaged  temperature  causes  a  reduction  in 
electrical  resistance,  and  thus  a  measureable  voltage  decrease.  The  resistance  response  of  the  probe 
wire  is  determined  from  the  voltage  response  since  the  heating  current  amplitude  is  known. 

The  probe  wire®  temperature  response  is  first  calibrated  by  measuring  the  electrical 
resistance  versus  power  input  in  the  absence  of  the  sample.  Subsequently,  once  a  sample  is  attached, 
the  change  in  the  measured  electrical  resistance  versus  power  input  of  the  probe  wire  is  correlated  to 
a  decrease  in  the  probe  wire®  average  temperature  rise  due  to  the  heat  flux  into  the  sample.  Since  the 
opposing  end  of  the  sample  is  maintained  at  the  base  (ambient)  temperature,  in  effect,  the 
temperature  drop  and  heat  flux  into  the  sample  can  be  estimated,  and  the  sample®  thermal  resistance 
is  determined. 
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Figure  4:  Schematic  of  the  Wollaston  probe  wire.  The  probe  wire  is  Joule  heated  with  a  power  of  Qrms, 
by  the  low  frequency  current,  Iiwrms ■  C Q(xJ  is  the  temperature  response  of  the  probe  wire  prior  to  making 
contact  with  the  sample  of  unknown  thermal  conductivity.  Following  contact  with  the  sample,  the 
temperature  response  of  the  probe  wire  is  6B(x).  The  measured  third  harmonic  voltage,  f';„  rms,  is  a 
function  of  the  probe  wired  thermal  resistance,  zero  current  electrical  resistance  Reo,  Qrms,  and  the 
thermal  resistance  of  the  sample.  The  sampled  thermal  conductivity  is  determined  by  measuring  the 
sampled  thermal  resistance  and  dimensions. 


Instead  of  DC  Joule  heating  and  voltage  measurements,  the  probe  wire  is  Joule  heated 
with  a  low  frequency  sinusoidal  AC  current,  I(t)=Iicosin(cot),  where  Ii(0  is  the  current  amplitude. 
The  initial  electrical  resistance  of  the  probe  wire  is  Reo.  Joule  heating  at  2co,  Q(t),  drives  in  phase 
temperature  oscillations,  0  (t),  which  are  proportional  to  Q(t),  by  a  thermal  transfer  function  Z0. 
The  2a)  temperature  oscillations  create  oscillations  in  electrical  resistance  by 
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Re  (/ )  =  Reo  [j  +  <2#  (7 )] ,  where  U  is  the  measurable  temperature  coefficient  of  resistance  of  the 

probe  wire.  Since  the  heating  current  is  I(t)=Iimsin(a>t),  the  voltage  response  of  the  probe  wire  is 
expected  to  be  a  combination  of  both  the  lco  and  3co  frequency  harmonics  of  the  initial  heating 
current  and  is  given  by, 


=  h.K*\aZcR 2JI 


sin  cot  - 


\aZjlRl 


sin  3  cot 


(2) 


A  Lock-in  Amplifier  is  used  to  measure  the  specific  RMS  value  of  the  three-omega 
component.  The  measurable  3  co  RMS  voltage  response  from  Equation  (2)  may  be  re-written  as 


^3 a>,RMS  ~  2  a^oha,RMS^eoQRMS 


(3) 


where  Qrms  f  I2  iw,rms  Reo. 

If  Re3co,RMS  =  V3(0, rms/1 i co,  rms,  Z0  is  determined  from  the  measured  slope  of  SRe3o)  RMS  /  5Qms 

.  The  thermal  resistance  ratio  of  probe  to  the  sample,  Rth,p/Rth.s,  is  embedded  inside  the  thermal 
transfer  function  Z0.  By  modeling  the  temperature  response  of  the  probe  wire  as  a  function  of 
RthiP/Rth,s,  a  theoretical  value  of  Z0  is  compared  to  the  experimental  value.  Initial  calibration  in  the 
absence  of  the  sample  is  used  to  first  determine  Rth,p,  and  a  subsequent  measurement  in  the  presence 
of  the  probe  wire  is  used  to  determine  Rth,s.  The  samples  thermal  conductivity  is  calculated  by  ks  = 
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Ls/Rth,s  As,  where  the  sample  cross  sectional  area  As  and  Length,  Ls  is  measured  using  the  SEM 


beam. 


2.3.  Three  Omega  Method  for  Thin  Films  and  Composites 

Similar  to  the  measurements  performed  on  individual  MWCNTs,  the  3 co  method  is 
performed  by  passing  a  current  with  angular  frequency  co  through  a  heater/sensor  micro- 
fabricated  on  the  sample.  The  heating  current  creates  Joule  heating  and  thus  a  temperature 
oscillation  at  a  frequency  2<x>.  Subsequently,  the  metal  heater/sensor Q  electrical  resistance  also 
oscillates  at  2  co,  causing  a  measureable  3  co  voltage  oscillation  detected  using  a  lock-in  amplifier. 
The  3  co  voltage  oscillations  are  used  to  determine  the  sample©  thermal  properties  [33,  34],  The 
frequency  range  for  the  experiment  is  chosen  by  examining  the  relationship  between  the  film 
thickness,  the  thermal  penetration  depth  of  interest,  and  heater  width.  Large  heater  widths  (when 
compared  to  sample  or  film  thickness)  produce  a  one-dimensional  heating  profile,  thereby 
providing  information  about  the  cross-plane  thermal  conductivity. 

For  electrically  conductive  samples  a  thin  insulating  film  must  be  deposited  prior  to  the 
metal  line  deposition  to  provide  electrical  isolation.  For  this  work,  a  thin  film  of  350  nm  Si02 
was  deposited  using  low  temperature  plasma  enhanced  chemical  vapor  deposition  (PECVD). 
The  Si02  serves  as  an  insulator  for  all  thin  films  and  composite  samples  evaluated  in  the  present 
study.  For  conducting  the  experiments,  an  approach  similar  to  the  differential  three  omega  [35- 
37]  method  was  utilized  to  measure  the  thermal  conductivity  of  a  -500  nm  thin  Sn  film,  bulk 
EPON-862  epoxy  film,  and  VA-MWCNT-epoxy  composite  samples,  with  the  same  heater 
material  (Aluminum)  and  width  on  them.  Additionally,  in  order  to  restrict  the  error  introduced 
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by  various  approximations,  the  numerically  obtained  exact  solution  of  the  film  on  substrate 
system,  as  described  by  Kim  et  al.  [38]  ,  Borca-Tascuic  et  al.  [35]  and  more  recently  by  Tong  et 
al.  [37],  has  been  used. 

The  analytical  procedure  for  interpreting  the  three  omega  experiments  is  based  on  an 
earlier  work  by  Cahill  and  Pohl  [39],  and  requires  the  assumptions  that  the  heater  be  considered 
of  infinite  length  and  that  the  sample  be  effectively  semi-infinite.  Additionally,  the  method 
requires  the  thermal  waves  to  penetrate  a  depth  greater  than  the  heater  width  into  the  sample, 
satisfying  the  approximation  of  an  infinite  line  source.  The  three-omega  method  is  applied  to  a 
film  on  a  substrate  type  configuration,  by  employing  either  the  slope  method  or  the  differential- 
3  co  approach.  For  these  methods,  the  thin  film  is  assumed  thin  enough  such  that  essentially 
one-  dimensional  Fourier  conduction  model  can  be  applied  for  analysis  and  interpretation  of  the 
experimental  data.  Specific  conditions  that  must  hold  to  ensure  accurate  analysis  are  that 
df  >  2b  »  Lf  ,  i.e.,  the  thermal  penetration  depth,  d f ,  must  be  greater  than  the  heater  width  and 

be  much  larger  than  the  film  thickness,  L  .  For  the  experimental  investigation  discussed  herein, 

the  thermal  penetration  depth  and  the  heater  width  were  greater  by  at  least  an  order  of  magnitude 
than  the  film  thickness  to  ensure  the  validity  of  these  approximations.  Additionally,  the 
thickness  of  the  underlying  substrate  L  must  be  such  that  the  condition  L  »  df  »  2b  is 

maintained  to  avoid  back  surface  reflection  and  assure  the  validity  of  the  line  source  and  semi¬ 
infinite  substrate  approximations.  The  thermal  conductivity  of  the  3  mm  epoxy  (EPON-862) 
sample  is  measured  using  the  analytical  solution  mentioned  in  Cahill  (1990)  for  a  narrow  line 
source  on  a  semi-infinite  substrate.  The  thermal  conductivity  of  a  500  nm  thin  tin  film  is 
measured  by  comparing  the  total  amplitude  and  phase  signal  of  temperature  oscillations  with  a 
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theoretically  obtained  solution  in  thin  layer  limit  using  Tong®  (2006)  [37]  two  layer  model. 
However,  for  the  experimental  work  on  VA-MWCNT  composites,  a  three-layer  model  based  on 
Feldman®  algorithm  (Feldman,  1999)  has  been  derived  and  is  used  to  interpret  the  experimental 
results.  Additional  details  about  the  three  omega  technique  can  be  found  from  various  references 
[33-39],  and  are  therefore  not  discussed  in  detail  here. 

The  main  components  of  the  3  co  set-up  used  for  this  work  are,  (i)  a  cryostat  (Janis 
Research  Model:  CCS-400H/204);  (ii)  a  temperature  controller  that  serves  to  regulate  the  sample 
temperature;  and  (iii)  a  lock-in  amplifier  to  detect  the  voltage  response  of  the  heater  (Fig.  5).  The 
cryostat  is  capable  of  operating  between  the  temperatures  10K  T  500K.  An  RV-8  rotary  vane 
pump  capable  of  developing  a  vacuum  equal  to  10‘4  torr  or  less  is  utilized.  A  silicon  shadow 
mask  is  batch  fabricated  using  conventional  lithography  technique.  The  aluminum  metal 
heater/sensor  line  is  magnetron  sputtered  through  this  micro-fabricated  shadow  mask.  The 
temperature  coefficient  of  resistance  of  the  heater/sensor  line  is  measured  prior  to  every 
experiment.  The  set-up  has  been  automated  using  LabVIEW  8.5.  The  bits  of  the  multiplying 
DAC  (AD  7541  KN)  are  set  GDnOand  GDffOto  balance  the  first  harmonic  signals  from  the  sample 
and  reference  resistor.  Moreover,  the  LabVIEW  program  helps  in  balancing  the  first  harmonic 
voltage  at  the  reference  resistor  and  the  sample  resistor  to  a  greater  degree  of  precision  and 
subsequently  extract  third  harmonic  voltage  signals  in  order  to  determine  the  temperature  rise  in 
the  heater/sensor. 
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Figure  5:  (a)  Schematic  of  Cryostat  Setup @  Nanomechanics  Lab,  Case  Western;  (b)  Composite  sample 
loaded  on  a  cryostat  after  laying  microfabricated  heaters  on  top;  (c)  Microfabricated  25pm  wide  heater 
with  250pm  x  250pm  pads  attached  for  wire-bonding;  (d)  Schematic  of  conventional  3-Omega 
measurement  for  film  on  substrate  configuration  used  for  measuring  tin  films  (e)  VA-MWCNT  array  on 
silicon  wafer;  and  (f)  Schematic  of  3-Omega  configuration  utilized  for  measuring  VA-MWCNT 
composite. 


3.  RESULTS  AND  DISCUSSION 

An  understanding  of  thermal  transport  in  a  VA  MWCNT-epoxy  based  TIM  in  an 
adhesively  jointed  structure  must  involve  thermal  conductivity  measurements  in  the  individual 
components,  including  the  epoxy  matrix,  the  thin  tin  capping  layer,  and  the  VA-MWCNT-epoxy 
composite.  Additionally,  where-ever  possible,  the  thermal  measurements  are  performed  at  lower 
temperatures  so  as  to  investigate  the  effects  of  increased  phonon  mean  free  path  lengths  on 
thermal  conductivity. 
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3.1.  Thermal  Conduction  in  individual  Carbon-nanotube 


In  the  present  study,  thermal  conductivity  of  individual  multiwalled  carbon  nanotubes 
taken  from  the  same  sample  batch  array  UC01,  are  measured  at  room  temperature.  The  average 
thermal  conductivity  of  an  individual  nanotube  from  the  UC01  batch  is  ~  60  W/m-K.  In  a  recent 
study  [40],  the  authorsO  have  measured  the  thermal  conductivity  of  individual  multiwalled 
nanotubes  from  both  heat-treated  and  non-heat-treated  sample  groups  showing  an  approximate  5 
fold  increase  in  thermal  conductivity  with  heat  treatment.  Heat-treatment  was  performed  at  3000 
°C  for  20  hours  in  an  argon  gas  atmosphere.  Moreover,  for  one  of  the  heat-treated  samples 
thermal  conductivity  as  high  as  730  ±  153  W/m-K  was  obtained. 

Qualitative  assessments  of  residual  amorphous  carbon  and  defects  were  made  by 
observing  the  Raman  peak  intensities  of  the  D  and  G-bands.  The  D  band  (defect  band)  is  known 
to  be  activated  by  both  carbonaceous  impurities  with  sp3  bonding  and  fragmented  sp2  bonds  [41, 
42],  both  being  features  of  MWCNT  defects.  The  G  band  is  activated  by  contiguous  sp2 
bonding,  /.  e. ,  a  high  degree  of  crystallinity  [41,  42],  Therefore,  when  comparing  the  two  sample 
groups,  a  smaller  DIG  band  ratio  is  representative  of  a  sample  with  fewer  defects  and  a  higher 
degree  of  graphitization.  The  authors  have  shown  that  significant  increase  in  thermal 
conductivity  had  a  strong  correlation  to  the  reduced  DIG  Raman  peaks.  The  average  Raman  D/G 
ratio  for  heat-treated  samples  was  0.21  ±  0.04,  while  the  non-heat-treated  group  had  an  average 
D/G  ratio  equal  to  0.69  ±  0.15.  Note  that  the  UC01  samples  are  not  heat-treated  and  Raman 
scattering  reveals  a  D/G  ratio  of  0.84  (Figure  6b).  The  reduced  value  of  thermal  conductivity 
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measured  in  the  UC01  sample  correlates  well  with  the  non-heat-treated  sample  group  with  higher 


D/G  Raman  ratios. 


Thermal  conductivity  is  undoubtedly  expected  to  improve  with  increased  crystallinity  and 
a  decreased  density  of  phonon  scattering  sites.  Andrews  et  al.  (2001)  [43]  had  observed  marked 
structural  and  chemical  improvements  in  MWCNT  sample  quality  with  heat-treatment.  It  was 
observed  that  annealing  temperatures  reduced  the  residual  Fe  catalyst  (essential  for  CVD  growth 
of  MWCNTq)  from  7.1%  by  wt.  to  0.1%.  Moreover,  it  was  reported  that  the  temperatures  as 
high  as  3000°C  healed  the  sidewall  defects  evident  from  the  reduced  inter-layer  spacing  between 
grapheme  shells  confirmed  by  x-ray  diffraction  (XRD)  and  transmission  electron  microscope 
(TEM). 


Wavenumber  (cm'1) 


Figure  6:  (a)  SEM  micrograph  of  thermal  conductivity  experiment  on  an  individual  MWCNT  selected 
from  the  UC01  sample  group,  (b)  Raman  intensity  with  excitation  wavenumber  for  D/G  Raman  ratio  of 
the  UC01  non-heat-treated  sample. 
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3.2.  Thermal  transport  in  Sn  thin  film 


In  the  present  study,  a  transition  zone  comprising  of  a  Sn  thin  film  is  introduced  at  the 
interface  between  the  MWCNTs  and  the  surrounding  material  (capping  layer)  to  minimize 
thermal  impedance  mismatch  including  thermal  resistance  due  to  interfacial  roughness.  The 
capping  layer,  which  is  generally  a  metal  with  high  thermally  conductivity  and  low  yield 
strength,  aids  in  filling  the  surface  asperities  on  both  mating  solids,  thereby  reducing  the 
interface  thermal  resistance  between  the  CNT-polymer  composite  and  the  bounding  surfaces. 
Considering  diffuse  interfaces,  even  in  the  regions  where  the  tin  layer  wets  the  CNT  tips 
completely  (no  air  gaps),  phonons  from  the  CNT  tips  are  expected  to  be  readily  transferred  to  the 
Sn  layer  because  of  the  relatively  wide  phonon  spectrum  available  to  the  phonons  in  Sn  when 
compared  to  the  CNTs  (phonon  spectra  in  a  material  is  proportional  to  the  density  of  states, 
which  is  proportional  to  the  inverse  of  the  cube  of  the  sound  speed  in  the  material).  Additionally, 
the  capping  layer  is  kept  relatively  thin  to  help  minimize  its  thermal  resistance.  For  this  purpose, 
a  thin  500  ±  50  nm  layer  of  tin  (46  W/m-K  and  melting  point  of  230  °C)  is  expected  to  serve  as 
an  effective  capping  layer  material.  As  previously  mentioned,  the  top  surface  of  the  composite  is 
polished  down  to  -100  nm,  as  verified  by  AFM.  Therefore,  a  500  nm  thick  tin  layer  should  be 
sufficient  to  fill  the  asperities.  The  thin  layer  of  tin  is  deposited  on  the  composite  sample  using 
sputtering  —  the  sample  is  then  placed  in  a  vacuum  oven  at  near  melting  point  of  tin  so  that  it 
forms  a  weld  contact  with  the  exposed  CNT  tips  on  the  top  surface. 

Recently,  Prasher  et  al.  (2009)  [44]  confirmed  that  the  weld  contact  formed  by  Sn-like 
soft  (and  compliant)  metals  tend  to  increase  the  thermal  conductance  by  an  order  of  magnitude 
compared  to  more  common  Van  der  Waal  contacts  at  the  interfaces.  In  a  recent  study  by  the 
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authors, [45]  the  thermal  conductivity  of  a  500+50  nm  tin  film  was  determined  to  be  46  +  4  W/m- 
K  at  room  temperature  (see  APPENDIX  II  for  more  details).  The  thermal  conductivity  was 
found  to  remain  nearly  constant  in  the  temperature  range  of  240  K  I  300  K  (Figure  7).  This  is 
understood  to  be  due  to  the  specific  heat  being  constant  in  the  temperature  range  with  the 
Dulong-Petit  limit  being  at  210  K.  Moreover,  the  phonon  and  electron  mean  free  paths  at  room 
temperature  are  relatively  small,  50  nm  and  <  5nm,  respectively,  compared  to  the  boundary 
dimensions  [45].  Therefore,  boundary  scattering  is  not  expected  to  play  a  role  in  the  thermal 
conductivity  of  the  sample.  Furthermore,  even  though  the  reduction  in  thermal  conductivity  in 
metals  is  known  to  occur  at  elevated  temperatures  due  to  the  increased  phonon-phonon  inelastic 
scattering,  the  test  temperatures  in  our  present  case  are  not  high  enough  for  inelastic  phonon- 
phonon  (Umklapp)  interactions  to  occur,  and  thermal  conduction  continues  to  be  dominated  by 
the  specific  heat.  The  estimated  thermal  resistance  of  the  Sn  thin  film  over  the  measured 
temperature  range  is  nearly  constant  at  1.09  x  10'8  K-m2/W. 
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Figure  7:  Thermal  Conductivity  of  500  nm  Sn  thin  film  (Replotted  from  [45],  Kaul  and  Prakash,  201 1). 

3.3.  Thermal  Transport  in  Epoxy  Polymer 

The  thermal  conductivity  of  EPON-862  epoxy  is  measured  from  100  K  to  340  K  (Fig.  8), 
using  the  three-omega  technique.  The  measured  thermal  conductivity  of  EPON-862  at  room 
temperature  is  0.224+0.02  W/m-K.  The  slight  increase  in  thermal  conductivity  with  temperature 
is  understood  to  be  a  result  of  the  increase  in  specific  heat  with  test  temperature.  The  thermal 
transport  in  amorphous  polymers  (Reese  et  al.)  [46],  is  phonon-dominated  with  the  mean  free 
path  of  phonons  remaining  a  constant  due  to  scattering  by  the  amorphous  disordered  regions. 
For  the  case  of  EPON-862,  the  phonon  mean  free  path  is  estimated  to  be  0.32  nm  based  on  the 


25 


sound  velocity  (dilatational  wave  velocity)  of  1949  m/s  (corresponding  to  the  elastic  modulus 
value  of  2.58  GPa  [47]  and  Poisson©  ratio  of  0.35)  and  specific  heat  equal  to  1060  J/Kg-K  [48]. 
The  glass  transition  temperature  of  EPON-862  is  408  K  [47],  and  the  thermal  conductivity  of 
Epon-862  is  expected  to  increase  continuously  until  it  reaches  the  transition  point.  A  small  jog  in 
the  temperature  dependence  of  thermal  conductivity  is  observed  at  around  225  K,  which  is 
consistent  with  the  observation  made  earlier  in  many  polymers  for  specific  heat  dependency  with 
temperature,  and  can  be  attributed  to  the  transition  of  the  local  amorphous  regions  within  the 
polymer  from  glassy  to  the  rubbery  state.  The  rotational  and  vibrational  motions  of  atomic 
groups  associated  with  the  polymer  chains  in  the  rubbery  state  require  a  lesser  amount  of  energy 
for  reorganization  and  mobility  when  compared  to  the  glassy  state. 
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Figure  8:  Thermal  Conductivity  vs  Temperature  of  Epoxy  EPON-862  from  100  K-340  K. 
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3.4.  Thermal  transport  in  VA-MWCNTA-Epoxy  Composite 


In  order  to  characterize  the  thermal  conductivity  of  the  YA  MWCNT-epoxy  composite  the 
three-omega  method  outlined  in  Section  2.3  is  employed.  A  thin  layer  (~  500  nm)  of  tin  with  a 
Si02  insulating  layer  on  top  is  deposited  on  the  composite  substrate  to  prevent  direct  contact  with 
the  aluminum  heater  line  and  contact  pads.  The  estimation  of  both  thermal  conductivity  and 
thermal  interface  resistance  between  the  layers  is  important  to  obtain  the  through-thickness 
thermal  conductivity  of  the  VA  MWCNT-epoxy  composite.  Tong  et  al.  in  an  earlier  work  [37] 
derived  a  two-layer  model  based  on  FeldmanQ  notation.  A  three-layer  1-D  model  including  the 
finite  interface  resistance  between  layers  (based  on  FeldmanQ  notation)  has  been  derived  for  this 
work. 
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Figure  9:  Schematic  representation  of  1-D  thermal  transport  model  for  a  three- layered  structure  with  a 
heater  on  top. 
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The  main  steps  of  the  ID  thermal  model  are  provided  below.  The  heat  equation  for  a  1-D 


multilayer  problem  can  be  written  as, 


y-j 


d2Tj{y,t) 

¥ 


=  PfP, 


dTjiy,t) 

dt 


(4) 


where  j=l,  2  and  3  represent  the  individual  layers  (l=Si02,  2=Sn,  3=CNT-epoxy  composite);  7} 
is  the  temperature  at  a  point  y/  in  an  individual  layer;  and  kyj,  p.  and  Cp  represent  the  thru¬ 
thickness  thermal  conductivity,  density,  and  the  specific  heat  of  the  individual  layers. 


Under  periodic  excitation  with  angular  time  frequency,  CO  ,  the  general  solution  to  Eq.  (4)  can  be 
expressed  as 


Tj(yp  t)  =  + 


(5) 


I  ionp.C  , 

Where  rij  = - - — —  is  the  inverse  of  the  penetration  depth.  The  two  parts  of  the  solution 

V  kyJ 

represent  thermal  waves  propagating  to  the  positive  and  negative  -directions  with  two  complex 
coefficients,  A  and  B,  to  be  determined  by  boundary  conditions. 


Let  us  assume  for  this  ID  material,  a  periodic  heat  flux  with  amplitude  qo  and  frequency  co  is 
injected  at  the  surface  y1  =  0 
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q(y]  =  0,t)  =  q0e~Uo‘. 


(5) 


Following  Feldman  Q  notation, [27]  we  can  express  the  temperature  at  any  position  y  as  a  vector, 


Tj(yrt)  = 


A. e(njyj-iM) 

B, e~njyri<a,) 


(6) 


The  matrix  formulation  aids  in  finding  the  correlation  between  the  temperatures  at  two  different 
locations  within  the  medium  as  well  as  that  across  the  interface  between  the  layers.  The 
boundary  condition  is  based  on  the  assumption  that  with  no  heat  generated  or  absorbed  at  the 
interfaces,  the  heat  flux  is  continuous  across  the  interfaces,  and  there  is  discontinuity  in 
temperature  across  a  j  and  j+1  interface  due  to  a  finite  interface  thermal  resistance  R],j+i. 
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Sy3 


yj+\  =° 
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T3-^3 


where  Lj  is  the  thickness  of  the  j  layer. 


The  total  thermal  impedance  Z  of  the  three  layered  structure  neglecting  the  interface  thermal 
resistance  can  be  obtained  as 
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(8) 


where  y  =  ky  rj  and  q0  is  the  amplitude  of  the  input  heat  flux.  Note  that  the  exponential  terms 

in  the  solution  above  have  been  simplified  using  Taylor  expansion  assuming  the  films  are  very 
thin  compared  to  the  penetration  depth.  Since  the  Si02  and  Sn  films  deposited  on  top  of  the 
composite  are  very  thin  (in  nanometer  scale)  compared  to  the  penetration  depth  of  heat,  i.e. 

\r/xLx\ « 1,  |  ?]2L2\  « land  |/73L3|  ~  1  for  the  range  of  frequencies  used,  makes  this  approximation 
valid. 


Similarly,  the  total  thermal  impedance  Z  for  the  three-layered  structure  with  interface 
resistances  Ri=  R/j  (between  Si02  and  Sn)  and  AS  =  AS,.?  (between  Sn  and  MWCNT),  can  be 
expressed  as 
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The  model  is  verified  by  comparing  them  with  the  results  obtained  for  the  two  layer 
model  as  described  by  Tong  et  al.  (2006)  [37].  Moreover,  using  the  condition,  y3  =0  in  the 
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absence  of  a  third  layer,  and  \rj2L2\  ~  1  such  that,  the  heat  wave  dies  down  within  the  second  layer 

that  now  acts  as  a  substrate.  By  substituting  these  conditions  in  the  Equation  (10),  we  get  the 
same  result  as  in  Tong  et  al.  (2006),  i.e. 


„  1  „  L 

Z  =  —  +  i?j  +— - 
y2  kx 


1  -^-^Rl-yl2Rx 

7i  Yi 


(10) 


Table  1.  Model  Parameters  that  fit  the  experimental  amplitude  and  phase  of  temperature 
oscillations 
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Figure  10:  Comparison  of  Thermal  Model  with  and  without  interface  resistances  with  experimental  data: 
(a)  Temperature  oscillations  vs  frequency;  (b)  Phase  of  oscillations  vs  frequency. 
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The  experimental  data  on  temperature  oscillations  and  phase  are  compared  with  the 
predictions  of  the  thermal  model  with  and  without  thermal  interface  resistance  between  the 
layers.  The  thermal  model  is  based  on  several  parameters  (shown  in  Table  1),  e.g.,  thermal 
properties  of  the  various  layers,  and  the  unknown  thermal  interface  resistances.  An  iterative 
algorithm,  as  described  in  Tong  et  al.  (2007),  [25]  was  employed  to  obtain  the  best  fit  values  for 
the  unknowns.  The  comparison  of  the  experimental  results  and  model  predictions  based  on  the 
best  fit  values  for  the  unknown  parameters  is  shown  in  Fig  10.  From  these  comparisons  the 
thermal  conductivity  of  the  CNT-Epoxy  composite  is  estimated  to  be  5.8  ±  0.64  W/m-K  at  room 
temperature.  The  measured  amplitude  as  well  as  the  phase  of  the  temperature  oscillations  are 
found  to  be  sensitive  to  the  thermal  conductivities  of  tin,  MWCNT-epoxy  composite,  and  the 
thermal  interface  resistances.  A  ±20%  uncertainty  in  thermal  conductivity  of  either  tin  or  the 
MWCNT-epoxy  composite,  results  in  a  corresponding  ±  5-10%  shift  in  the  amplitude  of  the 
temperature  oscillations  from  the  measured  values.  The  thermal  conductivity  of  the  composite 
exhibits  a  slight  increase  with  temperature  between  240  K  and  300  K  (Fig.  11),  and  is  understood 
to  be  dominated  by  the  specific  heat  dependence  with  temperature  of  the  Epon-862  epoxy  [9, 
49], 
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Figure  11:  Thermal  Conductivity  Vs  Temperature  of  VA-MWCNTA-Epoxy  (EPON-862) 

Nanocomposite  from  240K-300K. 

The  best  fit  estimates  for  the  interface  thermal  resistances  Rj.2  (between  SiC>2  and  Sn)  and 
R2-3  (between  Sn  and  MWCNT)  are  5  x  10'5  m2-K/W  and  8  x  10'6  T  8.5  x  10'7  m2-K/W, 
respectively.  The  amplitude  or  phase  of  temperature  oscillations  is  especially  sensitive  to  the 
thermal  resistance  of  the  rougher  interface  (i.e.,  the  SiC^T  Sn  interface)  in  our  case.  However,  it 
is  relatively  less  sensitive  to  the  interface  thermal  resistance  of  the  relatively  smooth  Sn- 
MWCNT  interface,  and  therefore  by  keeping  all  other  parameters  constant  we  can  estimate  the 
interface  resistance  of  the  smooth  Sn-MWCNT  interface  to  within  an  order  of  magnitude  range. 
The  interface  resistance  between  the  layers  is  a  combination  of  Kapitza  resistance  based  on 
diffuse  mismatch  model  (DMM)  and  constriction  resistance  at  the  contacts  of  carbon  nanotubes 
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with  the  Sn  layer.  Recently,  Prasher  et  al.  (2005,  2008)  [50,  51],  suggested  that  for  nano-sized 
constrictions  (which  the  vertically  aligned  MWCNTs  make  with  the  layer  above  it),  the  effects  of 
mismatch  in  acoustic  properties  are  more  dominant  than  due  to  constriction  of  the  heat  flux  lines. 
Therefore,  vibrational  spectra  mismatch  between  Sn,  a  soft  metal  with  relatively  low  sound 
speed,  and  VA  MWCNTs  is  expected  to  restrict  the  thermal  conductance  through  the  MWCNT 
composite  when  bounded  by  surfaces  on  both  sides.  The  Debye  temperature  of  Sn  is  170  K, 
while  the  Debye  temperature  for  multi-walled  carbon  nanotubes  is  much  higher,  960  K-2500  K, 
indicating  a  large  mismatch  in  vibrational  spectra  of  the  two  materials.  Interface  thermal 
resistance  between  two  highly  dissimilar  materials  with  large  differences  in  phonon  spectra,  e.g  , 
interfaces  involving  metal  and  dielectric  solids,  as  pointed  out  by  Lyeo  and  Cahill  (2006)  [52]  is 
not  only  due  to  the  coupling  of  electrons  in  a  metal  and  phonons  in  a  dielectric  substrate,  but  also 
by  anharmonic  processes,  namely  three-phonon  processes  that  can  contribute  a  significant 
additional  channel  for  the  transport  of  heat  by  altering  phonon  frequency  by  inelastic  scattering 
processes.  For  these  reasons,  carbon  nanotubes  may  be  particularly  attractive  for  exchange  of 
thermal  energy  with  Sn  across  the  CNT-Sn  interfaces  involving  both  phonons  and  electrons. 

Interestingly,  even  after  a  spate  of  research  on  VA  MWCNT  composites,  there  are  only 
few  published  works  in  the  literature  that  reported  thermal  conductivities  greater  than  1  W/m-K 
(Figure  12).  Sihn  et  al.  (2008)  [26]  have  reported  the  highest  thermal  conductivity  so  far  for  a 
30pm  VA  MWCNT-epoxy  composite.  Using  a  simple  1-D  thermal  resistance  network  model  to 
estimate  the  thermal  conductivity  of  their  MWCNT  composite  from  their  reported  data,  we 
estimate  the  room  temperature  thermal  conductivity  of  their  VA  MWCNT-epoxy  composite  to  be 
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137  W/m-K.  Considering  a  10  vol.  %  of  MWCNTs  in  their  composite,  the  thermal  conductivity 


of  individual  MWCNTs  is  in  excess  of  1000  W/m-K. 
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Figure  12:  Thermal  Conductivity  vs  Volume  Fraction  of  Nanowires  and  Nanotubes. 

The  primary  reasons  for  a  less  than  1  W/m-K  thermal  conductivity  of  the  MWCNT 
composites  reported  in  the  literature  may  be  due  to  the  low  volume  fraction  of  CNTs  used  in  the 
composites  as  well  as  their  defect  state  and  purity.  Additionally,  in  order  to  improve  their 
performance  as  TIM  material,  thermal  conductance  at  CNT-capping  layer  interfaces  need  to  be 
addressed.  However,  results  of  the  present  study  suggest  that  the  thermal  interface  resistance  for 
relatively  smooth  and  polished  surfaces  may  contribute  to  only  second  order  effects  in  the 
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thermal  performance  of  the  TIMs  at  room  temperature.  In  a  parallel  study  [40],  the  authors 
investigated  thermal  conductivity  in  annealed  and  as  received  thermal  CVD  grown  individual 
multiwalled  carbon  nanotubes.  They  reported  a  five-fold  increase  in  thermal  conductivity  of  the 
MWCNTs  from  the  annealed  batch.  Moreover,  it  was  also  noted  that  depending  on  the  degree  of 
inherent  physical  deformities/defects  in  the  CNTs,  thermal  conductivity  of  even  the  annealed 
samples  could  be  lower  than  those  obtained  for  the  as  received  batch  (Figure  13).  High 
resolution  SEM  micrographs  of  CNT  arrays  (Figure  14)  from  the  same  batch  as  the  array  that 
was  used  in  the  fabrication  of  the  composite  indicate  physical  deformities  and  defects,  including 
entanglement  of  carbon  nanotubes  and  fused  contacts,  that  are  understood  to  be  the  reason  for 
reduced  thermal  conductivity  of  the  composites  in  general.  It  is  to  be  noted  that  most  previous 
works  reported  in  the  literature  on  VA  MWCNT-polymer  composites  [15,  17,  26,  53]  show 
room-temperature  thermal  conductivity  as  high  as  metal  nanowire  array-polymer  composites  [54, 
55]  or  the  semiconductor  nanowire  array-polymer  composites  [56,  57]  with  nanowire  volume 
fractions  in  the  same  range  as  the  carbon  nanotubes.  Nevertheless,  the  five-fold  increase  in 
thermal  conductivity  in  the  post-annealed  MWCNTs  suggests  that  the  thermal  conductivity  of  VA 
MWCNT-epoxy  composites  can  be  designed  to  be  as  high  as  25  W/m-K,  which  can  potentially 
have  a  major  impact  in  the  design  of  multifunctional,  load-bearing,  high-thermal  performance 
structural  interfaces  for  a  variety  of  thermal  management  applications. 
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Figure  13:  Thermal  Conductivity  vs  Raman  Ratio  of  individual  multi- walled  nanotubes.  [Replotted  from 
Ref  [40],  Bifano,  Kaul  and  Prakash,  2011] 


Figure  14:  SEM  micrographs  (9  kV,  3.5  mm  working  distance,  and  spot  size,  2.5)  of  MWCNT  array 
batch  of  which  the  composite  was  made  and  tested:  (a)  Entanglements  of  Individual  Nanotubes  (87,000X) 
and  (b)  Individual  samples  are  often  attached  as  if  they  are  fused  together  (100,000X). 
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SUMMARY 


This  paper  reports  development  of  vertically  aligned  multi-wall  carbon  nanotube  (VA- 
MWCNT)  array  composites  for  thermal  energy  management  in  load-bearing  structural 
applications.  Unlike  previous  studies  on  the  characterization  of  VA-MWCNT-based  thermal 
interface  materials  (TIMs)  for  primarily  non-load  bearing  applications,  the  material  systems  of 
interest  here  involve  the  use  of  VA  MWCNTs  in  an  epoxy  matrix.  The  epoxy  matrix  imparts 
mechanical  strength  to  these  systems  while  the  VACNTs  provide  avenues  for  high  thru-thickness 
thermal  conductivity  across  a  typical  material  interface.  In  order  to  obtain  the  thermal 
characteristics  of  these  multifunctional  TIMs,  we  report  measurements  of  thermal  conductivity  in 
the  Sn-capped  VA-MWCNT-epoxy  composites  as  well  as  in  its  individual  constituents,  i.e.,  bulk 
EPON-862  (matrix  material)  and  Sn  thin  film,  in  the  temperature  range  240  K  T  300  K,  and 
individual  multi-wall  carbon  nanotubes  at  room  temperature  taken  from  the  same  VA-MWCNT 
batch  as  the  one  used  to  fabricate  the  CNT-epoxy  TIM.  The  thermal  conductivity  of  the  epoxy 
and  Sn  thin  film  was  obtained  as  a  function  of  temperature  by  using  a  cryostat  in  conjunction 
with  the  three  omega  method.  The  thermal  conductivity  of  individual  free  standing  MWCNT 
samples  was  obtained  by  employing  the  Wollaston  T-type  three-omega  probe  method  inside  a 
high  resolution  scanning  electron  microscope  equipped  with  nanomanipulators  and  a  gas 
injection  system  for  electron  beam  induced  platinum  deposition.  The  thermal  conductivity  of 
bulk  EPON-862  epoxy  was  observed  to  increase  gradually  from  0.1  W/m-K  at  100  K,  to  0.24 
W/m-K  at  340  K,  and  is  understood  to  be  dominated  by  the  heat  capacity  of  the  epoxy.  The 
thermal  conductivity  of  Sn  capping  layer  (500  micron  film  thickness)  was  observed  to  remain 
nearly  constant  at  about  46  W/m-K  over  the  above  mentioned  temperature  range.  The  thermal 
conductivity  of  individual  freestanding  carbon  nanotubes  was  measured  to  be  about  60  W/m-K  at 
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room  temperature.  A  1-D  multilayer  thermal  model  that  includes  effects  of  thermal  interface 
resistance  and  the  thermal  conductivity  of  the  CNT-epoxy  composite  and  its  constituents  was 
used  to  interpret  the  experimental  results.  The  thermal  conductivity  of  the  carbon  nanotube- 
epoxy  composite  was  estimated  to  be  about  5.8  W/m-K  and  exhibits  a  slight  increase  with 
temperature  in  the  range  of  240  K  to  300K.  The  best  fit  estimates  for  the  interface  thermal 
resistances  between  Si02  and  Sn  and  between  Sn  and  MWCNT  were  5  x  10'5  m2-K/W  and  8  x 
10'  I  8.5  x  10'  m  -K/W,  respectively.  The  results  of  this  study  suggests  that  the  inclusion  of  a  Sn 
thin  layer  on  the  VA-MWCNT  array,  as  well  as  the  morphological  structure  of  the  individual 
MWCNT&  are  dominating  factors  that  control  the  overall  thermal  conductivity  of  the  TIM.  These 
results  are  encouraging  in  light  of  the  fact  that  the  thermal  conductivity  of  a  VA-MWCNT  array 
can  be  increased  by  an  order  of  magnitude  by  using  a  standard  high-temperature  post- annealing 
step.  In  this  way,  multifunctional  (load  bearing)  TIMs  with  effective  through  thickness  thermal 
conductivities  as  high  as  25  W/m-K,  can  be  potentially  fabricated. 
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APPENDIX  I 


E  ect  o  Heat  Treatment  on  T  erma  Con  uct  t  o  In  ua 
Mu  t  -  a  e  Carbon  Nanotube 


ABSTRACT 

Thermal  conductivity  measurements  in  commercially  available  chemical  vapor 
deposition  (CVD)  grown  heat-treated  and  non-heat  treated  multiwalled  carbon  nanotubes 
(MWCNTs)  are  reported.  The  thermal  conductivity  of  individual  samples  is  measured  using  a 
suspended  platinum  (Pt)  wire  as  a  thermal  resistance  probe  in  a  “T-type”  configuration.  Changes 
in  third  harmonic  voltage  are  measured  across  the  heated  suspended  platinum  wire  as  a  specimen 
is  attached  to  the  Pt  wire’s  mid-point.  An  analytic  model  is  used  to  correlate  the  reduction  in  the 
average  temperature  of  the  probe  wire  to  the  thermal  resistance  (and  thermal  conductivity)  of  the 
attached  sample.  Experiments  are  implemented  inside  a  scanning  electron  microscope  equipped 
with  nanomanipulators  and  electron  beam  platinum  deposition  to  improve  thermal  contact 
resistances.  The  results  indicate  a  nearly  5  fold  increase  in  the  average  thermal  conductivity  of 
MWCNT  samples  annealed  with  a  20  hour  3000°C  annealing  heat  treatment  compared  to  the  as- 
grown  samples.  However,  specimen  specific  morphological  defects  are  found  to  negate,  to  a 
large  degree,  the  advantage  of  the  heat  treatment  process.  The  average  thermal  conductivity  of 
the  heat-treated  MWCNT  specimens  was  measured  to  be  216  W/m-K,  with  the  highest  measured 
thermal  conductivity  being  730  ±153  W/m-K.  The  results  highlight  the  importance  of  sample 
quality  of  MWCNTs  for  use  in  thermal  management  applications. 


1.  INTRODUCTION 


There  is  little  agreement  amongst  the  few  reported  thermal  conductivity  measurements  [1-5] 
in  MWCNTs.  This  likely  due  to  the  vast  structural  and  chemical  differences  between  the  samples 
measured.  Room  temperature  thermal  conductivity  measurements  in  MWCNTs  have  been  reported 
to  be  as  high  as  3000  W/m-K  [1]  and  as  low  as  17  W/m-K  [5].  The  chemical  vapor  deposition 
(CVD)  process  commonly  used  in  the  synthesis  of  SWCNT  and  MWCNTs  is  known  to  produce  a 
variety  of  CNT  qualities  depending  on  growth  temperature,  catalyst  and  carbon  source  gasses  used 
[6], 


A  2001  study  [7]  showed  many  structural  and  chemical  improvements  in  CVD  grown 
MWCNTs  with  high  temperature  annealing.  Annealing  temperatures  that  exceeded  2250°C  were 
found  to  produce  x-ray  diffraction  patterns  similar  to  those  of  highly  graphitic  MWCNTs  grown  by 
arc-discharge  methods.  Smaller  interlayer  wall  spacing,  indicative  of  highly  graphitic  structure, 
were  also  possible  with  higher  annealing  temperatures  up  to  3000°C.  Residual  Fe,  the  primary 
chemical  vapor  deposition  growth  catalyst  for  MWCNTs,  was  also  found  to  reduce  from  7.1  %  by 
wt.  after  CVD  growth  to  less  than  0.01  %.  Major  Fe  reduction  occurred  both  at  the  ends  and  in  the 
core  of  the  MWCNTs  once  the  annealing  temperature  surpassed  the  vaporization  temperature  of  Fe 
(>1800°C).  While  heat  treatment  improved  the  microstructural  quality  of  the  samples,  larger  defects 
that  created  bends  and  branching  in  the  tubes  were  not  removed  in  the  heat-treatment  process.  More 
continuous  crystalline  planes  and  a  reduction  in  residual  catalyst  particles  are  expected  to  increase 
phonon  mean  free  path  (MFP),  and  thus  increase  the  lattice  (phonon)  thermal  conductivity  in 
MWCNTs. 
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More  recently,  Jin  et  al.  [8]  reported  that  annealing  of  macroscopic  MWCNT  bundles  at 
2800°C  for  4  hours  resulted  in  a  9  fold  increase  in  thermal  conductivity  at  300  K  (~2.5  to  23  W/m- 
K).  However,  from  these  measurements  an  improvement  in  intrinsic  thermal  conductivity  of 
individual  MWCNT  could  not  be  unambiguously  ascertained  due  to  the  possibly  dominating  effects 
of  thermally  resistive  junctions  between  tubes,  which  may  have  also  been  affected  by  the  heat- 
treatment  process.  To  the  best  of  the  authors’  knowledge,  heat-treatment  effects  on  the  thermal 
conductivity  of  individual  MWCNTs  have  not  been  characterized. 

Over  the  past  50  years,  thermal  properties  of  electrically  conductive  samples  and  their 
surroundings  have  been  investigated  by  measuring  the  3  harmonic  voltage  signal  produced  when 
the  sample  experiences  joule  heating  via  an  electric  alternating  current  [9-11].  Comparison  of  the 
measured  3rd  harmonic  voltage  across  the  sample  with  an  analytic  model  permits  the  thermal 
properties  of  the  sample  to  be  determined.  The  most  simplistic  steady-state  thermal  conductivity 
measurement  method  is  the  “self-heating”  technique  [12].  In  general,  the  self-heating  method 
requires  samples  to  be  electrically  conductive,  transport  heat  diffusively,  and  the  temperature 
coefficient  of  resistance  of  the  sample  to  be  known.  The  self-heating  method  was  first  applied 
to  a  bundle  of  CNTs  [12],  and  subsequently  to  individual  MWCNTs  [3]. 

An  extension  of  the  self-heating  technique  is  the  “T-Type  Probe  method,”  which  replaces  the 
sample  with  a  pre-calibrated  suspended  wire  of  known  electrical  resistivity  and  temperature 
coefficient  of  resistance  (Fig.  1).  In  this  way,  by  attaching  a  free  standing  specimen  to  the  probe 
wire  and  measuring  the  reduction  in  average  temperature  (inferred  from  the  change  in  voltage)  of 
the  suspended  probe  wire,  the  thermal  resistance,  and  thus,  thermal  conductivity  of  the  specimen 
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is  determined.  To  date,  the  T-type  method  has  been  used  to  measure  a  variety  of  microscale  samples 
[13,  14].  In  2005  Fujii  et  al.  [2]  developed  a  micro-fabricated  Pt  thermal  resistance  T-type  probe  that 
used  3co  voltage  detection  to  measure  the  thermal  conductivity  of  individual  MWNCTs. 

In  this  study,  a  T-Type  probe  composed  of  Wollaston  wire  is  implemented.  The  use  of 
commercially  available  Wollaston  wire  was  first  introduced  in  2007  in  both  lco  and  3co 
configurations  for  T-Type  CNT  thermal  conductivity  measurements  of  individual  MWCNTs  at  room 
temperature  [5],  However,  only  one  measurement  of  17  W/m-K  was  reported.  The  low  value  of  the 
thermal  conductivity  was  attributed  to  high  thermal  contact  resistance  at  the  probe-sample  interface 
and/or  the  possibility  of  ballistic  transport  in  the  rather  short  length  sample  (380  nm)  tested. 

Many  experimental  techniques  to  measure  thermal  conductivity  in  individual  nanoscale 
samples  currently  employ  a  micro-fabricated  sensing  device  [1,  2,  15,  16].  Use  of  commercially 
available  Wollaston  wire  has  the  advantage  of  being  extremely  cost  effective  compared  to 
conventional  microfabrication  methods,  and  allows  for  a  large  volume  of  samples  to  be 
characterized  in  a  short  span  of  time.  In  the  present  investigation,  all  thermal  property 
characterization  studies  are  conducted  inside  a  high-resolution  scanning  electron  microscope 
equipped  with  nanomanipulators  and  gas  injection  system  for  high  throughput  analysis.  Moreover, 
appropriate  choice  of  materials  and  package  design  allows  for  multiple  Wollaston  T-Type  probes  to 
be  prepared  with  varying  thermal  resistances,  and  thus  measurement  sensitivities. 

2.  METHODS 

2.1.  Three  Omega  Measurements  for  the  T-Type  Probe  Method 
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In  order  to  conduct  the  three  omega  measurements,  the  probe  wire  is  heated  by  a  low 
frequency  alternating  current  I(t)=Iiwsin(cot),  where  Ii0J  is  the  current  amplitude.  The  electrical 
resistance  of  the  probe  wire  at  zero  current  is  Reo.  For  sufficiently  low  heating  currents  I(t), 
Joule  heating  is  assumed  to  be  of  the  form, 

Q(t)  =  I2(t)Reo  =  ^  /  lo  Reo  (1  -  cos(  2ojt))  (1) 

While  conducting  the  experiments,  attention  must  be  paid  to  maintain  the  validity  of  Eq.  (1).  The 
temperature  rise  of  the  probe  wire  must  be  limited  to  only  a  few  degrees  above  ambient  temperature. 
If  care  is  not  taken,  the  average  heat  input  would  be  better  modeled  using  an  iterative  Q(t)  =  I  (t)R(t) 
instead  of  Q(t)  =  I2(t)Reo,  resulting  in  a  non-linear  dependence  of  the  change  in  probe  wire  electrical 
resistance  to  Joule  heating.  The  non-linearity  is  carefully  monitored  through  comparison  of  linear 
and  quadratic  regression  models. 

The  spatially  averaged  temperature  response  of  the  wire  6(t)  is  proportional  to  the  joule 
heating  by  a  thermal  transfer  function  Z0  as, 

S(t)  =  Z0Q(t)  (2) 


As  long  as  the  heating  current  is  driven  at  a  sufficiently  low  frequency  the  thermal  response  of  the 
probe  wire  remains  in  phase  with  the  heating  current  and  the  thermal  transfer  function  is  largely 
frequency  independent  (See  Appendix  I).  Since  Z0  is  a  function  of  the  thermal  resistance  of  the 
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probe  or  the  probe-sample  combination,  experimental  determination  of  ZQ  allows  the  sample’s 
thermal  resistance  and  thus  thermal  conductivity  to  be  obtained. 

Using  Eq.  (2)  along  with  the  electrical  resistance  response  of  the  probe  wire, 
Re  (t)  =  Reo  [1  +  a  0  (t)] ,  where  a  is  the  temperature  coefficient  of  resistance,  the  electrical  resistance  of 
the  probe  wire  Re(t)  will  experience  oscillations  at  2co  given  by 

Re  (0  =  Reo+\  ccZJ^Rj  - i  aZJ^Rl,  cos(2 cot)  .  (3) 

Using  Ohm’s  law  along  with  Eq.  (3),  the  voltage  across  the  probe  wire  can  be  written  as 

=  KK+\ocZ0Rlll  sin«a-  Uz/,!X 

Equation  4  shows  that  the  measured  voltage  is  the  combination  of  the  Ico  and  3co  frequency 
harmonics  of  the  initial  heating  current.  A  Lock-in  Amplifier  is  used  to  selectively  measure  the  RMS 
voltage  corresponding  to  either  component.  The  reported  MWCNT  data  in  this  study  utilizes  the 
3a)  detection  since  the  constant  ohmic  signal  contained  within  the  lco  signal  is  expected  to  lead 
to  a  greater  degree  of  scatter.  In  terms  of  the  RMS  voltage,  the  voltage  at  the  third  frequency 
harmonic,  V3(0  RMS ,  can  be  written  as 
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(5) 


^3w,RMS  ~  2  a^ohm,RMS^eoQRMS 


2 

where  Qrms  =  I  ho .rms  Reo-  Dividing  both  sides  of  Eq.  (5)  by  Ii<o,rms  ,  the  3co  RMS  electrical 
resistance,  defined  as  Re3(0:  rms  =  V3c0i  rms/I im.  rms,  is  written  as, 


Reioj.RMS  2  a^eo^oQ-RMS  '  (^) 

Next,  the  thermal  transfer  function  Z0  is  experimentally  determined  by  obtaining  the  slope  of  the 
Re3w,RMs  vs.  Qrms  plot.  Since  Za  is  a  function  of  the  sample  thermal  resistance,  Rti,  s,  and  the  probe- 
wire  thermal  resistance,  Rti,,p,  comparing  Z0  with  model  predictions  (see  Sec.  2.2)  for  joule  heating  in 
the  probe  wire  with  and  without  the  sample  allows  both  Rt/,  p  and  Rth,s  to  be  determined. 

2.2.  Analytical  Model  for  the  Thermal  Transfer  Function  Z0 

The  steady  state  thermal  response  of  the  probe  wire  is  modeled  using 


/  d26{x)  _  _  Qrms 
P  dx 2  ApLp  ’ 


(7) 


where  x  is  the  coordinate  along  the  length  of  the  probe  wire  (Fig.  1),  kp  is  the  thermal  conductivity  of 
the  probe  wire,  LP  is  the  length  of  the  probe  wire,  and  6(x)  is  the  spatially  dependent  temperature 
rise.  Convection  and  radiation  heat  losses  are  neglected  by  performing  the  measurement  with  very 
small  heating  amplitudes  and  under  vacuum  conditions  within  the  SEM  chamber  [17]. 
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The  thermal  resistance  of  the  sample,  Rth,s,  is  incorporated  into  the  model  by  a  flux  boundary 


condition  at  the  point  of  sample  attachment,  i.e., 


m  _k  A 

d6(x) 

_d6{x) 

\ 

-  H p/ip 

, th,S 

dx 

V 

x=0+  SX 

X=0~  y 

(8) 


The  boundary  conditions  at  the  ends  of  the  probe  wire  are 

6(+Lp/  2)  =  6{-Lp  /  2)  =  0 


(9) 


The  piecewise  parabolic  solution  for  positive  and  negative  x  satisfying  Eqs.  (7),  (8),  and  (9)  is 


0(x,?))-QrmsLp 


8  kPSp 


1- 


'  x 


L!  2 


+ 


'jC 

1  +  r/ 


LI  2 


-1 


(10) 


where  the  parameter  //  is  the  ratio  of  the  thermal  resistance  of  the  probe  wire  to  the  thermal 
resistance  of  the  sample  R,i,,s  and  is  defined  by  //  =  Rth,p/4Rth,s  [5].  The  thermal  resistance  of  the 
probe  wire  and  the  sample  is  given  by  Ru,,p  =  Lp/kpAp  and  Rth,s  =  Ls/ksAs  ,  respectively.  In  the 
absence  of  the  sample  Rth.s  =  oo,  rj  =  0,  and  therefore  the  well-known  inverted  parabolic  solution 
for  a  Joule  heated  suspended  wire  is  recovered. 
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Integrating  over  the  length  of  the  probe  wire,  the  spatially  averaged  temperature  rise  0  over 


the  length  of  the  sample  can  be  written  as 


s=v?‘ 


U 

RMSth,P 


l-|(l  +  77  l)  1 


(11) 


Comparison  of  Eq.  (11)  and  (2)  reveals  the  theoretical  thermal  transfer  function  to  be, 


z  —  1  /? 

°  12  th>p 


i-4(i +rlrl 


(12) 


2.3.  Experimental  Procedure 

Inserting  Eq.  (12)  into  Eq.  (6)  and  taking  rj  =  0  (absence  of  sample),  the  thermal  resistance  of 
the  probe  wire  is  first  estimated  using 


24 

RfhJ>  =  ~aR~, 


AR 


e3co,RMS 


RMS 


J  No  Sample 


(13) 


As  aside,  replacing  the  probe  wire  with  a  sample  of  unknown  thermal  resistance  permits  the 
thermal  resistance  and  therefore  thermal  conductivity  of  the  sample  to  be  determined  using  Eq.  (13). 
The  caveat  is  that  the  sample’s  resistance  and  temperature  coefficient  of  resistance  must  be  known. 
This  method  is  termed  as  the  “Self-Heating  Method”. 
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Next,  the  sample  is  attached  and  the  new  slope  of  Resm,ms  vs.  Qrms  is  measured.  By 


comparison  of  the  ratio  of  the  slopes,  <p,  where  cp  is  defined  as  [5], 


* 


RMS  )with  Sampie 
(A  R  eilB, rms  /AQrmS  )No  Sample 


(14) 


The  thermal  resistance  of  the  sample  is  found  using, 


D  _  J_  D 

^th,S  4  ^th,P 


f i  ^ 

4 


(15) 


where  s  =  \  —  <f>. 


The  sample’s  thermal  conductivity  is  determined  from  its  thermal  resistance  by  ks=Ls/Rth,sAs , 

2  2 

where  the  sample’s  cross-sectional  area  is  As  =  n(R0  -Ri).  In  the  present  study,  the  inner  radius  is 
approximated  using  the  ratio  of  the  average  inner  to  outer  diameter  supplied  per  batch  by  both 
samples  manufacturers. 
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2.4.  Heater/Sensor  Device  for  in-situ  SEM  testing 


The  heater/sensor  probe  wire  device  used  in  the  study  is  fabricated  from  commercially 
available  Wollaston  wire  (Sigmund-Cohn  Corp.,  Mt.  Vernon,  NY);  the  wire  is  composed  of  a  75pm 
silver  sheath  surrounding  an  approximate  500-750  nm  platinum  core.  A  concentrated  nitric  acid 
solution  is  used  to  selectively  etch  the  silver  sheath  exposing  the  Pt  core  (Fig.  2a). 

To  increase  throughput,  the  entire  measurement  system  is  composed  of  four  etched 
Wollaston  wires,  numbered  1  to  4  in  Fig.  2a.  While  only  one  probe  wire  is  utilized  per 
experiment,  due  to  the  extreme  delicacy  of  the  fine  platinum  wires  a  custom  switch-box  is  used 
to  independently  activate  each  wire  (not  shown).  The  device  comprises  four  probe  wires,  with 
nearly  the  same  average  diameters  but  with  varying  lengths,  and  thus  with  different  thermal 
resistances.  Figure  2b  shows  the  four  probe  wires  mounted  inside  the  SEM  chamber  (FEI  FE 
NanoSEM  600).  All  4  probe  wires  are  reachable  by  nanomanipulators  inside  the  SEM  chamber. 

2.5.  Probe  Wire  Thermal  Resistance  and  Sensitivity 

The  length  of  the  exposed  platinum  core,  and  thus  the  thermal  resistance  of  the  probe  is 
readily  controlled  by  carefully  maneuvering  the  droplet  of  nitric  acid  dining  the  etch  process.  This 
permits  the  sensitivity  of  the  probe  to  be  maximized  relative  to  the  sample  being  tested. 

If  the  relative  uncertainty  in  the  measured  temperature  is  5016  (where  56  is 
approximately  0.05K  (as  measured  from  the  linear  regression  of  Re3m,RMs  vs.  Qrms  using  the 
experimentally  determined  temperature  coefficient  of  resistance  of  the  probe  wire  a  ~  0.0033,  and 

maintaining  the  average  temperature  rise  of  the  probe  wire,  6  ~  10  K),  then  the  ratio  of  5616  to 
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the  maximum  relative  uncertainty  in  the  sample  resistance,  dRth,s/Rth,s,  is  a  quantifiable  sensitivity 


factor  [5] 


sei_e 

3Rth,S  /  R„,s 


(16) 


If  6Rth,s/Rth,s  =0.10,  Eq.  (16)  is  readily  solvable  for  maximum  and  minimum  values  of  //.  While 
maximum  probe  sensitivity  occurs  when  the  probe  and  sample’s  thermal  resistance  is  rj  =  1,  the 
uncertainty  in  the  measured  thermal  resistance  of  the  sample  will  be  within  10  %  of  its  true  value 
if  the  actual  value  of  rj  is  between  0.077  and  12.923.  The  average  thermal  resistance  of  the  all 
prepared  Wollaston  wire  probes  used  in  this  study  is  4.02  x  107  ±  1.15  x  107  W/K.  Therefore, 

o  n 

samples  having  thermal  resistances  as  high  as  5.20  x  10  W/K  and  as  low  as  3.11  x  10  can  be 
effectively  measured. 


2.6.  Method  Validation 

The  analysis  of  the  experiment  and  equipment  setup  is  verified  by  both  a  lco  and  3co  self¬ 
heating  experiment  on  99.99%  Au  wire.  The  Au  wire  is  chosen  for  its  uniformity  of  diameter 
(3.4%  variation)  over  the  etched  Wollaston  probe  wire  that  has  an  average  of  a  17%  variation  in 
diameter.  The  lco  and  3co  Au  self-heating  validation  experiments  result  in  a  2.3%  and  3.6% 
difference  from  literature  value  (306.6  W/m-K),  respectively.  Because  the  Wollaston  T-type 
probe  method  measures  thermal  resistance,  diameter  variation  leads  to  twice  the  uncertainty  in 
cross-sectional  area  resulting  in  a  variation  in  thermal  conductivity  of  no  less  than  34%. 
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Therefore,  comparison  of  the  measured  thermal  conductivity  of  the  Wollaston  probe  wire  with 
that  of  bulk  platinum  is  infeasible. 

The  use  of  the  suspended  Pt  wire  as  a  thermal  resistance  sensor  is  also  validated  (Fig.  3). 
A  separate  length  of  etched  Wollaston  wire  is  soldered  to  a  nanomanipulator  and  treated  as  a 
sample  of  unknown  thermal  conductivity.  To  maximize  the  sensitivity  of  the  experiment,  the 
sample  length  is  chosen  such  that  r\~  1 . 

Although  the  lengths  of  etched  Wollaston  wire  are  sequentially  cut  from  the  same  batch, 
sample  and  probe  diameters  both  exhibit  variations  in  diameter.  Instead  of  comparing  thermal 
conductivity  of  the  probe  wire  to  the  sample,  the  thermal  resistance  per  unit  length  of  the  Pt 
probe  wire  is  compared  to  that  of  the  thermal  resistance  per  unit  length  of  the  sample  Pt  wire. 
lco  and  3co  experiments  measure  the  Pt  sample’s  thermal  resistance  per  unit  length  to  be  within 
5.7%  and  4.1%,  respectively,  of  the  Pt  probe.  Pt-EBID  is  deposited  at  the  contact  point  to 
increase  the  contact  area  by  more  than  75%  and  thus  similarly  reducing  thermal  contact 
resistance  (Fig.  3a).  Figure  3b  depicts  the  change  in  slope  of  Re3m>RMs  vs.  Qrms  caused  by  the 
reduction  in  the  average  temperature  of  the  Wollaston  probe  wire  due  to  contacting  the  sample,  and 
the  subsequent  reduction  in  the  probe  wire’s  average  temperature  due  to  the  reinforcement  of  the 
sample-probe  interface  with  Pt-EBID.  Prior  to  applying  the  Pt-EBID,  the  measured  sample’s 
thermal  resistance  is  2.9  times  greater,  which  emphasizes  the  importance  of  improving  the 
thermal  contacts,  especially  as  sample  diameters  decrease. 

3.  MWCNT  Samples 
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Two  sample  groups  are  investigated  in  the  present  study:  1)  MWCNT  samples 
synthesized  via  thermal  CYD;  and  2)  MWCNT  samples  synthesized  via  thermal  CYD  and 
subsequently  heat-treated  at  3000°C  for  20  hours  (thermal  annealed).  The  heat-treated  thermal 
CVD  MWCNT  samples  are  acquired  from  US  Nanomaterials  Research  (US4400,  US4406, 
US4412).  The  second  set  of  samples,  termed  “non-heat  treated”,  are  comprised  of  thermal  CVD 
grown  MWCNTs  purchased  from  US  Nanomaterials  Research  (US4309,  US4315)  and  NanoLab 
(PD15L5-20).  The  NanoLab  samples  were  synthesized  at  600  ±  15°C  and  all  US  Nanomaterials 
Research  samples  were  grown  at  665  ±  15°C. 

The  samples  tested  have  a  high  degree  of  variation  in  diameters,  shapes,  lengths,  and 
straightness.  Comparison  of  Fig.  4a  to  4c  and/or  4e  indicates  that  the  high  temperature 
graphitization  process  shortens  the  lengths  of  the  MWCNT  samples  by  more  than  75%. 
Adjacent  nanotubes  appear  to  become  ‘seared’  together  (Fig.  4b).  A  greater  majority  of  heat- 
treated  samples  have  debris  attached  to  the  outside  sample  walls  when  compared  to  the  non-heat 
treated  MWCNTs  (Fig.  4f).  The  non-heat  treated  batches  (US4309,  US4315)  have  a  higher 
degree  of  abrupt  changes  in  angle  (Fig.  4c  and  d)  compared  to  the  PD15L5-20  batch  (Fig.  4e  and 
f).  The  source  of  these  angle  changes  are  likely  due  to  defects  in  the  lattice  structure  of  the 
CNTs  [7], 

3.1.  Raman  Spectroscopy 

Raman  Spectroscopy  (excitation  wavelength  of  514  nm)  is  used  to  validate  the 
improvement  in  graphitization  and  purity  for  the  heat-treated  sample  batches  compared  to  the 
non-heat  treated  sample  batches.  An  assessment  of  both  residual  amorphous  carbon  and  defects 
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is  made  by  examining  the  ratio  of  the  D  band  peak  (-1350  cm'1)  intensity  to  the  G  band  peak 
(-1580  cm"1)  intensity  [18].  The  D  band  (defect  band)  is  activated  by  both  the  carbonaceous 
impurities  with  sp3  bonding  and  fragmented  sp2  bonds.  These  features  are  indicative  of 
MWCNT  defects  [19].  The  G  band  is  activated  by  contiguous  sp2  bonding,  i.e.,  a  high  degree  of 
crystallinity.  Therefore,  when  comparing  two  sample  groups,  a  smaller  DIG  band  ratio  is 
representative  of  a  sample  with  fewer  defects  and  a  higher  degree  of  graphitization.  To  date, 
standards  to  quantify  purities,  diameters,  wall  numbers,  and  most  importantly  growth  conditions 
of  MWCNT  samples  are  not  well  established.  Lacking  a  standardized  MWCNT  sample  set, 
Raman  spectroscopy  remains  a  qualitative,  yet  effective,  assessment  of  purity.  It  is  important  to 
note  that  the  Raman  laser  used  to  activate  the  phonon  modes  within  the  samples  has  a  spot  size 
of  -  1  pm.  Therefore  the  Raman  response  is  an  average  response  of  not  one  specific  sample,  but 
a  representative  group  of  samples  exposed  to  the  laser. 

Figure  5  compares  the  Raman  intensity  of  3  non-heat  treated  sample  batches  (US4309, 
US4315,  and  PD15L5-20)  to  the  Raman  intensity  obtained  from  3  heat-treated  sample  batches 
(US4400,  US4406,  US4412).  The  Raman  intensity  is  normalized  with  the  maximum  G  band 
peak.  Comparison  of  Fig.  5a  to  Fig  5b  suggests  that  a  significant  amount  of  defect  healing  and 
graphitization  has  occurred  in  the  heat-treated  samples  during  to  the  3000°C  annealing  process. 

.  RESULTS  AND  DISCUSSION 

4.1.  Thermal  Conductivity  vs.  Sample  Quality 

On  a  sample  specific  level,  SEM  images  (including  Figs.  6  and  7)  are  used  to  qualitatively 
categorize  both  the  non-heat  treated  and  heat-treated  samples  into  three  main  categories:  samples 
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with  1)  minor  defects;  2)  moderate  defects;  and  3)  severe  defects.  CNT  side  growth  (branching) 
at  angles  to  the  main  CNT  tube  axis  have  been  previously  shown  to  be  caused  by  wall  defects 
[7],  Besides  branching,  the  specimens  in  the  present  study  also  show  abrupt  changes  in  diameter 
as  well  as  the  presence  of  excessive  surface  debris.  None  of  the  selected  specimens  from  either 
manufacturer  are  considered  ‘flawless.’  Not  surprising,  the  mean  thermal  conductivity  in 
specimens  from  the  severe  defects  category  is  statistically  less  than  those  of  the  specimens  from 
the  moderate  and  minor  categories.  Samples  described  by  the  minor  and  moderate  category  are 
subsequently  divided  into  heat-treated  (Fig.  6)  and  non-heat  treated  (Fig.  7)  sample  groups.  The 
severe  category  comprises  of  rejected  sample  group.  In  this  manner,  a  filtering  process  was 
performed  to  remove  samples  having  gross  morphological  deformities  from  the  study. 

A  total  of  29  measurements  of  MWCNT  conductivity  are  performed.  6  non-heat  treated  and 
15  heat-treated  MWCNTs  are  measured  to  determine  the  effect  of  high  temperature  graphitization 
on  thermal  conductivity.  8  samples  comprised  the  rejected  group.  The  mean  thermal 
conductivities  and  Raman  DIG  ratios  of  the  non-heat  treated  MWCNT  samples  provided  by  US 
Nanomaterials  Research  and  the  non-heat  treated  MWCNT  samples  provided  by  NanoLab  are 
found  not  to  have  a  statistically  significant  difference.  For  this  reason,  both  sample  sets  are 
placed  into  the  same  non-heat  treated  group. 

The  non-heat  treated  sample  batches  produce  D/G  ratios  that  are  approximately  3.3  fold 
greater  than  the  heat-treated  batches  indicating  a  significant  presence  of  average  defects  and 
amorphous  carbon  content.  Figure  8  suggests  that  samples  selected  from  the  heat-treated  batches 
produce  a  mean  thermal  conductivity  nearly  5  fold  times  (4.88x)  that  of  the  non-heat  treated 
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group,  216  ±  149  W/m-K  and  44  ±  29  W/m-K,  respectivley.  Note  the  large  standard  deviations 
are  due  to  the  variations  in  thermal  conductivity  within  the  group  and  do  not  reflect  measurement 
uncertainty.  Comparison  of  the  non-heat  treated  and  heat-treated  samples  using  a  one-way 
analysis  of  variance  (ANOYA)  resulted  in  a  p-value  of  0.015  indicating  that  the  graphitization 
process  has  a  statistically  significant  effect  on  the  thermal  conductivity  of  the  samples. 

Due  to  the  unique  morphology  of  each  sample  tested,  MWCNT  thermal  conductivity  may  be 
better  evaluated  on  a  more  sample  specific  level.  However,  the  mean  thermal  conductivity  of  all 
non-heat-treated  samples  is  found  to  be  similar  to  other  CVD  grown  MWCNTs  measured  via  the 
‘Self-Heating’  of  1-2  mm  long  MWCNT  bundles  with  estimated  phonon  mean  free  paths  of  a  few 
nanometers  [20]  and  vapor  grown  carbon  fibers  having  diameters  greater  than  10  microns  [21], 

In  the  present  study,  the  highest  measured  thermal  conductivity  is  730  ±  153  W/m-K  for 
a  heat-treated  sample  having  a  measured  OD  of  35  ±  4  nm,  estimated  ID/OD  ratio  of  0.29,  and  a 
length  of  10.77  ±  0.04  pm.  Not  surprising,  this  sample  is  the  straightest  of  all  samples  measured 
and  has  the  least  number  of  visible  physical  defects  (Fig.  6a).  A  small  amount  of  surface  debris 
is  present  on  the  sample  and  seems  fused  to  the  samples  outer  wall.  When  considering  its 
straightness,  the  sample  does  not  represent  the  majority  of  tested  specimens. 

The  greater  number  of  ‘kinks’  in  the  rejected  sample  group  likely  reduce  phonon  mean 
free  path  accounting  for  the  reduced  thermal  conductivity.  Interestingly,  all  8  rejected  samples 
came  from  the  heat-treated  sample  group.  The  rejected  sample  group  is  found  to  have  a  thermal 
conductivity  of  50  ±  28  W/m-K,  where  the  standard  deviation  is  due  to  the  large  variation  of 
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thermal  conductivities  within  the  group  and  do  not  reflect  the  measurement  uncertainty.  Since 
the  rejected  group  consists  of  all  heat-treated  samples  and  has  a  mean  conductivity  similar  to  that 
of  the  non-heat  treated  group,  the  reduction  in  thermal  conductivity  of  these  rejected  samples 
suggest  that  the  gross  physical  deformities  negate  the  effect  of  graphitization.  As  described  by 
Ref.  [7],  heat  treatment  even  at  high  temperatures  of  3000°C  is  unable  to  heal  major  defects. 
Therefore  it  is  likely  these  specific  samples  contained  severe  defects  prior  to  heat-treatment,  but 
were  not  representative  of  the  majority  of  heat-treated  MWCNTs  as  shown  by  their  decreased 
D/G  Raman  peak  ratio. 

The  measured  thermal  conductivities  do  not  correlate  to  the  sample  outer  diameter  (Fig. 
9a).  Outer  diameter  measurements  may  correlate  to  thermal  conductivity,  since  smaller  diameter 
tubes  generally  have  fewer  wall  numbers  permitting  less  inter-wall  phonon  scattering  [2], 
Following  this  reasoning,  a  single  wall  carbon  nanotube  has  the  highest  theoretical  thermal 
conductivity.  If  a  correlation  of  thermal  conductivity  to  diameter  of  the  samples  tested  in  the 
present  study  exists,  the  correlation  is  likely  to  be  unnoticeable  due  to  the  variations  in  the 
structural  quality. 

It  has  been  proposed  that  shorter  sample  lengths  may  be  responsible  for  a  reduction  in 
phonon  mean  free  paths  and  thus  apparent  thermal  conductivities  [22],  Similar  to  diameter,  the 
measurement  of  these  effects  places  a  strict  demand  on  sample  quality. 

The  magnitude  of  the  thermal  conductivities  reported  in  the  present  study  are  below  the 
highest  measured  MWCNT  of  approximately  3000  W/m-K  [1],  However,  the  variation  of 
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specimen  quality  makes  comparison  to  other  literature  values  especially  difficult.  Very  little 
information  is  available  on  the  quality  and  processing  techniques  of  the  tested  MWCNTs 
reported  in  the  literature.  The  heat-treated  sample  group  measured  a  mean  thermal  conductivity 
in  agreement  with  the  individual  MWCNT  measurement  measured  via  self-heating  [4],  The 
variation  in  an  individual  MWCNT’ s  crystallinity  and  amount  of  residual  catalyst  material 
between  samples  is  likely  to  account  for  the  wide  variations  in  the  reported  literature  values  of  an 
isolated  CNTs  thermal  conductivity  (Fig  9b). 

Comparing  the  average  20  %  relative  uncertainty  in  the  thermal  conductivity  of  each  sample 
tested  to  the  71%  variation  about  the  mean  of  the  heat-treated  sample  group  and  the  180%  variation 
about  the  mean  of  the  non-heat  treated  sample  group  suggests  that  variations  in  sample  quality 
within  each  sample  group  are  significant.  The  large  variations  are  likely  due  to  wall  defects  in  the 
heat-treated  group  and  both  wall  defects  and  residual  Fe  catalyst  in  the  non-heat  treated  group. 

4.2.  Contact  Resistance  and  Other  Sources  of  Uncertainty 

Thermal  contact  resistance  at  the  sample-probe  and  sample-manipulator  interfaces  is  said 
to  stem  from  two  possible  sources,  heat  flow  restriction  due  to  surface  roughness  of  the 
contacting  surfaces  and  mismatch  in  the  phonon  density  of  states  at  the  interface  between  the 
MWCNT  sample  and  Pt-EBID.  Since  the  probe  wire  is  composed  of  Pt,  the  contact  thermal 
resistance  is  expected  to  be  minimized  by  the  use  of  Pt-EBID.  This  is  assuming  that  vaporized 
residual  carbon  is  not  excessively  deposited  on  the  interface  with  the  SEM  electron  beam  prior  to 
the  Pt-EBID  deposition.  Deposition  of  residual  carbons  is  avoided  by  limiting  unnecessary 
imaging  at  the  EBID  site  prior  to  deposition.  Also,  increasing  the  length  of  the  Pt-EBID  pads  on 
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the  specimen  should  improve  thermal  contact  resistance  by  increasing  the  contact  area.  If 
thermal  contact  resistance  were  explicitly  known,  the  magnitudes  of  the  reported  thermal 
conductivities  would  increase.  Currently  the  probe-sample  thermal  contact  resistance  is 
attempting  to  be  measured  by  changing  the  length  of  the  sample  while  maintaining  the  original 
probe-sample  contact.  Regardless  of  the  thermal  contact  resistance  values,  comparisons  between 
sample  groups  using  the  Wollaston  T-type  probe  method  are  highly  effective. 

Excluding  the  effect  of  thermal  contact  resistance,  the  largest  source  of  experimental 
error  lies  in  the  measurement  of  sample  diameter  and  the  definition  of  an  effective  sample 
diameter  if  the  diameter  lacks  uniformity  along  the  sample  length.  The  SEM  is  capable  of 
making  repeatable  diameter  measurements  with  a  conservative  uncertainty  of  ±  4  nm.  Any 
variation  in  sample  diameter  will  add  to  this  uncertainty.  Since  thermal  conductivity  is  a 
function  of  the  sample  cross-section  area,  the  relative  uncertainty  in  sample  diameter,  8D0/D0 ,  is 
doubled.  The  average  relative  uncertainty  in  sample  diameter  for  all  measured  samples  in  this 
study  is  9  %,  and  is  dominated  mostly  by  non-uniformity  of  sample  diameter.  Other 
uncertainties  include  sample  lengths,  SL/L  (~l  %)  [23],  room  temperature  probe  electrical 
resistance,  5Reo/Reo  (<  0. 1  %),  and  temperature  coefficient  of  resistance  of  the  probe,  da/a  (2.2 
%).  These  uncertainties  in  conjunction  with  the  error  in  neglecting  radiation  heat  loss  (~  0.5  %) 
and  the  aforementioned  sensitivity  of  the  Wollaston  T-Type  probe  results  in  a  relative  uncertainty  of 
all  measurements  to  be  approximately  20  %. 
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5.  CONCLUSION 


This  study  highlights  the  beneficial  effects  of  high  temperature  heat  treatment  (~3000°C) 
on  the  thermal  conductivity  of  CVD  grown  MWCNTs  when  compared  to  non-heat  treated  CVD 
grown  MWCNTs.  A  4.9  fold  increase  in  thermal  conductivity  is  obtained  for  MWCNT  samples 
which  undergo  a  20-hour  3000°C  post  annealing  heat  treatment.  However,  in  most  samples 
defects  are  found  to  negate  any  advantage  of  the  heat  treatment  process.  Of  the  29  samples 
tested,  the  highest  measured  thermal  conductivity  is  730  ±  153  W/m-K.  Variations  in  the 
physical  structures  of  the  sample  are  seen  within  each  sample  group,  and  between  each  sample 
group  with  the  highest  values  of  thermal  conductivity  being  measured  in  the  straightest  heat- 
treated  samples. 
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APPENDIX  I 

The  heating  current  used  to  Joule  heat  the  Pt  probe  is  driven  at  a  sufficiently  low 
frequency  to  prevent  a  phase  shift  between  the  heating  frequency  and  the  temperature  response. 
This  is  achieved  by  choosing  a  heating  frequency  whose  period  is  less  than  or  equivalent  to  the 
thermal  diffusion  time  x  of  a  suspended  wire,  where  x  =  L2 /a.  To  verify  the  operational 
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frequency  range  of  the  experiment,  the  Pt  probe  wire  is  modeled  as  a  ID  wire  having  constant 
cross-section  Area,  A,  Length  L,  density  p,  and  heat  capacity  C.  Homogeneous  temperature 
boundary  conditions  0(0,t)  =  6(L,t)  =  0 ,  and  initial  condition  0{x,  0)  =  0 ,  along  with  Joule 
heating  (current  Iimeia,t )  is  solved  together  with  the  time  dependent  governing  equation 


d0{x,t)_d26{x)  Ilco2Reoei2a>t 
dt  dx2  pCLA 


(17) 


Using  the  method  of  eigenfunction  expansion,  the  spatially  average  series  solution  of  Eq.  (17)  is 
written  as 


-  1  L 

6{cd)  -  —  j  0(x,  CO)dx 


_  ^  4lJ%o  (1  -  (-!)»  )  e^_  - 
n= i  pCLA  [njp)  (i2co+(n7t)2  /  rj 


(18) 


The  frequency  dependent  thermal  transfer  function  is  then 

„  ^  _y  4  ((-y-l)(g^-^)2?-l) 

Q(t)  TApCLA  m )2  i2a>  +  (im)2  /  t 

Eq.  (19)  highlights  the  importance  of  minimizing  the  frequency  to  maintain  the  real  part  of  the 

rr\ 

thermal  transfer  function  and  thus  the  measured  3  harmonic  voltage  signal.  The  denominator  of 
the  last  term  of  Eq.  (19)  validates  that  longer  Wollaston  probe  wires  (larger  r  =  L  /a)  require 
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lower  frequencies  to  minimize  the  phase  difference  between  the  heating  current  and  the  thermal 
response.  As  frequency  approaches  zero,  the  steady  state  result  utilized  in  Eq.  (2)  is  recovered. 


Z(co^0)  =  4RthPYJ 

n=  1 


a+c-ir1) 

( ri7r )4 


^tk,P  _  ry 

12  ~  ‘ 


(20) 
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Figure  1.  Schematic  of  Pt  probe  wire  (red  line),  and  attached  sample  (/ horizontal  black  line).  The 
thermal  resistance  of  the  sample  is  incorporated  into  the  analytical  model  using  a  flux  boundary 
condition  at  x  =  0.  The  parabolic  dashed  line,  6i(x),  represents  the  increase  in  temperature  of  the 
probe  prior  to  coming  into  contact  with  the  sample.  The  solid  line,  fb(x),  represents  the  increase  in 
temperature  of  the  probe  wire  following  the  contact  with  the  sample.  The  manipulator  tip  and  each 
end  of  the  probe  wire  are  assumed  to  remain  at  ambient  temperature  conditions,  6  =  0. 
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Figure  2.  Images  of  sensor  package  (A)  during  the  Wollaston  etching  process.  The  package 
contains  two  removable  chipsets.  Each  chipset  has  two  Wollaston  probe  wires  marked  1-4.  Each 
Wollaston  wire  is  labeled  (a).  The  Wollaston  wires  are  shown  connected  to  copper  contact  pads 
using  an  electrically  conductive  silver  epoxy.  Label  (b)  marks  the  droplet  of  nitric  acid  on  the 
Teflon  spoon  used  to  etch  the  Ag  clad  from  the  Wollaston  to  reveal  the  Pt  core.  (B)  The  sensor 
package  mounted  inside  the  SEM  chamber  with  micromanipulators  marked  (c)  and  SEM  stage 
marked  (d). 
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Figure  3.  (A)  SEM  image  of  the  validation  experiment  where  a  similar  length  of  etched  Wollaston 
Pt  wire  is  used  as  both  the  T-type  probe  and  the  unknown  sample.  Image  shows  the  use  of  Pt-EBID 
deposition  to  decrease  thermal  contact  resistance.  (B)  ResmtRMs  vs.  Qrms  plot  prior  to  sample 
connection,  following  sample  connection,  and  with  the  sample  connection  improved  using  Pt-EBID. 
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Figure  4.  Scanning  electron  microscope  (SEM)  images  of  MWCNTs  batches  from  which 
samples  are  selected,  (a)  and  (b)  US  Nanomaterials  Research:  heat-treated  batch  representative 
of  all  batches  US4400,  US4406,  and  US4412;  (c)  and  (d)  US  Nanomaterials  Research:  non-heat 
treated  batch  representative  of  batches  US4309  and  US4315;  (e)  and  (f)  Nanolab:  non-heat 
treated  batch  (PD15L5-20). 
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Figure  5.  Raman  intensity  versus  wavenumber  (514  nm  excitation)  of  (a)  sample  batches  PD15L5- 
20,  US4309,  US4315  of  which  non-heat  treated  specimens  are  selected;  (b)  sample  batches  US4400, 
US4406,  US4412  of  which  heat-treated  specimens  are  selected.  The  Raman  intensity  is  normalized 
by  the  maximum  value  of  the  G  peak  for  clarity. 
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Figure  6.  SEM  micrographs  of  experiments  from  the  heat-treated  sample  set.  In  each  case  the  Pt 
probe  wire  is  pictured  above  the  manipulator  tip  with  sample  joining  the  probe  wire  to  the 
manipulator.  (A)  Highest  measured  thermal  conductivity  (730  ±152  W/m-K),  statistically  an 
outlier;  (B)  and  (C)  the  two  highest  measured  values  of  the  heat-treated  samples  that  are  within 
the  20%  confidence  interval  about  the  groups  mean,  290  ±  28  W/m-K  and  254  ±  42  W/m-K, 
respectively;  (D)  and  (E)  the  two  lowest  measured  values  of  the  heat-treated  samples,  80  ±  16 
W/m-K  and  67  ±11  W/m-K,  respectively. 
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Figure  7.  SEM  micrographs  of  experiments  resulting  in  the  two  highest  (A  and  B)  and  two 
lowest  (C  and  D)  measured  thermal  conductivity  values  of  the  non-heat  treated  sample  group. 
(A)  102  ±  17  W/m-K;  (B)  61  ±  6  W/m-K;  (C)  24  ±  4  W/m-K;  (D)  21  ±  9  W/m-K.  In  each  case 
the  Pt  probe  wire  is  pictured  above  the  manipulator  tip  with  sample  joining  the  probe  wire  to  the 
manipulator. 
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Figure  8.  Mean  thermal  conductivity  of  the  three  measured  sample  groups  plotted  versus  the 
ratio  of  the  D  to  G  Raman  spectroscopy  peaks  as  measured  on  each  respective  sample  batch 
using  a  514  nm  excitation  source.  Error  bars  correspond  to  the  standard  deviation  calculated 
from  the  variance  of  all  samples  within  each  sample  group  and  do  not  reflect  the  measurement 
uncertainty. 
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Figure  9.  (A)  Measured  thermal  conductivities  of  the  heat-treated,  non-heat  treated,  and  rejected 
sample  groups;  (B)  Thermal  conductivity  measurements  of  both  heat-treated  and  non-heat 
treated  sample  groups  plotted  against  references  [1-4]. 
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APPENDIX  II 


THIC  NESS  AND  TEMPERATURE  DEPENDENT  THERMAL  CONDUCTIVITY  OF 

NANOSCALE  TIN  FILMS 


ABSTRACT 

Thin  films  in  general  exhibit  different  thermal  properties  compared  to  bulk  due  to  size  effect  [1-3].  In 
this  study,  the  thermal  conductivity  of  sputtered  Sn  films  of  thickness  500  nm  ±  50  nm  and  100  nm  ±  20 
nm  are  obtained  from  55K  to  300K  and  from  40K  to  31  OK,  respectively,  using  the  three  omega  method. 
The  thermal  conductivity  of  500  nm  thin  film  at  room  temperature  is  46.2  ±4.2  W/m-K,  which  is  lower 
when  compared  to  its  bulk  value  of  63  W/m-K,  and  increases  gradually  as  the  temperature  is  lowered  to 
55K.  In  contrast,  the  thermal  conductivity  of  the  100  nm  thin  film  exhibits  even  reduced  thermal 
conductivity,  36  ±  2.88  W/m-K  at  300K,  when  compared  to  the  500  nm  film,  and  decreases  as  the 
temperature  is  lowered.  The  reduction  in  thermal  conductivity  of  Sn  thin  film  may  be  due  to  the 
pronounced  effects  of  electron  scattering  at  the  grain  boundaries  as  well  as  the  twin  boundaries  in  addition 
to  the  scattering  from  the  boundary  surface  at  lower  temperatures.  These  experimentally  determined 
thermal  conductivities  are  compared  to  models  that  take  into  account  size  effects  on  thermal  conductivity 
of  metallic  films  based  on  electronic  scattering  as  proposed  by  Fuchs-Sondheimer  (FS),  Mayadas- 
Shatzkes  (MS)  and  Qiu  and  Tien  (QT).  The  experimentally  measured  thermal  conductivity  of  Sn  films  is 
in  good  agreement  with  the  MS  model  indicating  the  importance  of  the  grain  boundary  scattering. 
Thickness  measurements  are  obtained  by  ellipsometry  and  profilometer.  The  estimation  of  the  mean  grain 
size  in  both  films  and  the  evidence  of  twin  boundaries  are  obtained  by  Atomic  Force  Microscopy. 


INTRODUCTION 

Thin  films  (metal  or  dielectric)  are  essential  constituents  in  current  functional  devices  having  myriad 
of  applications  relevant  to  energy  transport  and  energy  conversion  technologies  such  as  thermoelectrics, 
Nano-Electro-Mechanical  Systems  (NEMS),  optoelectronics,  solar,  photovoltaics  etc.  Understanding  the 
thermal  transport  and  more  fundamentally  the  behavior  of  heat  carriers  (electron,  photon  and  phonon)  in 
these  materials  is  critical  to  their  reliable  and  robust  performance,  thus  accommodating  denser  circuits 
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and  higher  switching  speeds.  Nanoscale  thin  films  of  tin  (Sn),  are  important  as  a  solder  material  in 
microelectronics  applications  and  as  a  thermal  interface  material  to  improve  thermal  conductance  between 
two  non-conforming  surfaces.  This  is  due  to  its  properties  of  reasonably  good  thermal  conductivity  (64 
W/m-K  at  room  temperature)  [4],  lower  melting  point  (505  K)  [5]  and  its  ability  to  conform  well  to 
surface  irregularities.  The  total  thermal  resistance  of  a  thermal  interface  material  (TIM)  depends  on  its 
thermal  conductivity,  bond  line  thickness  (BLT)  and  the  contact  resistances  of  the  TIM  with  the  two 
bounding  surfaces.  In  this  regards,  a  high  thermal  conductivity  material  and  lower  bond  line  thickness 
sufficient  to  fill  the  interfacial  gap  between  the  asperities  is  desirable  to  improve  the  thermal  contact 
conductance. 

Reduced  thermal  conductivity  at  nanoscale  due  to  size  effects  was  initially  reported  [6]  in  the 
resistivity  measurements  of  metallic  thin  films  specifically  when  the  film  thickness  was  comparable  to  the 
mean  free  path  of  electrons.  This  anamoly  was  first  explained  by  Fuchs-Sondheimer  theory  [7], 
attributing  it  to  the  scattering  of  electrons  at  the  surface.  However,  the  Mayadas-Shatzkes  theory  [8], 
proposed  much  later,  was  in  better  agreement  with  the  experimental  results  that  showed  increase  in 
resistivity  with  decreasing  film  thickness.  Even  though  both  theories  are  fundamentally  derived  from  the 
well-known  Boltzmann  transport  equations,  the  MS  theory  also  accounts  for  the  grain  boundary  scattering 
in  addition  to  the  surface  boundary  scattering  of  electrons.  Recent  experimental  [9-10]  and  modeling  [11- 
12]  studies  on  nanoscale  polycrystalline  metallic  thin  films  reported  reduced  thermal  conductivities 
compared  to  their  bulk  value  and  additionally  in  some  cases,  the  violation  of  Wiedemann-Franz  law.  The 
most  widely  studied  metallic  films  at  nanoscale  are  aluminum,  copper  and  gold.  With  regards  to  Sn,  the 
electrical  resistivity  of  tin  films  from  2  mm  down  to  3  pm  thickness  have  been  investigated  thoroughly  up 
to  very  low  temperatures  [6].  The  electrical  resistivity  has  been  observed  to  be  increasing  with  decreasing 
tin-foil  thickness  and  with  increasing  temperature.  The  shortening  of  the  mean  free  path  due  to  inelastic 
collision  of  the  conduction  electrons  with  the  boundary  (lateral)  surfaces  of  the  metal  have  been 
understood  to  increase  the  resistivity  at  thicknesses  comparable  to  the  electron  mean  free  path.  However, 
there  has  been  limited  investigation  into  the  thermal  properties  of  tin  particularly  as  a  function  of 
temperature.  Previous  studies  have  reported  the  measurements  of  thermal  conductivity  of  pure  tin 
(99.999%,  bulk)  either  at  very  low  temperatures,  3.5K  -  25K  [13],  or  at  very  high  temperatures  373K- 
873K  [14].  Recently,  Choi  et  al  [15],  have  reported  the  thermal  diffusivity  measurement  of  2pm  Sn  film. 

Transient  methods  for  the  measurement  of  thermal  transport  properties  are  widely  accepted  and  can 
be  divided  into  frequency-domain  and  time-domain  techniques.  The  most  popular  methods  under  these 
domains  are:  (1)  three  omega  method  that  employs  an  electrically  heating  source,  a  metal  line 
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heater/sensor,  and  a  moderate  frequency  modulation  range  (0.1Hz-200KHz);  and  (2)  thermoreflectance 
methods  (TDTR  and  FDTR)  that  use  ultra-fast  laser  pulses  with  modulating  frequency  up  to  80  MHz.  The 
high  frequency  range  is  advantageous  for  the  measurement  of  heat  diffusion  in  very  thin  (angstrom  scale) 
films  and  across  interfaces.  The  three  omega  method,  developed  by  Cahill  et  al  [16].,  is  the  most  reliable 
and  versatile  technique  for  the  measurement  of  thin  films  and  bulk  materials.  The  method  circumvents 
errors  due  to  radiation  and  convection  and  its  implementation  is  relatively  inexpensive.  Recent 
modifications  to  Cahill’s  original  three  omega  method  [16]  has  extended  its  applicability  to  electrically 
conducting  and  semiconducting  films  [17-20],  multi-layer  structures  including  composites  [2,  21-23,  24- 
25]  in  addition  to  estimating  the  anisotropic  thermal  conductivities  of  a  thin  film  [  21,23]  and  the 
interfacial  thermal  resistance  [  2,  21-23]. 

In  this  work,  the  differential  three  omega  method  [21-23]  is  employed  for  the  measurement  of  thermal 
conductivity  of  soft  nanoscale  tin  (Sn)  films  of  thickness,  500nm,  from  55K  to  300K  and  lOOnm,  from 
40K  to  31  OK  respectively.  The  thermal  conductivity  of  these  films  is  compared  to  a  previous  study  to 
highlight  the  influence  of  thickness  on  thermal  transport.  The  experimental  data  is  then  compared  to 
previously  proposed  models  based  on  the  electron  mean  free  path  comparable  to  thickness  or  average 
grain  size.  The  temperature  dependent  thermal  conductivity  of  these  films  is  discussed  by  examining  the 
electron  and  phonon  contribution  separately.  Additionally,  a  brief  overview  of  three  omega  method,  the 
sample  preparation  including  the  shadow  mask  microfabrication  and  the  calibration  of  deposited  metal 
line  is  presented. 

THREE  OMEGA  TECHNI  UE:  A  BRIEF  OVERVIE 

In  this  section,  a  brief  overview  of  the  general  three  omega  technique  is  presented  to  better 
illustrate  the  required  changes  for  accommodating  soft  metallic  films.  The  three  omega  method  generally 
proceeds  by  applying  an  alternating  current  of  frequency  through  a  metal  heater  line  (see  Fig.  1)  that 
has  been  directly  deposited  on  an  electrically  insulating  sample.  This  current  heats  the  sample  at  a 
frequency  of  2  co  due  to  Joule  heating,  producing  temperature  oscillations  also  at  a  frequency  of  2&>with 
an  amplitude,  say,  A T1(a ,  and  phase  difference  c/>  .  Since  the  resistance  of  pure  metals  increases  linearly 

with  temperature,  the  temperature  oscillations  introduce  resistance  oscillations  in  the  metal  line  at  a 
frequency  of  2  co .  The  resistance  oscillations  at  2  co ,  with  the  source  current  at  frequency  co ,  generate  a 
small  oscillating  voltage  signal  across  the  metal  line  at  3  co .  The  amplitude  of  the  3  co  voltage  V3co  is  given 

in  terms  of  the  voltage  amplitude  at  the  first  harmonic,  V(0  ,  by 
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K  =  —I.R  aAO  =  —  V  aAT2 

Here,  iXm  is  the  current  sent  to  the  metal  line,  R  is  the  average  resistance  and  a T2co  represents  the 

temperature  oscillations  of  the  metal  line  at  frequency  2  co  .  The  temperature  coefficient  of  resistance  a 
needs  to  be  calibrated  for  every  sensor  microfabricated  on  the  sample  at  the  start  of  experiment. 


Fig.l  Block-diagram  representing  the  three  omega  electrical  circuit. 


In  the  present  work,  in  order  to  determine  the  temperature  rise  in  the  deposited  metal  line  the  third 
harmonic  voltages  were  isolated  using  a  lock-in  amplifier  (Stanford  Research  Systems  SR850)  which 
operates  over  a  wide  range  of  frequencies  (1  mHz- 102  kHz)  in  conjunction  with  a  custom  designed 
circuitry  based  on  Fig.l.  Unity  gain  differential  amplifiers  were  employed  to  acquire  the  voltage  signal 
across  the  potentiometer  and  the  metal  line  (that  acts  as  a  heater/  sensor)  on  the  sample  and  a  12  bit 
multiplying  digital-to-analog  converter,  DAC  (Analog  Devices  No.  7541KN)  in  conjunction  with  an 
operational  amplifier  (Analog  Devices  No.  549KH)  so  that  the  first  harmonic  differential  signal  could  be 
balanced  out.  The  set-up  has  been  automated  using  Lab  VIEW  8.5.  The  Lab  VIEW  program  also  helps  in 
balancing  the  first  harmonic  voltage  at  the  reference  resistor  and  the  sample  resistor  to  a  greater  degree  of 
precision.  The  bits  of  the  multiplying  DAC  are  set  4 on’  and  4 off’  to  balance  the  first  harmonic  signal  from 
the  sample.  This  step  is  equivalent  to  making  adjustment  on  the  potentiometer  manually.  However,  when 
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automated,  the  digital  signal  precision  is  close  to  4  decimal  places,  reading  VA-b  on  the  lock-in  equal  to 
±0.000  IV.  A  multi-DAQ  card  (NI-6025E)  is  used  for  communicating  with  the  multiplying  DAC  through 
Lab  VIEW.  The  third  harmonic  voltages  at  a  specified  temperature  were  measured  by  performing  a 
logarithmic  frequency  sweep  typically  from  50  Hz  to  10  kHz.  Finally,  the  thermal  conductivity  of  the 
sample  was  evaluated  by  measuring  the  frequency  dependence  of  the  temperature  oscillations  in  the  metal 
line.  By  correlating  the  experimentally  determined  voltage  to  the  analytically  derived  temperature 
changes,  the  thermal  conductivity  of  the  sample  was  determined  without  the  knowledge  of  the  film’s 
density  and  specific  heat. 

The  analytical  derivation  based  on  an  earlier  work  by  Cahill  and  Pohl  [26]  requires  the 
assumptions  that  the  heater  be  considered  of  infinite  length  and  that  the  sample  be  effectively  semi¬ 
infinite.  Additionally,  the  method  requires  the  thermal  waves  to  penetrate  a  depth  greater  than  the  heater 
width  into  the  sample,  satisfying  the  approximation  of  an  infinite  line  source.  The  three  omega  method 
can  also  be  applied  to  a  film  on  a  substrate  type  configuration,  by  employing  either  the  slope  method  or 
the  differential-3  co  approach.  For  these  methods,  the  thin  film  is  assumed  thin  enough  such  that 
essentially  one  dimensional  Fourier  conduction  model  can  be  applied  for  analysis  and  interpretation  of 
experimental  data.  Specific  conditions  that  must  hold  to  ensure  accurate  analysis  are  that  df  >  2 w  » tf , 

i.e.,  the  thermal  penetration  depth,  d f ,  must  be  greater  than  the  heater  width  and  be  much  larger  than  the 
film  thickness,  tf .  For  the  experimental  investigation  discussed  herein,  the  thermal  penetration  depth  and 
the  heater  width  were  greater  by  at  least  an  order  of  magnitude  than  the  film  thickness  to  ensure  the 
validity  of  these  approximations.  Additionally,  the  thickness  of  the  underlying  substrate  ts  must  be  such 
that  the  condition  ts  »  df  »  w  is  maintained  to  avoid  back  surface  reflection  and  assure  the  validity  of 

the  line  source  and  semi-infinite  substrate  approximations.  Additional  details  about  the  three  omega 
technique  can  be  found  from  various  references  [1-3,  16-18,  21-23,  26]  and  are  therefore  not  discussed  in 
detail. 

EXPERIMENT 

For  the  experimental  procedure,  the  differential  three  omega  technique  [21-23]  was  utilized  to 
measure  the  thermal  conductivity  of  a  (-500  nm)  thick  Sn  film  and  a  (-100  nm)  relatively  thin  Sn  film, 
each  deposited  on  a  silicon  substrate.  The  assumptions  described  above  were  validated  by  selecting  an 
appropriate  frequency  range  to  ensure  sufficient  penetration  depth  and  a  suitable  heater  dimension  to 
satisfy  the  inequality  conditions.  Additionally,  in  order  to  restrict  the  error  introduced  by  various 
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approximations,  the  numerically  obtained  exact  solution  of  the  film  on  substrate  system,  as  described  by 
Kim  et  al  [2],  Borca-Tascuic  et  al  [21]  and  more  recently  by  Tong  et  al  [23],  has  been  utilized. 

A.  SAMPLE  PREPARATION 

Thin  films  of  Sn  (99.999%  pure,  Target  from  Kurt  Lesker)  of  required  thickness  was  deposited  on  a 
silicon  substrate  by  sputtering  (Denton  Vacuum,  Discovery  18)  at  a  base  pressure  of  2  x  10  7  Torr.  Since 
tin  has  lower  melting  point,  the  target  was  procured  with  a  copper  back  plate  for  dissipating  heat  during 
the  deposition  process.  The  thickness  of  the  films  was  measured  using  a  Dektak  stylus  profilometer.  Tin 
being  electrically  conducting,  an  insulating  layer  was  needed  so  that  heater/  sensor  metal  line  does  not 
come  in  direct  contact  with  the  underlying  tin  layer  for  measurement.  In  order  to  prevent  heating  of  the  tin 
thin  film,  low  temperature  deposition  of  -300  nm  Si02  was  achieved  using  a  Plasma  Enhanced  Chemical 
Vapor  Deposition  (PECVD)  process.  A  typical  PECVD  process  operates  at  300  °C  or  higher,  but  tin  melts 
at  a  temperature  of  230  °C,  lower  than  this.  Therefore,  a  procedure  described  elsewhere  (for  polymer 
samples)  [20]  to  deposit  the  Si02  using  a  tetraethylorthosilicate  (TEOS)  precursor  was  employed  at  a 
processing  temperature  of  150°C.  Before  every  deposition  an  autoclean  procedure  and  outgassing  steps 
were  executed  to  remove  any  source  of  contamination  in  the  machine.  The  final  condition  of  the  sample 
was  confirmed  after  the  metal  line  deposition  revealed  no  continuity  between  heaters  located  on  the  same 
sample  surface.  A  resistance  value  of  the  order  106  indicated  that  the  oxide  layer  provided  sufficient 
insulating  capacity  to  run  the  experiment. 

The  aluminum  metal  heater  line  deposition  was  achieved  by  magnetron  sputtering  using  a 
microfabricated  shadow  mask  as  shown  in  Fig. 2.  The  shadow  mask  is  used  in  order  to  avoid  the 
conventional  photolithography  procedure  every  time  in  addition  to  making  it  available  for  a  wide  variety 
of  materials  that  cannot  survive  the  wet  etching  process.  Aluminum  heaters  of  widths  25pm,  50pm  and 
100pm  were  deposited  using  these  shadow  masks.  After  heater  line  deposition,  the  diced  samples  were 
placed  in  a  24-pin  side  braze  package  so  that  wires  could  be  bonded  to  the  heater  pads  using  either 
ultrasonic  bonding  or  epoxy  bonding.  The  heater  width  and  length  were  verified  by  the  non-contact 
Nanovea  optical  profilometer  (Fig. 3). 
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Fig.  2: (a) A  Schematic  representation  of  the  footprint  of  stencil  mask  on  a  4”  Si  wafer.  The  design 
accommodates  16  individual  re-usable  Si  shadow  masks  with  heaters  of  widths  25pm  to  100pm  etched  on  it 
using  Deep  Reactive  Ion  Etching  (DRIE)  process,  (b)  Optical  microscopy  image  of  an  etched  25pm  heater  on 
the  shadow  mask. 


Fig.3:  The  intensity  scan  image  (from  Nanovea  profilometer)  showing  a  25  pm  wide  heater  wire- 
bonded  on  the  Sn-thin  film  (500  nm)  sample.. 
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B.  SENSOR  CALIBRATION 


At  the  start  of  each  experiment,  the  deposited  heater/sensor  wire  lines  on  the  S  thin  film  samples  were 
calibrated  (as  shown  in  Fig.  4).  In  the  experiments,  liquid  helium  low  temperatures  are  preferred  because 
of  the  longer  mean  free  path  of  the  heat  carriers,  i.e.  the  electrons  and  phonons,  and  hence  a  lucid 
explanation  to  phenomena  impeding  the  thermal  transport.  At  higher  temperatures,  these  effects  may  be 
masked  by  Umklapp  scattering.  Low  temperatures  for  this  set-up  are  achieved  by  means  of  a  cryostat 
(Janis  Research  Model:  CCS-400H/204).  The  cryostat  is  capable  of  operating  between  the  temperature 
range  10K-500K,  and  with  special  thermometry  can  go  up  to  800K.  It  is  critical  to  reduce  the  heat  load  to 
achieve  greater  cooling  in  less  time  especially  when  working  at  temperatures  below  100  K.  The  heat  load 
is  reduced  by  keeping  the  vacuum  pressure  on  the  order  10'3-  10'5  torr  while  the  He  circulates  through  the 
cryostat.  This  also  helps  in  analysis  by  minimizing  the  convection  effects  that  would  possibly  introduce 
errors  during  measurement.  An  RV-8  rotary  vane  pump  is  used,  which  is  capable  of  developing  vacuum 
equal  to  10'4  torr  or  less. 

The  variability  in  temperature  at  the  sample  mount  is  accomplished  by  a  control  heater  at  the  cold 
finger  by  means  of  a  calibrated  silicon  diode  DT-670B  sensor  and  a  Lakeshore  332S  temperature 
controller.  The  silicon  diode  sensor  is  used  for  both  controlling  and  monitoring  the  temperature  at  the 
sample  mount.  The  temperature  controller  is  operated  in  a  closed  control  loop  mode  with  PID 
(Proportional  Integral  Derivative)  tuning  of  the  heater  power  based  on  the  feedback  signal  from  the 
silicon  diode  sensor. 
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Fig.4.  Resistance  Vs  Temperature  calibration  curve  of  a  25  jam  wide  heater  on  500nm  Sn  thin  film. 


For  the  sensor  calibration,  the  samples  were  placed  in  the  cryostat  and  a  probe  signal  from  the 
lock-in  amplifier,  at  a  relatively  small  amplitude  voltage,  is  sent  to  avoid  significant  heating  of  the 
sample.  At  each  temperature  excursion,  the  resistance  of  the  metal  line  is  measured  and  finally  plotted 
against  the  measured  temperature.  The  slope  of  the  measurement  yields  the  temperature  coefficient  of 
resistance  (TCR)  of  the  deposited  metal  line  sensor.  The  repeatability  of  the  sensor  was  assured  by  the 
fact  that  the  change  in  resistance  for  a  unit  temperature  change  varied  within  -0.1%  between  the  cooling 
and  heating  excursions  for  a  heater  with  25pm  heater  width  and  length  2.5mm.  Following  the  procedure 
described  above,  the  TCR  of  a  25pm  width  aluminum  metal  line  on  a  (~500nm  )  Sn  thin  film  is  1/44.708 
or  0.02237  K'1,  with  slight  deviation  from  linearity  at  -  60  K. 

RESULTS  AND  DISCUSSION 

Thermal  conductivity  as  a  function  of  temperature  of  (~500nm)  thick  and  (-lOOnrn)  thin  Sn  films  was 
measured  using  the  three  omega  method.  An  insulating  layer  of  ~300nm  Si02  was  deposited  using  a  low 
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temperature  PECVD  process.  Metal  heaters  of  aluminum  were  deposited  on  both  samples  through  shadow 
mask  using  magnetron  sputtering.  The  samples  were  finally  packaged  and  wire  bonded  using  ultrasonic 
bonder  or  epoxy  bonding.  The  resistance  value  in  mega-ohms  between  any  two  deposited  heaters  ensured 
the  avoidance  of  pin  holes  and  therefore  reliable  measurement.  Metal  line  deposition  was  calibrated  as  a 
function  of  temperature  after  placing  the  samples  in  the  cryostat. 

The  thermal  conductivity  in  a  polycrystalline  metallic  film  has  contribution  from  both  electrons  and 
phonons.  At  very  low  temperatures,  electrons  are  the  dominant  heat  carriers.  A  finite  mean  free  path  of  the 
carriers  is  caused  by  the  thermal  vibrations  of  the  lattice  of  the  order  of  the  wavelength  of  the  sound 
waves  in  the  metals.  The  mean  free  path,  larger  than  the  interatomic  distances,  is  increased  by  lowering 
the  temperature.  However,  the  mean  free  path  does  not  increase  indefinitely  as  the  electrons  experience 
inelastic  scattering  by  the  lattice  imperfections,  impurities  and  the  boundary  surface.  Phonon  and  electron 
mean  free  paths  are  therefore  estimated  for  bulk  tin  in  order  to  understand  the  thermal  conductivity  of 
(~500nm)  thick  and  (~100nm)  thin  Sn  films. 

A.  PHONON  AND  ELECTRON  MEAN  FREE  PATHS  IN  TIN  FILMS 

The  average  phonon  mean  free  path  (Table  1)  based  on  the  kinetic  theory  and  using  the  available 
literature  values  for  bulk  thermal  conductivity,  density  and  specific  heat  is  ~  65  nm.  The  average  random 
velocity  of  phonons  is  taken  to  be  the  dilatation  wave  speed  in  Sn.  The  electron  mean  free  path  in  bulk  tin 
(Table  2),  is  estimated  as  a  product  of  Fermi  velocity  and  the  relaxation  time  -  the  average  time  a  heat 
carrier  travels  before  it  is  scattered  and  changes  its  direction,  since  electrons  at  the  Fermi  level  are 
generally  responsible  for  conduction.  The  relaxation  time  in  bulk  tin  is  estimated  by  an  empirical  relation 
that  uses  experimentally  determined  resistivity  values  at  77  K,  273  K  and  373  K.  Additionally,  Andrews 
[4]  reported  the  experimentally  observed  electron  mean  free  path  in  bulk  tin  as  100  nm  at  79  K  and  70  pm 
near  4  K.  However,  at  77  K,  the  electron  mean  free  path  is  relatively  lower  -21  nm  and  this  is  mostly  due 
to  the  fact  that  the  value  of  radius  ratio  rs  /  a0  and  Fermi  velocity  for  bulk  tin  is  available  in  literature  only 

at  298K. 
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Table  1.  Phonon  mean  free  path  in  bulk  tin  (at  298K) 
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7365b 

2261 

202d 

55- 

2064 

210e 

OS 

5760° 

2557 

202 

60- 

2334 

210 

71 

Elastic  Modulus,  E  =  32  GPafand  i>  =  0.25 
a,c,d[4]  b[5]  e[27]  f[28]. 

Table  2.  Electron  mean  free  path  in  bulk  tin  (at  4K,  77K,  273K,  373K) 


Temperature 

373Kb 

4Ka 

77Kb 

273Kb 

Pm 

2.1 

10.6 

15.8 

(pQ-cm) 

1.1 

0.23 

0.15 

T  (xlO'14 

sec) 

70pm 

20.68 

4.3 

2.82 

mean 

100  nmc 

free  path 

/  nm 

^0.22^ 


V  Pm  J 


r  V 


IQ-14 ;  /  =  vft  lfs/a0  =2.22  at  300K  (for  tin);  vF  =1 .9xl06  m/s;  a[6] 


[4]  c[6]  and  at  79K. 
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B.  ATOMIC  FORCE  MICROSCOPY  OF  TIN  FILMS 


The  500  nm  and  the  100  nm  thick  polycrystalline  Sn-thin  films  were  analyzed  using  atomic  force 
microscope  to  estimate  the  grain  size.  A  rough  estimate  of  the  size  scale  can  also  be  made  from  the  height 
scan  of  films  with  trace-retrace  profiles  of  the  AFM  probe,  as  shown  in  Figs  5c  and  5d  for  the  500  nm  and 
100  nm  tin  films,  respectively.  For  the  same  scan  area  10pm  x  10pm  and  resolution  (Fig.5c,  d),  the  grain 
density  is  higher  in  100  nm  Sn  films  compared  to  500  nm  Sn  films.  The  grain  size  range  for  the  500  nm 
thick  Sn  film  is  450  to  2000  nm,  while  the  average  grain  size  for  the  100  nm  thick  Sn  film  is  relatively 
lower,  and  is  in  the  150  -  600  nm  range.  Interestingly,  a  careful  examination  of  AFM  images  (Fig.5a,b 
and  Fig.5c,d)  also  reveal  the  presence  of  twin  boundaries  in  both  Sn  films  more  often  than  grain 
boundaries.  Moreover,  the  density  of  twin  boundaries  in  100  nm  thin  film  is  much  greater  than  those  in 
500  nm.  The  twin  boundaries  are  estimated  to  be  of  the  order  lower  than  75  nm.  The  presence  of  twin 
boundaries  has  been  reported  more  commonly  for  the  group  IV  elements  like  silicon  having  fcc/diamond 
structure. 
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Fig.  5:  AFM  scan  images  of  nanoscale  Sn  films:  (a)  Amplitude  scan  of  500nm  Sn  film;  (b)  Amplitude 
scan  of  lOOnm  Sn  film;  (c)  Height  scan  of  500nm  Sn  film  and  (d)  Height  scan  of  lOOnm  Sn  film.  It  is 
evident  that  lOOnm  Sn  film  has  significantly  higher  density  of  grain  boundaries  and  twin  boundaries  than 
the  500nm  Sn  film. 


C.  TEMPERATURE  DEPENDENT  THERMAL  CONDUCTIVITY  OF  TIN  FILMS 

For  the  work  presented  in  here,  the  thermal  conductivity  of  the  500  nm  ±  50  nm  thick  films  was 
measured  from  55  K  to  300  K  while  the  thermal  conductivity  of  100  nm  ±  20  nm  thin  films  was  measured 
from  40K  to  31  OK  (Fig.6). 


Fig. 6.  Thermal  conductivity  of  a  500nm  Sn  film  (55K-300K)  and  a  lOOnm  Sn  film  (40  K-310  K) 
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The  thermal  conductivity  of  500  nm  Sn  film  is  46  ±  4.2  W/m-K  at  300K  and  increases  gradually 
as  the  temperature  is  lowered  down  to  100  K.  As  the  temperature  is  lowered  further,  the  thermal 
conductivity  exhibits  a  sharp  increase  and  gradually  rolls  off  around  50  K  at  which  the  thermal 
conductivity  is  131  ±  17  W/m-K.  The  electron  mean  free  path  in  bulk  tin  at  room  temperature  is  very  low 
(less  than  20  nm)  as  explained  earlier  so  it  may  not  contribute  to  the  heat  conduction  markedly.  The 
Debye  temperature  of  Sn  is  around  170  K,  beyond  which  the  specific  heat  reaches  its  Dulong-Petit  limit  at 
~  210  K.  This  indicates  that  the  gradual  increase  in  thermal  conductivity  may  be  solely  due  to  the  gradual 
increase  in  the  mean  free  path  of  the  phonons  and/or  electrons  as  the  temperature  is  lowered  while  the 
specific  heat  does  not  change  (decrease)  much  up  to  120  K  (The  specific  heat  decreases  by  less  than  8% 
of  its  value  at  170  K).  In  order  to  understand  the  relative  contribution  of  the  phonons  and  electrons  to  the 
thermal  conductivity, of  Sn,  the  temperature  at  which  the  phonon  contribution  begins  to  exceed  the 
electronic  contribution  can  be  estimated  using  using  the  relation  [4], 


T  =  0.145 


/  \l/2 

\^F  J 


(0) 


where  0D  is  the  Debye  temperature,  TF  is  the  Fermi  temperature,  and  z  is  the  nominal  valence.  This 

temperature  for  Sn  can  be  estimated  as  188.7  K.  Lower  than  170K,  even  though  the  specific  heat 
decreases  gradually,  the  electron  mean  free  path  increases  at  a  high  enough  rate  to  contribute  to  heat 
conduction.  As  described  previously,  the  electron  mean  free  path  is  100  nm  around  80K  and  increases 
sharply  to  70  pm  at  4  K.  Therefore,  the  sharp  increase  in  thermal  conductivity  can  be  attributed  to  the 
jump  in  electron  mean  free  path,  and  gradually  rolls  off  around  50  K  at  which  point  the  electronic  heat 
carriers  experience  inelastic  scattering  from  the  lateral  boundary  surface  (500nm  ±  50nm),  and 
consequently  the  electron  mean  free  path  decreases. 

In  a  sharp  contrast  to  the  500  nm  thick  Sn  film,  the  thermal  conductivity  of  the  100  nm  Sn  thin  film  is 
determined  to  decrease  gradually  as  the  temperature  is  lowered.  The  thermal  conductivity  of  the  100  nm 
thin  film  is  36  ±  2.88  W/m-K  at  room  temperature  and  decreases  with  temperature  up  to  200  K,  near  the 
Debye  temperature  of  Sn,  where  it  plateaus  and  then  decreases  all  the  way  down  to  40  K.  The  gradual 
decrease  in  thermal  conductivity  from  room  temperature  can  be  attributed  to  the  to  the  scattering  of 
phonons  by  the  grain  boundaries  and  the  twin  boundaries  as  well  as  the  inelastic  scattering  at  the 
boundary  surface  (lOOnm  ±  20nm)  since  the  phonon  mean  free  path  -  65  nm  is  comparable  to  the  twin 
boundaries  more  than  the  grain  boundaries  in  addition  to  the  thin  film  thickness.  It  is  to  be  noted  that 
density  of  the  twin  boundaries  is  significantly  lower  and  the  grain  size  is  noticeably  higher  in  the  500  nm 
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thick  Sn  film  (as  evidenced  by  the  AFM  images  described  earlier),  and  therefore  their  presence  does  not 
markedly  affect  the  phonon  mean  free  path  as  in  the  case  of  100  nm  thin  film.  Moreover,  the  scattering 
effect  of  boundary  (lateral)  surface  compounds  the  effect  of  scattering  from  the  twin  boundaries  and  grain 
boundaries  in  the  case  of  the  Sn  100  nm  thin  film.  Consequently,  as  the  temperature  is  lowered  from  the 
Debye  temperature,  the  thermal  conductivity  decreases  both  due  to  lower  specific  heat  and  the  effect  of 
mean  free  path.  At  around  80  K,  the  electron  mean  free  path  is  also  long  enough  to  experience  the 
scattering  effects  due  to  twin  boundaries  and  the  boundary  surface.  In  this  way,  the  thermal  conductivity 
in  the  case  of  the  100  nm  Sn  film  decreases  continuously  due  to  the  combined  effects  of  scattering  of 
majority  carriers  i.e.  electrons,  and  the  lower  lattice  specific  heat  as  not  many  phonon  modes  are  excited 
at  temperatures  very  low  from  the  Debye  temperature  limit . 
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Fig.7.  Complete  thermal  conductivity  profile  of  Sn  including  thin  film  and  bulk  from  10K  to  1473K. 
Note  the  phase  change  from  solid  Sn  to  liquid  Sn  occurs  at  505K. 


Most  of  the  previous  thermal  conductivity  studies  [13,  14,  29,  30]  performed  on  tin  have  been  on  bulk 
tin  and  have  been  conducted  either  at  temperatures  near  the  liquid  helium  temperature  or  at  very  high 
temperatures  near  its  molten  state  i.e.  phase  transition  from  solid  tin  to  liquid  tin.  To  our  knowledge,  Choi 
et  al  (2006)  [15]  was  first  to  measure  thermal  properties  of  deposited  Sn  thin  film  of  thickness  ~  2000  nm. 
This  film  was  deposited  using  the  dc  magnetron  sputtering  method  at  150  W.  The  thermal  diffusivity  was 
measured  at  room  temperature  using  a.c.  calorimetry  and  mirage  method.  The  measured  thermal 
diffusivity  of  2000  nm  Sn  film  was  39  10'6  m2/s,  which  is  equivalent  to  the  thermal  conductivity  of  60.3 
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W/m-K,  which  is  closer  to  the  bulk  value  of  64  W/m-K  at  room  temperature.  In  the  present  study,  we 
report  the  first  measurements  of  thermal  conductivity  in  nanoscale  Sn  thin  films  500nm  and  lOOnm  as  a 
function  of  temperature.  The  current  work  is  compared  to  the  previous  studies  done  on  bulk  tin  and  a 
complete  thermal  conductivity  profile  of  Sn  as  a  function  of  temperature  is  presented  from  10  K  to  1473 
K  for  Sn  thin  film  and  bulk  (Fig. 7).  We  did  not  find  any  thermal  conductivity  measurements  of  bulk  tin  at 
intermediate  temperatures  from  22  K  to  293  K.  Note  that  the  thermal  conductivity  of  500  nm  thin  film  is 
lower  than  the  expected  bulk  values  if  continued  from  lower  to  higher  temperatures.  The  thermal 
conductivity  of  bulk  tin  is  1050  W/m-K  in  [0  1  0]  direction  and  650  W/m-K  in  [0  0  1]  direction  at  10  K 
and  supposedly  decreases  gradually  from  22  K  to  room  temperature,  after  which  as  the  temperature  is 
increased  the  thermal  conductivity  experiences  a  plateau  and  drops  precipitously  around  500  K,  at  which 
it  changes  from  solid  tin  to  liquid  tin,  following  which  it  rises  a  little  with  the  increase  in  temperature  and 
then  again  drops  down  to  its  lowest  value  of  10.1  W/m-K  at  1473  K.  The  increase  of  the  total  thermal 
conductivity  in  the  liquid  state  is  explained  by  the  fact  that  the  increase  of  the  kinetic  energy  of  the 
electrons  with  increasing  temperature  is  not  upset  by  an  increase  in  the  scattering  of  the  electrons  by  the 
thermal  vibrations  of  the  atoms  which  apparently  in  the  liquid  state  is  little  dependent  on  the  temperature. 

D.  THIC  NESS  DEPENDENT  THERMAL  CONDUCTIVITY  OF  TIN  FILMS 


The  measured  thermal  conductivity  of  500  nm  Sn  film  at  room  temperature  is  46±4.2  W/m-K,  higher 
than  the  measured  thermal  conductivity  of  36±2.88  W/m-K  in  100  nm  Sn  film.  Moreover,  Choi  et  al.[15] 
reported  thermal  diffusivity  equivalent  to  thermal  conductivity  of  60.3  W/m-K  for  2000  nm  thick  Sn  film 
at  room  temperature.  This  indicates  that  similar  to  other  thin  films,  size  effect  is  observed  in  the  thermal 
conductivity  of  Sn  at  room  temperature,  and  the  thermal  conductivity  decreases  with  thickness. 


The  size  reduction  factor  in  metallic  thin  films  was  first  explained  by  Fuchs-Sondheimer  (FS)  model 
[7]  by  solving  the  rigorous  Boltzmann  transport  equation.  In  this  model,  only  surface  boundary  scattering 
is  considered  through  boundary  conditions  for  solving  the  BTE.  The  size  effect  is  modeled  by  a  factor 
which  accounts  for  the  observed  property  in  thin  films  normalized  to  the  value  in  bulk  as  follows, 
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where  p  is  the  specularity  parameter  that  represents  electron  diffuse  scattering  at  the  boundary  surface 
with  a  probability  of  (1  -p)  (p=0.5  means  partially  diffuse  scattering  and  is  generally  used);  k0  is  the  ratio 
of  film  thickness  to  the  mean  free  path,  and  is  the  inverse  of  the  Knudsen  number. 

A  more  elaborate  model  derived  from  BTE  was  proposed  by  Mayadas  and  Shatzkes  (MS)  [8],  which 
includes  the  grain  boundary  scattering  along  with  the  surface  boundary  scattering  in  a  polycrystalline 
metallic  film.  The  size  reduction  factor  based  on  MS  theory  can  be  expressed  as  follows: 
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where  in  H{t ,  </>)  =  1  +  d  /  Costfi^jl  -  j/2  ,  k0  =  ^  and  a  =  ^  _  Ry  d  is  the  film  thickness  and  d  is 

the  average  grain  size,  p  is  the  specular  reflection  coefficient  and  r  is  the  electron  reflection  coefficient  at 
the  grain  and  twin  boundaries  combined  in  this  work. 


Based  on  the  above  model,  Qiu  and  Tien  (QT)  [31]  proposed  an  approximate  solution  when  the  film 
thickness  and  the  grain  size  are  not  too  small  compared  with  the  electron  mean  free  path,  as  follows 
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The  measured  thermal  conductivity  of  500nm  and  lOOnm  Sn  films  in  this  work  at  room  temperature 
and  the  measured  thermal  conductivity  of  2000nm  in  an  earlier  work  by  Choi  et  al.  [15]  are  compared 
with  the  predictions  of  the  Fuchs-Sonheimer  (FS)  model,  Mayadas  and  Shatzkes  (MS)  model  and  the  Qiu 
and  Tien  (QT)  model.  While  the  FS  model  is  not  in  agreement  with  the  measured  thermal  conductivity, 
the  QT  model  is  in  agreement  within  15%  error,  while  the  MS  model,  which  includes  both  the  grain 
boundary  scattering  of  electrons  as  well  as  the  inelastic  scattering  at  the  boundary  surface,  is  in  excellent 
agreement  with  the  measured  thermal  conductivity  of  Sn  films  of  thickness  100  nm,  500  nm  and  2000  nm 
at  room  temperature.  This  indicates  that  indeed  both  the  scattering  phenomena  are  responsible  for  the  size 
reduced  thermal  conductivity  in  the  polycrystalline  tin  (Sn)  thin  films. 
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Film  Thickness  jnm) 


Fig.  8  Comaprison  of  experimental  thermal  conductivity  values  of  Sn  films  of  thickness  lOOnm,  500nm  and 
2000nm  with  FS  model,  MS  model  and  QT  model.  The  values  are  normalized  in  the  form  of  a  size  reduction 
factor. 


SUMMARY 

We  report  the  thermal  conductivity  of  500nm  Sn  film  from  55K  to  300K  and  a  lOOnm  Sn  film 
from  40K  to  31  OK.  While  the  thermal  conductivity  of  500nm  thick  film  increases  as  the  temperature  is 
lowered  due  to  increase  in  the  mean  free  path  of  the  heat  carriers,  the  thermal  conductivity  of  lOOnm 
decreases  gradually  as  the  temperature  is  lowered.  The  decrease  in  thermal  conductivity  is  attributed  to 
the  scattering  of  electrons  by  the  grain  boundaries  and  twin  boundaries  as  evidence  by  the  AFM  images 
and  further  confirmed  by  the  excellent  agreement  with  the  Mayadas-Shatzkes  model  and  Qiu-Tien  model. 
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Abstract 

This  paper  reports  dependency  of  specific  heat  and  ballistic  thermal  conductance  on 
cross-sectional  geometry  (tube  versus  rod)  and  size  (i.e.,  diameter  and  wall  thickness),  in 
free-standing  isotropic  non-metallic  crystalline  nanostructures.  The  analysis  is  performed  using 
dispersion  relations  found  by  numerically  solving  the  Pochhammer-Chree  frequency  equation  for 
a  tube.  Estimates  for  the  allowable  phonon  dispersion  relations  within  the  crystal  lattice  are 
obtained  by  modifying  the  elastic  acoustic  dispersion  relations  so  as  to  account  for  the  discrete 
nature  of  the  material’s  crystal  lattice.  These  phonon  dispersion  relations  are  then  used  to 
evaluate  the  specific  heat  and  ballistic  thermal  conductance  in  the  nanostructures  as  a  function  of 
the  nanostructure  geometry  and  size.  Two  major  results  are  revealed  in  the  analysis:  increasing 
the  outer  diameter  of  a  nanotube  while  keeping  the  ratio  of  the  inner  to  outer  tube  radius  (y) 
fixed  increases  the  total  number  of  available  phonon  modes  capable  of  thermal  population. 

Secondly,  decreasing  the  wall  thickness  of  a  nanotube  (i.e.,  increasing  y)  while  keeping  its  outer 
diameter  fixed,  results  in  a  drastic  decrease  in  the  available  phonon  mode  density  and  a  reduction 
in  the  frequency  of  the  longitudinal  and  flexural  acoustic  phonon  modes  in  the  nanostructure.  The 
dependency  of  the  nanostructure’s  specific  heat  on  temperature  indicates  ID,  2D,  and  3D 
geometric  phonon  confinement  regimes.  Transition  temperatures  for  each  phonon  confinement 
regime  are  shown  to  depend  on  both  the  nanostructure’s  wall  thickness  and  outer  radius. 

Compared  to  nano  wires  ( y  =0),  the  frequency  reduction  of  acoustic  phonon  modes  in  thinner 
walled  nanotubes  (y  =  0.96)  is  shown  to  elevate  the  ballistic  thermal  conductance  of  the 
thin- walled  nanotube  between  0.2  and  150  K.  At  20  K,  the  ballistic  thermal  conductance  of  the 
thin-walled  nanotube  (y  =  0.96)  becomes  300%  greater  than  that  of  a  solid  nanowire.  For 
temperatures  above  150  K,  the  trend  in  ballistic  thermal  conductance  inverts.  The  greater  number 
of  phonon  modes  in  nanostructures  with  increased  outer  diameter  and  wall  thickness  is  shown  to 
have  a  larger  contribution  to  ballistic  thermal  conductance  when  compared  to  the  increased 
contribution  from  the  frequency  reduction  of  acoustic  phonon  modes  in  thinner  walled  nanotubes. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

Specific  heat  and  thermal  conductivity  in  bulk  crystalline 
solids  are  intrinsic  properties,  which  are  independent  of  the 

1  Author  to  whom  any  correspondence  should  be  addressed. 


crystal  dimensions.  This  is  understood  to  be  true  even  at  low 
temperatures  when  phonon  wavelengths  are  large.  Specific 
heat  and  thermal  conductivity,  however,  become  dependent  on 
the  geometry  of  the  structure  when  the  size  of  the  structure 
interferes  with  the  existence  and  propagation  of  energy 
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carriers,  predominantly  the  dominant  phonon  wavelength  in 
non-metallic  materials  [1-4]. 

The  mechanisms  governing  the  thermal  transport  in 
cylindrical  structures  are  of  particular  interest  in  the  design 
and  development  of  more  efficient  thermal  management/heat 
dissipation  technologies  [2,  3]  and  in  the  design  and 
development  of  next-generation  thermoelectric  materials 
requiring  low  thermal  conductivities  combined  with  high 
electrical  conductance  [5].  Undoubtedly,  both  applications  are 
expected  to  benefit  from  further  mechanistic  studies  of  thermal 
transport  through  cylindrical  nanostructures. 

Literature  is  ripe  with  various  methods  of  thermal 
conductivity  estimates  in  individual  nanowires  [3,  4,  6-8],  as 
well  as  carbon  nanotubes  [1,  9-12].  However,  to  date,  a 
direct  comprehensive  comparison  of  the  specific  heat  capacity 
and  ballistic  thermal  conductance  in  nanowires  with  those 
in  nanotubes  of  varying  geometry  (and  dimensions)  does 
not  exist.  For  example,  in  the  past,  the  thermal  properties 
of  carbon  nanotubes  (CNTs)  have  been  estimated  from 
phonon  dispersion  relations  calculated  from  spring  constant 
and  bond  angle  models  [13].  While  such  atomistic  studies 
have  adequately  accounted  for  diameter  dependencies  of 
the  CNTs,  these  models  become  increasingly  complex  with 
the  addition  of  atoms  to  the  wall  thickness;  the  use  of 
3D  cylindrical  continuum  wave  models  to  represent  phonon 
dispersion  is  an  attractive  alternative  to  investigating  the  effects 
of  geometry  and  size  on  specific  heat  and  ballistic  thermal 
conductance  in  nano  wires  and  nanotubes.  These  phonon 
dispersion  relationships  provide  insights  into  the  allowable 
phonon  frequencies  capable  of  thermal  population  for  a  given 
nanostructure,  and  are  an  invaluable  comparative  tool  for 
understanding  the  intrinsic  differences  between  the  thermal 
properties  of  nanoscale  wires  and  tubes. 

In  the  present  study,  variations  in  available  phonon 
mode  densities  as  a  function  of  the  nanostructure  geometry 
and  size  are  explored  using  acoustic  dispersion  relations 
derived  from  a  continuum-level  3D  elastic  wave  model. 
Estimates  for  the  allowable  phonon  dispersion  relations  are 
obtained  by  modifying  the  elastic  dispersion  relations  so  as 
to  account  for  the  discrete  nature  of  the  material’s  crystal 
lattice  [14].  These  phonon  dispersion  relations  are  then  used 
to  estimate  the  specific  heat  and  ballistic  thermal  conductance 
of  the  nanostructures.  Geometric  effects  on  wave  dispersion 
in  cylindrical  [15]  and  rectangular  nano  wires  [16]  have 
previously  been  studied  using  continuum-level  elastic  models. 
The  consequences  of  these  models  on  phonon  dispersion  and, 
therefore,  specific  heat  and  ballistic  thermal  conductance  has 
only  been  explored  in  the  case  of  cylindrical  nanowires  [14]. 
But,  to  date,  the  effects  of  phonon  dispersion  modification  on 
specific  heat  and  thermal  conductance  due  to  nanotube  wall 
thickness  and  outer  diameter  have  not  been  investigated. 

2.  Modeling 

2.1.  Elastic  dispersion  relations  as  a  function  of  the 
nanostructure  geometry 


tube.  To  derive  the  Pochhammer-Chree  frequency  equation 
from  the  complete  3D  cylindrical  elastic  wave  model,  we  apply 
the  method  outlined  by  Graff  [17].  The  frequency  equation  is 
derived  for  the  nanotube  using  the  assumption  of  stress-free 
boundary  conditions  on  both  inner  and  outer  free  surfaces. 
However,  it  must  be  noted  that  the  stress-free  boundary 
conditions  are  chosen  as  a  first  approximation  to  analyze  the 
effects  of  nanoscale  geometry  on  phonon  confinement  and 
its  effects  on  specific  heat  and  ballistic  thermal  conductance. 
In  some  recent  research  [18-20],  it  has  been  reported  that 
phonon  dispersion  is  altered  by  the  presence  of  surface  stress 
and/or  residual  stress  fields  within  the  nanostructures.  In  this 
regard,  the  existence  of  a  positive  hoop  stress  due  to  bond  order 
loss  of  surface  atoms  in  solid  nanowires  is  likely  to  decrease 
the  frequency  of  the  acoustic  phonon  modes  and  increase  the 
average  phonon  group  velocity  [19].  As  will  be  later  shown, 
these  effects  are  likely  to  elevate  the  specific  heat  as  well 
as  the  ballistic  thermal  conductance  of  the  nanostructures. 
Nonetheless,  the  role  of  cross-sectional  geometry  alone  on 
specific  heat  and  ballistic  thermal  conductance  has  not  been 
investigated  for  nanotubes  of  varying  wall  thickness  and 
diameters,  and  warrants  a  systematic  investigation.  For  this 
reason,  in  the  present  study,  stress-free  boundary  conditions 
are  chosen  to  simplify  the  analysis  and  at  the  same  time  isolate 
the  variables  in  question. 

Due  to  the  complexity  of  the  applied  method,  it  is 
advantageous  to  summarize  the  key  points  of  the  derivation. 
Using  a  Helmholtz  decomposition,  the  wave  displacement 
field  of  a  cylindrical  structure  can  be  described  using  a  scalar 
potential,  (p,  and  a  vector  potential,  H,  such  that 

u  =  V0  +  V  x  H.  (1) 


The  potentials  (p  and  H  satisfy  the  scalar  and  vector  wave 
equations, 

1  920 


V20  = 


v2h  = 


c2  dt2 

i  a2H 

c?  dt2 


(2) 


(3) 


where  c\  is  the  longitudinal  (compressive)  wave  velocity  and 
C2  is  the  transverse  (shear)  wave  velocity,  both  dependent  on 
the  elastic  properties  of  the  material 


c\  = 


E(l-v) 


p(  l  +  v)(l-2v) 


(4) 


<?2 


2p(l  +  v) 


(5) 


Equations  (4)  and  (5)  are  formulated  in  terms  the  density  p, 
Young’s  modulus  E ,  and  Poisson’s  ratio,  v.  For  the  purpose 
of  this  study,  E  =  1.3  TPa  (reflecting  the  value  of  Young’s 
modulus  reported  for  a  single- walled  CNT  in  [21]),  v  =  0.27 
and  p  =  2260  kg  m-3,  are  used  in  all  calculations. 

The  scalar  components  of  the  displacement  vector  u,  are 
described  in  the  cylindrical  coordinate  system  r,  0,  and  z,  and 
are  given  by 

_  90  1  dHz  dHe 

r  dr  +  r  dO  dz 


Dispersion  relations  are  found  by  numerically  solving  the 
Pochhammer-Chree  frequency  equation  for  an  infinitely  long 
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_  1  30  dHr  dHz 
r  dO  dz  dr 


30  i  a 

“z  =  -T-  +  -  —  (rHe) 
dz  r  dr 


1  dHr 
r  dO 


(7) 

(8) 


The  potentials  given  in  equation  (1)  have  general  solutions, 


O  =  f(r)  cos  nO  cos  (cot  +  &z)  (9) 

Hr  =  hr(r)  sinnO  sin(cot  +  kz)  (10) 

Hq  =  he  (r)  cos  nO  sin  (cot  +  kz)  (11) 

Hz  =  h3(r)  cosnO  cos(cot  +  &z)  (12) 

where  the  functions  /(r),  hr(r),  he(r ),  and  /i3(r)  involve 
solutions  to  the  Bessel  equation,  such  that 

f  =  AZn(\a\r)  +  BWn(\a\r)  (13) 

h3  =  A3Zn(\p\r)  +  B3Wnmr)  (14) 

2/1!  =*r-/i0  =2A1Zn+1(|i8|r)  +  251Wrn+1(|i8|r)  (15) 

2/22  =hr+he  =  2A2Zn_l(\p\r)  +  2£2W„_1(|£|r).  (16) 


In  equations  (13)— (16),  A,  #,  Ai,  B i,  A2,  #2,  A3, 
and  #3  are  unknown  constants,  and  Zn  and  represent 
nth  order  ordinary  Bessel  and  modified  Bessel  functions, 
respectively.  The  specifics  of  Bessel  function  choice  are  given 
in  appendix  A.  a  and  /3  contain  the  frequency  co,  and  the 
wavenumber  k,  and  are  given  by 

a  =  yjo)2/c\  —  k 2  (17) 

P  =  yj  co2  /  c?2  —  k2.  (18) 

There  are  four  functions  f,h\,h2,h3  containing  the  eight 
unknown  constants,  and  only  six  boundary  conditions  that  can 
be  applied.  The  property  of  gauge  invariance  [17]  allows  for 
any  one  of  h\,h2,h3  to  be  equated  to  zero  without  loss  of 
generality.  Because  stress-free  boundary  conditions  are  applied 
to  the  inner  and  outer  surfaces  of  the  tube,  stress  field  equations 
are  required  to  determine  all  the  remaining  unknown  constants. 

The  stress  fields  are  obtained  via  the  elastic  constitutive 
relation  (Hooke’s  law), 


&ij  —  hekk^ij  T  2fieij 


(19) 


where  Lame’s  constant  X,  and  the  shear  modulus  /x,  are 
functions  of  E ,  and  v.  The  components  of  the  strain  tensor, 
etj ,  (on  the  right-side  of  equation  (19)),  are  given  by 


err  = 


dur 

dr 


(20) 
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(22) 

(23) 
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(24) 


( due  1  duz\ 

Vdi  +  rWJ 


(25) 


and  are  determined  by  differentiation  of  equations  (6)-(8). 

Substituting  equations  (20)-(25)  into  (19)  yields  the  stress 
field  equations.  The  stresses  <rrr,  arz ,  <Jre  are  of  the  only 
concern  and  together  contain  the  six  unknown  constants,  Ai, 
B\ ,  A2,  B2 ,  A3,  B3,  where  A  and  B  have  been  eliminated  by 
the  application  of  gauge  invariance  [17].  The  tractions  at  the 
inner,  r  =  R[,  and  the  outer,  r  =  R0,  boundaries  are  zero  and 
therefore  the  stress  components  crrr,  arz ,  crre  at  the  inner  and 
outer  boundaries  are  set  to  zero. 

The  Pochhammer-Chree  frequency  equation  is  derived  by 
solving  a  matrix  determinant.  The  elements  of  the  matrix 
determinant  are  the  coefficients  of  the  six  unknown  constants 
(Ai,  B 1,  A2,  B2 ,  A3,  £3)  in  each  of  the  six  stress  field 
equations.  Because  the  unknown  constants  have  nontrivial 
solutions,  setting  the  six-by-six  matrix  determinant  to  zero 
satisfies  the  boundary  conditions.  Solutions  to  the  frequency 
equation  are  values  of  k  and  co  that  satisfy  the  six-by-six 
determinant 


\cij\  =  0  ij  =  1---6.  (26) 

Components  of  equation  (26)  are  given  in  appendix  A. 

Equation  (26)  is  solved  using  a  bisection  algorithm  written 
within  MATLAB®.  An  infinite  number  of  allowable  frequency 
solutions  of  the  type  co(n,m)  =  f(k)  are  obtained.  The  (n,  m) 
solutions  are  the  allowable  harmonics  or  modes  of  vibration 
of  the  structure.  Every  azimuthal  mode  number,  n ,  has  an 
infinite  number  of  corresponding  radial  solutions,  m.  A  total 
of  50  radial  modes,  m,  are  first  calculated  for  each  azimuthal 
mode  number,  n,  for  n  =  0-50.  The  class  of  longitudinal 
and  torsional  modes  are  formed  for  values  of  n  =  0.  Flexural 
modes  are  represented  by  values  of  n  >  1 .  Note  that  all  flexural 
modes  must  be  considered  twice  due  to  v  and  y  flexural  wave 
polarizations.  Normalized  frequencies  ( coR0/c0 )  are  calculated 
for  normalized  wavenumber  (kR0)  values  from  kR0  =  0  to 
kR0  =  (jrR0/a),  where  a  =  1.7835  A.  This  choice  of  the 
effective  lattice  constant,  a ,  will  be  discussed  in  more  detail  in 
section  2.2. 

Figure  1  compares  the  n  =  0  and  1  dispersion  curves 
calculated  via  equation  (26),  for  a  solid  rod  (y  =  0),  a 
thick-walled  tube  (y  =  0.5),  and  a  very  thin-walled  tube 
(y  =  0.96).  Examination  of  figure  1  reveals  the  dependency 
of  the  allowable  wave  modes  on  the  cross  section  geometry. 
First,  a  significant  decrease  in  density  of  the  allowable  modes 
of  vibration  is  observed  with  decreasing  tube  wall  thickness 
(increasing  y).  Secondly,  as  the  outer  radius  of  the  tube  is 
increased  the  density  of  the  available  modes  of  vibrations  is 
observed  to  increase  for  all  wall  thicknesses.  Lastly,  of  the 
four  acoustic  modes  (co  0  as  k  — >►  0),  one  torsional  mode 
(TA)  (figures  1(a),  (c)  and  (e)  long  dashed),  one  longitudinal 
mode  (LA)  (figures  1(a),  (c)  and  (e)  short  dashed),  and  one 
doubly  degenerate  flexural  mode  (FA)  (figures  1(b),  (d)  and  (f) 
short  dashed),  the  normalized  frequencies  and  group  velocities 
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Figure  1.  Elastic  dispersion  relations  calculated  by  solving  equation  (26):  (a)  torsional  and  longitudinal  modes  ( n  =  0)  of  a  solid  wire 
(y  =  0);  (b)  flexural  modes  (n  =  1)  of  a  solid  wire  (y  =  0);  (c)  torsional  and  longitudinal  modes  ( n  =  0)  of  a  thick-walled  tube  (y  =  0.5); 
(d)  flexural  modes  (n  =  1)  of  a  thick-walled  tube  (y  =  0.5);  (e)  torsional  and  longitudinal  modes  ( n  =  0)  of  a  thin-walled  tube  (y  =  0.96); 
(f)  flexural  modes  (n  =  1)  of  a  thin- walled  tube  (y  =  0.96).  The  acoustic  modes  (co  0  as  k  0)  are  marked  by  their  respective  dashed 
lines.  The  normalized  frequency  of  the  longitudinal  acoustic  (LA)  mode  and  the  flexural  acoustic  (FA)  mode  decrease  with  decreasing  wall 
thickness  (increasing  gamma).  The  slope  of  the  torsional  acoustic  (TA)  mode  remains  constant.  The  density  of  allowable  mode  vibrations 
increase  with  both  the  increase  in  wall  thickness  (decrease  in  gamma)  and  the  increase  in  outer  radius,  R0. 


corresponding  to  the  LA  and  FA  modes  are  reduced  with  a 
decrease  in  the  tube  wall  thickness  (increasing  y). 

Figure  2  is  a  schematic  representation  of  low  frequency 
TA,  LA,  and  FA  displacements  of  a  nanotube  having  a  wall 


thickness  characterized  by  y  =  0.96  and  R0  =  2.5  nm.  The 
displacement  fields  are  calculated  at  a  frequency  of  1  THz, 
corresponding  to  a  normalized  frequency  of  ~0.1.  The 
amplitudes  are  normalized  for  a  qualitative  comparison. 
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Figure  2.  Schematic  representation  of  acoustic  mode  displacements  in  y  =  0.96,  5  nm  outer  diameter  nanotube,  (a)  Torsional  acoustic  mode 
(TA);  (b)  longitudinal  acoustic  mode  (LA)  and  (c)  doubly  degenerate  flexural  acoustic  modes  (FA). 


Figure  3.  Comparison  of  n  =  0  longitudinal  acoustic  (LA)  mode 
dispersion  of  five  structures  having  a  varying  ratio  of  inner  to  outer 
diameter,  y.  y  =  0  is  the  solid  wire,  while  y  =  0.96  represents  a 
thin-walled  tube.  The  non-dispersive  co  =  C\k  behavior  is  shown 
(dashed  line)  as  a  reference. 


Figures  3  and  4,  respectively,  show  the  variation  in  the 
acoustic  LA  and  FA  modes  as  a  function  of  the  tube  wall 
thickness.  From  figure  3  we  can  infer  that  all  nanostructure 
geometries  studied  in  the  present  investigation  show  a  nearly 
non-dispersive  LA  mode  (similar  to  co  =  c\k)  for  coR0/c0  < 
1.8  and  kR0  <  1.  These  non-dispersive  waves  exist  only 
for  kR0  <  1  and  are  confined  to  the  structure’s  longitudinal 
z-axis,  since  the  group  velocity  of  these  modes  are  that 
of  a  P-wave,  c\  (dashed  line).  As  the  wavenumber  k  is 
increased,  for  a  fixed  outer  radius  of  the  nanostructures,  the 
thin- walled  tube  (y  =  0.96)  shows  the  maximum  deviation 
from  the  c\  velocity,  followed  by  structures  with  increasing 
wall  thickness.  This  behavior  is  as  expected  since  the 


Figure  4.  Comparison  of  n  =  1  flexural  acoustic  (FA)  mode 
dispersion  behavior  of  five  structures  having  a  varying  ratio  of  inner 
to  outer  diameter,  y.  y  =  0  is  the  solid  wire  while  y  =0.96 
represents  a  thin- walled  tube.  The  inset  shows  quadratic  dispersion 
behavior  for  kR0  <  1.  The  quadratic  low  frequency  dispersion 
behavior  is  nearly  identical  to  all  tube  geometries  and  responsible  for 
the  T0-5  specific  heat  dependence  at  low  temperatures  [14]. 


thinnest  wall  tube  has  the  greatest  average  circumference, 
allowing  preferential  circumferential  mode  population  as  the 
wavelength  is  decreased  (wavenumbers  increase).  In  summary, 
the  LA  mode’s  transition  from  non-dispersive  to  dispersive 
wave  behavior,  seen  in  figure  3  at  kR0  ~  1,  is  indicative  of 
a  transition  from  purely  longitudinal  wave  propagation  down 
the  axis  of  the  structure  to  a  coupling  of  longitudinal  waves 
and  azimuthal  waves  around  the  structure’s  circumference. 

The  FA  mode  exhibits  a  quadratic  behavior,  co  oc  k2, 
for  each  structure  for  kR0  <  1  (figure  4).  This  behavior  is 
in  agreement  with  [21-23],  and  is  unlike  the  linear  flexural 
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dispersion  relationship  obtained  from  the  graphene  zone¬ 
folding  method  for  thin- walled  tubes  [13].  Furthermore,  both 
figures  3  and  4  suggest  that  at  low  wavenumbers,  all  acoustic 
modes  are  nearly  independent  of  cross-sectional  geometry. 
This  behavior  is  to  be  expected,  since,  at  lower  frequencies  the 
waves  are  primarily  confined  to  the  longitudinal  z-axis.  For 
kR0  >  1,  the  frequency  of  the  LA  and  FA  modes  decrease 
with  a  reduction  in  wall  thickness  as  wavenumbers  increase 
past  the  point  for  which  the  wavelengths  solely  occupy  the 
axial  dimension  of  the  tube  or  rod  and  transition  to  a  coupled 
longitudinal  and  circumferential  combined  mode  (figures  3 
and  4). 

2.2.  Approximating  the  phonon  dispersion  relationship 

The  wavenumbers  in  the  elastic  dispersion  relation  represent 
waves  that  vary  continuously  along  the  tube  axis  and  are 
not  representative  of  a  periodic  Brillouin  zone.  To  better 
utilize  the  elastic  dispersion  model  as  an  approximation  for 
phonon  dispersion,  the  continuum-based  frequency  spectrum 
must  be  modified  to  account  for  the  discrepancy  between  the 
discreteness  of  the  crystalline  lattice  and  the  continuity  of  the 
3D  continuum  wave  model. 

The  material  used  in  the  description  of  all  nanostructures 
in  this  study  is  assumed  to  have  a  diamond  lattice.  However, 
it  is  not  possible  to  choose  a  conventional  unit  cell  of  an 
interpenetrating  face-centered  cubic  lattice  structure  so  that  the 
unit  cell  contains  one  atom  per  basis  [24] .  Similar  to  the  work 
described  in  [25],  an  effective  lattice  constant,  a,  is  employed 
by  dividing  the  conventional  unit  cell  lattice  constant  of  the 
diamond  lattice  by  one-half.  This  approximation  provides  an 
effective  simple  cubic  lattice  with  approximately  one  atom 
per  basis  and  is  appropriate  given  the  perspective  of  this 
investigation. 

A  Bom-von-Karman  boundary  condition  is  employed  to 
account  for  the  lack  of  periodicity  in  the  elastic  dispersion 
relation  by  forcing  the  group  velocity  of  the  lowest  frequency 
(n,  m)  modes  to  zero  at  the  far  right  of  the  dispersion  relation, 
&max  =  7 x/a.  To  apply  the  Born- von- Karman  boundary 
condition,  every  (n,m)  dispersion  mode  is  multiplied  by  [14], 

sin  (ka  /2)  /  (ha  /2) .  (27) 

Further  discussion  regarding  equation  (27)  is  provided  in 
appendix  B.  Figure  5  shows  a  comparison  between  a 
typical  elastic  dispersion  relation  (solid  lines)  derived  from 
the  continuum  wave  model,  and  the  phonon  dispersion 
approximation  (dashed  lines)  obtained  by  multiplying  every 
( n ,  m )  mode  by  equation  (27).  The  maximum  wavenumber 
represents  the  edge  of  the  first  Brillouin  zone  of  a  symmetric 
simple  cubic  lattice.  This  modification  creates  an  effective 
standing  wave  within  each  lattice  cell  for  the  lowest  frequency 
modes,  thus  allowing  the  elastic  dispersion  relations  to  better 
approximate  the  phonon  dispersion  relation. 

The  total  number  of  required  phonon  modes  is  approxi¬ 
mately  equal  to  three  times  the  total  number  of  available  lat¬ 
tice  points  in  the  cross-sectional  area  of  the  tube.  The  fac¬ 
tor  of  three  represents  the  three  available  degrees  of  freedom 


Figure  5.  Torsional  and  longitudinal  elastic  dispersion  ( solid  lines), 
and  the  torsional  and  longitudinal  phonon  dispersion  relations 
(< dashed  lines)  are  shown  for  the  case  of  a  thick- walled  nanotube 
(y  =0.5)  having  a  5  nm  outer  diameter.  An  approximation  of  the 
phonon  dispersion  relation  is  acquired  by  modifying  the  elastic 
dispersion  relation  using  equation  (27).  The  modification  of  the 
elastic  dispersion  relation  is  performed  for  each  value  of  y .  Group 
velocities  of  the  low  frequency  modes  approach  zero  at  maximum 
wavenumber.  The  elastic  dispersion  relation  and  the  newly  acquired 
phonon  dispersion  approximations  are  nearly  identical  at  low 
frequencies,  which  according  to  equation  (31)  are  of  greatest 
contribution  throughout  the  temperature  domain  of  0.1-600  K. 


for  each  mechanical  oscillator.  Because  the  system  is  approx¬ 
imated  as  a  simple  cubic  lattice,  the  total  number  of  (n,m) 
modes  available  for  thermal  population  in  the  cross  section  of 
the  structure  should  be  three  times  the  number  of  atoms  in  the 
cross-sectional  plane.  Estimation  of  the  exact  number  and  ar¬ 
rangement  of  atoms  is  avoided  by  limiting  the  investigation  to 
a  temperature  domain  below  600  K.  Since  the  purpose  of  this 
investigation  is  to  examine  phonon  confinement  effects,  the 
maximum  value  of  specific  heat  is  not  of  interest,  and  there¬ 
fore,  knowledge  of  an  exact  number  of  phonon  modes  is  not 
essential.  An  exact  number  of  phonon  modes  would  be  re¬ 
quired  if  complete  phonon  mode  saturation  were  to  take  place; 
this  would  occur  at  the  maximum  value  of  specific  heat.  As 
indicated  in  figure  6,  mode  saturation  would  occur  at  temper¬ 
atures  greater  than  600  K,  which  is  well  past  the  point  where 
specific  heat  values  for  each  tube  geometry  would  likely  con¬ 
verge.  Subsequent  to  the  marked  C  a  T 3  behavior,  the  short 
wavelength  phonons  are  unconfined  by  the  dimensions  of  the 
structure,  and,  therefore,  temperatures  exceeding  600  K  are  not 
of  interest. 

When  calculating  internal  energy  and  specific  heat  it  is 
important  to  have  an  ample  set  of  phonon  modes  available 
for  thermal  population.  As  temperature  increases  the  internal 
energy  of  the  system  activates  a  greater  number  of  phonon 
modes  that  are  higher  in  wavenumber  and  thus  frequency. 
During  the  present  calculations  the  upper  bound  of  ( n,m )  is 
chosen  such  that  the  internal  energy  and  thus  the  specific  heat 
converges  to  a  maximum  at  the  temperature  of  interest.  The 
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Figure  6.  Specific  heat  versus  temperature  for  five  inner  to  outer  tube 
radii  ratios,  y .  Each  cylindrical  geometry  has  a  5  nm  outer  diameter. 
As  tube  thickness  increases,  specific  heat  is  shown  to  collapse  onto 
that  of  a  solid  nanowire.  The  dashed  lines  are  representative  of  T0-5 
‘ID’,  T1-5  ‘2D’,  and  T 3  ‘3D’  phonon  confinement  behavior. 


convergence  indicates  that  a  proper  number  of  phonon  modes 
have  been  accounted  for  in  the  calculations. 


2.3.  Internal  energy 


The  volumetric  internal  energy  in  the  z -direction  is  calculated 
by  summing  the  contributions  from  every  ( n,m )  mode  pair. 


U(T)  = 


At 

1 

A~t 


dk 

2n  exp (Phco^n^m))  —  1 

71/(1  die  hcO(n,m) 

71  exp (Ph(0(n,m))  ~  1 


(28) 


where  At  is  the  cross-sectional  area  of  the  nanotube  wall 
thickness,  and  /3  =  l/k^T. 

The  upper  bound  on  the  integral  in  equation  (28),  n /a,  is 
the  Debye  wave  vector  cutoff  describing  the  outer  edge  of  the 
first  Brillouin  zone  in  the  &-space.  As  also  inferred  from  the 
denominator  of  the  integrand,  a  surplus  of  modes  beyond  the 
maximum  required  ( n,m )  have  negligible  contribution  since 
the  higher  order  modes  are  of  higher  frequency.  Therefore,  as 
long  as  the  internal  energy  converges  to  a  maximum  value,  a 
sufficient  number  of  ( n,m )  modes  have  been  considered. 


2.4.  Specific  heat  and  ballistic  thermal  conductance 


Differentiating  equation  (28)  with  respect  to  temperature  leads 
to  the  expression  for  specific  heat  capacity  per  unit  volume, 


1  d  C  dk 

CiT)=*.  -<£/- 


hco{ 


( n,m ) 


exp(phco(n,m))  -  1 

7T/a  dl  fe2~2 


hZMZ  h(D(n,m ) 


n  (exp (Pha)(n<m))  -  1) 


2  ' 


(29) 


Similar  to  equation  (28),  the  above  expression  does  not  require 
any  externally  imposed  mode  cut-offs  due  to  the  temperature 
dependent  nature  of  the  denominator. 

The  ballistic  conductance  in  the  z -direction  is  the  sum  of 
contributions  from  each  ( n ,  m)  mode, 


g(T)  =  EE  r.  C(n,m)V(ny 


m) 


n  m 

1  1  ^  ^  r,a  dk 


2  kBT2 


EE jf 

n  m 


0  TT  (exp (phco(n,m))  -  1): 


r  V(n,m)  • 


(30) 


In  equation  (30),  U(n>m)  is  the  group  velocity,  the  slope  of  the 
(n,  m)  dispersion  band,  at  wavenumber  k. 

When  hco/k^T  1,  all  three  integrands  in 

equations  (28)-(30)  have  appreciable  contributions  to  the  value 
of  their  respective  integrals.  Therefore,  phonon  frequencies 
satisfying 

co  «  k^T /h  (31) 


are  the  most  significant  contributors  to  the  internal  energy, 
specific  heat,  and  the  ballistic  thermal  conductance  [26].  Most 
importantly,  equation  (31)  provides  the  criteria  for  thermal 
mode  activation. 

Equations  (29)-(31)  reveal  that  specific  heat  and  ballistic 
thermal  conductance  are  primarily  controlled  by  three  phonon 
mode  characteristics  within  a  typical  nanostructure.  First,  of 
all  the  allowable  phonon  modes,  the  lowest  frequency  phonon 
modes  are  the  first  to  thermally  populate  with  increasing 
temperature.  Secondly,  lower  frequency  phonon  modes  are 
the  greatest  contributors  to  internal  energy,  specific  heat,  and 
the  ballistic  thermal  conductance.  Lastly,  the  total  thermal 
contribution  is  the  sum  of  all  the  thermally  activated  phonons 
at  a  particular  temperature. 


3.  Results  and  discussion 

Qualitatively,  geometric  phonon  confinement  effects  can  be 
explained  in  terms  of  the  dominant  phonon  wavelength  of 
the  system  [1].  The  dependence  of  the  dominant  phonon 
wavelength,  A.^  on  temperature  is  given  by  [26], 


~  2nhc^g/k^T  (32) 

where  cavg  is  an  averaged  phonon  group  velocity,  which 
can  be  calculated  from  any  one  or  group  of  (n,m)  modes 
for  any  range  wavenumbers  from  0  to  it /a.  Decreasing 
the  temperature  of  the  system  drives  the  dominant  phonon 
wavelength  of  the  system  higher.  Lower  temperatures 
will  eventually  cause  the  wall  thickness  and  circumferential 
dimensions  of  the  nanotube  to  become  4 frozen-out’  as 
the  dominant  phonon  wavelength  exceeds  the  respective 
characteristic  length  of  the  dimension.  Similar  to  the  elastic 
waves,  larger  phonon  wavelengths  (smaller  wavenumbers)  are 
confined  to  the  length  of  the  tube,  giving  rise  to  the  ID  regime. 
As  temperatures  increase,  the  dominant  phonon  wavelength 
becomes  smaller  and  begins  to  occupy  the  next  smaller 
dimensions  of  the  system,  starting  with  the  circumferential 
length  (2D  regime),  and  followed  by  the  radial  dimensions, 
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i.e.,  the  wall  thickness  or  the  diameter.  Phonon  population 
of  the  wall  thickness  marks  the  3D  regime.  Once  this  occurs, 
with  regards  to  phonon  confinement,  the  structure  is  thermally 
analogous  to  its  bulk  material. 

3.1.  Tube  thickness  variation 

Figure  6  shows  the  specific  heat  capacities  (per  unit  volume) 
for  five  cylindrical  geometries.  As  expected,  with  an  increase 
in  the  wall  thickness  the  thermal  behavior  converges  onto  that 
of  a  solid  nano  wire.  Evidence  of  ID  —>  2D  3D  transitional 
behavior  is  apparent  with  increasing  temperature  (figure  6). 

In  the  low  temperature  limit,  according  to  equation  (31), 
acoustic  phonon  modes  contribute  significantly  to  equations 
(29)  and  (30).  In  the  case  of  nanowires,  the  acoustic  dispersion 
relations  have  been  approximated  with  analytic  forms  for 
LA  (co  =  c\k ),  TA  (co  =  02k),  and  FA  (co  —  c0Dk 2) 
modes,  where  cQ  =  E/ p  and  D  is  the  diameter  of  the 
nano  wire  [14].  Contributions  to  specific  heat  from  each 
aforementioned  phonon  mode  can  be  independently  found 
using  equation  (29).  At  low  temperatures,  the  LA  and  TA 
phonon  mode  contributions  to  specific  heat  exhibit  a  C  a  T 
behavior,  while  the  FA  phonon  mode  contribution  shows  a  C  a 
T0-5  behavior  [14].  In  this  study,  FA  phonon  mode  dominance 
in  the  low  temperature  limit  is  present  for  all  geometries  and 
shown  by  a  C  a  T° 5  behavior  (figure  6).  FA  mode  dominance 
is  justified  since  the  quadratic  mode  behaviors  (figure  4  inset) 
are  of  lower  frequency,  and  thus  are  higher  contributing  when 
compared  with  the  LA  and  TA  modes  between  kR0  =  0  and  1, 
which  are  shown  in  figure  3.  Figures  3  and  4  also  suggest 
that  the  available  phonon  modes  at  low  frequencies  are  nearly 
independent  of  the  nanostructure’s  cross  section  geometry, 
especially  for  values  of  kR0  <  1.  Because  the  wavenumber  is 
2n/X,  phonon  modes  in  this  regime  have  wavelengths  greater 
than  the  radial  dimension  of  the  nanostructure.  Consequently, 
low  temperature  specific  heat  and  ballistic  thermal  conductance 
for  both  the  nanowire  and  nanotube  systems  are  driven  by  a 
similar  flexural  phonon  mode  mechanism  that  is  independent 
of  cross-sectional  geometry. 

With  increasing  temperatures,  the  T 0,5  dependent  specific 
heat  regime  is  followed  by  C  a  T 1,5  dependence  (figure  6). 
According  to  [1],  the  2D  mid-temperature  regime  is  a  result 
of  thermally  populated  regions  of  the  dispersion  curves 
kR0  >  1  and  kt  <  1,  where  t  is  the  thickness  of  the 
nanotube.  Figures  3  and  4  indicate  that  this  regime  of  the 
dispersion  curve  is  difficult  to  approximate  with  analytical 
forms.  Past  literature  has  described  the  2D  specific  heat 
regime  as  T 2  dependent  [27];  but  the  described  study  was 
based  on  approximated  phonon  dispersion  relations.  As 
expected,  kR0  >  1  and  kt  <  1  correspond  to  dominant 
phonon  wavelengths  capable  of  populating  the  length  and 
circumferential  dimensions  of  the  structure,  but  too  long  to 
populate  the  wall  dimensions. 

Figure  6  shows  that  the  thin- walled  tube  (y  =  0.96) 
is  the  first  to  exhibit  2D  specific  heat  dependence  as  the 
temperature  is  increased  above  0.1  K.  Because  y  =  0.96 
has  a  larger  mean  circumference  compared  to  the  thicker 
walled  tubes,  its  circumferential  dimension  is  more  readily 


populated  with  phonon  wavelengths  as  the  wavelengths  shorten 
with  increasing  temperature.  Moreover,  the  elastic  dispersion 
relations  shown  in  figures  3  and  4  reveal  that  the  y  =  0.96 
tube  supports  lower  frequency/longer  wavelength  LA  and  FA 
modes  over  a  wider  range  of  wavenumbers. 

Similar  to  the  work  in  [14,  25],  above  20  K,  the 
nano  wire’s  specific  heat  (y  =  0)  transitions  to  a  C  a  T 3 
dependent  regime.  Past  literature  has  fit  specific  heat  to  a 
C  a  T 3  dependence  by  applying  the  commonly  used  Debye 
approximation  [14,  25].  The  3D  Debye  approximation  for 
specific  heat  is  acquired  by  replacing  the  entire  dispersion 
relation  with  one  LA  mode,  co  =  c\k ,  and  two  transverse 
phonon  modes,  co  =  c^k ,  which  is  widely  accepted  as  the 
phonon  behavior  in  a  bulk  crystal.  The  3D  specific  heat  of  a 
nanowire  at  elevated  temperatures  is  derived  by  substituting  the 
aforementioned  bulk  phonon  relations  into  equation  (29).  As 
in  the  described  literature,  a  C  oc  T 3  specific  heat  suggests  that 
the  nanostructure  does  not  experience  phonon  confinement. 

The  results  of  the  present  study  suggest  that  the 
temperature  at  which  the  transition  from  2D  to  3D  behavior 
occurs  is  dependent  on  the  wall  thickness  of  the  nanotubes. 
Figure  6  shows  that  for  each  of  the  nanotube  geometries 
investigated,  except  for  the  case  of  the  thinnest  wall  tube 
(y  =  0.96),  the  specific  heat  converges  to  the  T 3  behavior  of 
the  nanowire.  Interestingly,  the  y  =  0.96  nanotube  is  observed 
to  be  exempt  from  the  T 3  transition.  At  600  K  the  dominant 
phonon  wavelength,  k&,  is  estimated  using  equation  (32)  to  be 
approximately  1.67  nm.  Because  the  y  =  0.96  nanotube  has 
a  1  A  wall  thickness,  the  phonons  in  the  radial  dimension  of 
the  nanotube  (at  600  K)  are  incapable  of  thermal  activation 
because  their  wavelengths  are  larger  than  the  wall  thickness 
of  the  tube. 

Figure  7  shows  the  calculated  ballistic  conductance  of 
cylindrical  geometries  y  —  0, 0.5, 0.75, 0.9,  and  0.96 

normalized  by  the  value  of  quantum  conductance  gQ  = 
Ttk^T /6h  [28].  These  geometries  correspond  to  wall 
thicknesses  of  a  solid  nanowire,  1.25  nm,  0.625  nm,  0.25  nm 
and  1  A,  respectively.  The  similarity  in  conductance  at  low 
temperatures  confirms  that  the  same  transport  mechanisms 
exist  in  nano  wires  and  nanotubes  regardless  of  wall  thickness. 
Figures  3  and  4  verify  that  each  structure’s  geometry  supports 
comparable  mode  vibrations  at  low  frequencies.  Therefore,  the 
low  temperature  thermal  behavior  of  nanowires  and  nanotubes 
are  conclusively  identical. 

According  to  figure  7,  at  temperatures  between  0.3  and 
150  K,  ballistic  conductance  is  seen  to  increase  with  decreasing 
wall  thickness.  At  20  K,  the  thermal  conductance  of  the 
thinnest  walled  tube  (y  =  0.96)  is  300%  greater  than  that 
of  the  solid  nanowire  (y  =  0).  Reducing  the  nanotube  wall 
thickness  decreases  the  frequency  of  the  LA  and  FA  phonon 
modes,  thereby  increasing  the  integrand  of  equation  (29),  and 
increasing  the  value  of  thermal  conductance.  Interestingly,  as 
the  temperature  increases  past  150  K,  the  trend  in  conductance 
changes.  At  300  K,  the  ballistic  thermal  conductance  of  the 
thinnest  walled  nanotube  shows  a  50%  decrease  compared 
to  the  solid  nanowire  (figure  7).  Plausibly,  because  the 
number  of  thermally  populated  phonon  modes  increases  with 
increasing  temperature  (via  equation  (30)).  It  is  possible 
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Figure  7.  Normalized  ballistic  conductance  versus  temperature  for 
five  cylindrical  geometries.  The  conductance  is  normalized  by  the 
quantum  of  thermal  conductance,  gQ  =  n2k^T /3h  [28]. 
Conductance  values  converge  to  4 g0  at  low  temperatures  due  to 
contributions  from  the  four  acoustic  modes. 


Figure  8.  Specific  heat  versus  temperature  for  thick- walled 
nanotubes  (y  =0.5)  having  varying  outer  diameter.  The  dashed  lines 
are  representative  of  T0-5  ‘ID’,  TL5  ‘2D’,  and  T 3  ‘3D’  phonon 
confinement  behavior. 


that  the  temperature,  above  150  K,  is  high  enough  so  that 
the  effect  of  the  increased  number  of  contributing  phonon 
modes  supersedes  the  effect  of  the  highly  contributing  lower 
frequency  phonon  modes  which  were  understood  to  play  a 
dominant  role  at  lower  temperatures. 

3.2.  Outer  diameter  variations 

Figure  8  shows  the  calculated  specific  heat  for  a  thick-walled 
nanotube  (y  =  0.5)  as  a  function  of  outer  diameter.  Specific 
heat  is  shown  to  decrease  with  increasing  outer  diameters  at 
low  temperatures.  This  reduction  of  specific  heat  in  nanotubes 
with  increasing  outer  diameter  follows  a  trend  similar  to  those 
observed  in  nano  wires  [14,  25]. 

With  increasing  temperatures,  larger  diameter  nanotubes 
having  the  same  value  of  y  are  first  to  deviate  from  a  ID 
specific  heat  regime  (figure  8).  A  similar  phenomenon  is  seen 
by  increasing  y  for  a  nanotube  having  a  fixed  outer  diameter 
(figure  6).  These  results  suggest  that  increasing  the  outer 
diameter  of  a  fixed  y  nanotube  is  mechanistically  similar  to 
keeping  the  outer  diameter  of  the  tube  fixed  and  decreasing 
the  wall  thickness.  In  both  cases,  the  mean  circumference  of 
the  nanostructure  increases.  As  previously  discussed  in  the 
case  of  the  elastic  dispersion  curves  (section  2.1),  a  cylindrical 
nanostructure  with  an  increased  mean  circumference  is  more 
apt  to  coupled  longitudinal  and  circumferential  phonon 
population  with  increasing  temperatures,  and  will  therefore  be 
the  first  structure  to  enter  a  2D  phonon  confinement  regime. 

Figure  9  indicates  that  increasing  the  diameter  of 
a  nanotube  with  fixed  y  increases  its  ballistic  thermal 
conductance.  When  compared  to  small  outer  diameter 
nanotubes,  the  larger  20  nm  outer  diameter  nanotube  with  the 
same  value  of  y  (y  =  0.5)  begins  to  show  elevated  thermal 


conductance  values,  at  approximately  2  K.  An  approximate  one 
order  of  magnitude  increase  in  thermal  conductance  between 
the  20  nm  outer  diameter  nanotube  and  the  5  nm  outer 
diameter  nanotube  is  observed  at  200  K.  The  inversion  of 
the  temperature  trend  in  conductance  shown  in  figure  7  is 
not  present  in  figure  9.  According  to  figure  1,  increasing 
the  outer  radius  of  the  nanostructure  creates  a  greater  number 
of  available  (n,m)  phonon  modes  within  a  set  frequency 
range  for  a  fixed  value  of  y  (y  =  0.5).  Thereby,  increased 
temperatures  increase  the  number  of  contributing  phonon 
modes  and  provide  a  greater  number  of  (j n,m )  contributions 
to  the  summations  in  equations  (28)-(30).  This  explains  the 
elevation  in  thermal  conductance  with  temperature  for  the 
larger  radius  nanotubes. 

It  is  worth  discussing  the  absence  of  a  trend  inversion 
in  ballistic  conductance  in  figure  9,  which  is  present  in 
figure  7.  Unlike  y  variations,  changing  R0  does  not  alter 
the  phonon  group  velocities  within  the  nanostructure,  or 
lower  the  available  ( n,m )  phonon  frequencies.  According 
to  equation  (31),  of  the  few  thermally  populated  phonons 
in  the  thin-walled  structure,  the  lower  frequency  acoustic 
phonons  have  a  higher  contribution  to  internal  energy  when 
compared  to  the  acoustic  phonons  in  the  relatively  thick- walled 
nanostructures.  Therefore,  thinner  walled  nanostructures  of  the 
same  outer  diameter  are  shown  to  have  a  higher  conductance 
between  2  and  150  K.  As  temperature  increases,  the  lower 
frequency  phonons,  which  have  a  greater  energy  contribution 
per  mode,  must  compete  with  the  sum  contribution  from 
the  growing  number  of  total  (, n,m )  phonon  modes,  which, 
although  having  less  contribution  per  mode,  are  more 
numerous  in  thicker  walled  nanostructures.  Because  increasing 
R0  does  not  lower  the  frequency  of  available  phonons,  the 
contribution  of  each  phonon  to  the  nanostructure’s  thermal 
conductance  remains  the  same,  therefore,  variations  in  R0 
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Figure  9.  Normalized  ballistic  conductance  versus  temperature  of  a 
thick- walled  nanotube  (y  =  0.5)  having  varying  outer  diameter.  The 
conductance  is  normalized  by  the  quantum  of  thermal  conductance, 
g0  =  7 i2k^T /3 h  [28].  Conductance  values  converge  to  4 gQ  at  low 
temperatures  due  to  contributions  from  the  four  acoustic  modes. 


while  maintaining  a  fixed  value  of  y  (figure  9)  do  not  produce 
the  trend  inversion  seen  in  figure  7. 

4.  Summary 

Complete  dispersion  relations  obtained  from  solving  the 
Pochhammer-Chree  frequency  equation  are  shown  to  be  a 
useful  tool  for  investigating  the  thermal  property  dependence 
on  geometry  and  size  of  elastic  non-metallic  crystalline 
nanotubes.  While  past  literature  has  extensively  evaluated 
ID  to  3D  transitions  in  nanowires  with  temperature,  this 
work  captures  similar  trends  in  cylindrical  nanotubes  with 
varying  geometry.  The  results  of  the  study  suggest  that  at 
low  temperatures  similar  thermal  transport  mechanisms  exist  in 
nanowires  and/or  nanotubes  regardless  of  their  wall  thickness 
or  outer  diameter.  Elastic  dispersion  relations  confirm  each 
structure’s  capability  to  support  comparable  low  frequency/low 
wavenumber  modes.  The  transitions  unique  to  nanotubes  in  the 
mid-temperature  range  are  found  to  be  dependent  on  geometry, 
and  exhibit  an  initial  T 1,5  dominance  evolving  towards  bulk- 
T 3  behavior  as  the  temperature  is  increased.  An  elevation  in 
the  thermal  conductance  of  thin- walled  nanotubes,  compared 
to  the  same  size  outer  diameter  nanowire,  is  possible  between 
2  and  150  K,  and  is  attributed  to  the  reduction  in  the  frequency 
of  acoustic  phonon  modes.  At  room  temperature,  the  same 
thin-walled  nanotube  shows  a  decrease  in  ballistic  thermal 
conductance  when  compared  to  the  nanowire  comprised  of  the 
same  material. 
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Appendix  A 

c\\  =  {2  n(n  —  1)  —  (/32  —  k2)a2}Zn{a\a) 

+  2kiotiaZn+i(ot\a) , 

C\2  =  2kpia2Zn{P\a)  —  2  ka(n  +  \)Zn+i(pia), 

C13  =  -2  n(n  l)Zn(p\a)  +  2X2n/3iaZn+i(/3ia), 

c\4  =  {2n(n  —  1)  —  (p2  —  k2)a2}Wn(a\a) 

+  2a\aWn+\{a\d), 

c15  =  2k2kP\a2Wn(P\d)  —  2(n  +  \)ka  Wn+i(fi\a), 
c i6  =  —2n(n  -  l)Wn(Pia)  +  2n/3ia Wn+i(Pia), 
c2 1  =  2n(n  —  l)Zn(a\a)  —  2X\na\aZn+\(a\a) , 

C22  =  —kp\a2Zn(pia)  +  2  ka(n  +  l)Zn+\(/3\a), 

C23  =  -{2 n(n  -  1)  -  p2a2}Zn{P\a)  -  2X2piaZn+i(Pla), 
c24  =  2  n(n  —  \)Wn{(X\a )  —  2na.\aWn+\{(X\d) , 
c25  =  ~ X2kPia2Wn(Pia )  +2  ka(n  +  l)W„+iOSia), 
c26  =  -{2 n(n  -  1)  -  p2a2}Wn(Pia)  -  2plaWn+l{pla) , 

C31  =  2nka\Zn{a\d)  —  2X\koi\a2  Zn+\{a\a) , 

C32  =  nf5\aZn(f}\a)  -  (P2  -  k2)a2 Zn+x{Pia) , 
c33  =  —nkaZn(j3\a), 

C34  =  2nkaWn(a,\a)  —  2ka\a2Wn+i{a\a) , 

c35  =  X2nPiaWn(Pla)  -  (p2  -  k2)a2Wn+l(Pia), 

c36  =  —rikaWn(P\a), 

c4i  =  {2 nin  -  1)  -  ( P 2  -  k2)b2}Zn(a\b ) 

+  2\\U\bZn+\((x\b) , 

c42  =  2kP\b2Zn(P\b)  -  2 kb(n  +  ^Zn+^b), 

C43  =  -2 n(n  -  l)Zn(pYb)  +  2X2nplbZn+l(plb), 

C44  =  {2 n(n  -  1)  -  (p2  -  k2)b2}Wn(aib) 

+  2a\bWn+\{a\b), 

c45  =  2)-2kp\b2Wn(P\b)  -  2(n  +  \)kbWn+i(Pib) , 
c46  =  -2 n(n  -  l)Wn(plb)  +  2nplbWn+1(plb), 

C51  =  2n(n  —  \)Zn{ot\b)  —  2X\naibZn+\(otib) , 

C52  =  —kp\b2Zn{P\b)  +  2  kb(n  +  \)Zn+i{Pib), 

c53  =  -{2 n(n  -  1)  -  p2b2}Zn(pxb)  -  2X2P\bZn+\(P\b), 

C54  =  2  n(n  —  l)Wn{pi\b)  —  2noi\bWn+\(oiib) , 

C55  =  -X2kp1b2Wn(plb)  +  2kb(n  +  1  )Wn+1(plb), 
c56  =  -{2 n(n  -  1)  -  pWWnip^b)  -  2 pxbWn+l(Pib), 

C6i  =  2nka\Zn{a\b)  —  2X\ka,\b2  Zn+\(a\b) , 
c62  =  npibZ^Pib)  -  (P2  -  k2)b2Zn+i(Pib), 
c63  =  -nkbZn(pib), 

cm  =  2nkbWn(ct\b)  —  2ka\b2Wn+\(a\b) , 
c65  =  X2nP\bWn(P\b)  -  (p2  —  k2)b2Wn+i  (P\b), 

C66  =  —nkbWn(p\b), 
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where  a  =  R0  and  b  =  Rj.  If  a 2  or  /3 2  is  greater  than  zero, 
Zn  is  the  nth  order  Bessel  Function  of  the  first  kind  Jn.  If  a2 
or  ft2  is  less  than  zero,  Zn  corresponds  to  the  modified  Bessel 
function  of  the  first  kind  In .  Similarly,  Wn  corresponds  to  the 
Bessel  function  of  the  second  kind,  Yn,  or  modified  Bessel 
function  of  the  second  kind,  Kn ,  depending  on  the  sign  of  a2 
or  ft2.  The  constants  X\,  and  A,  2  are  +1  if  Bessel  functions  J 
and  Y  are  used,  and  —  1  if  I  and  K  are  used.  The  reader  is 
referred  to  [17],  for  further  specifics  regarding  the  derivation 
of  the  elastic  frequency  equation. 

Appendix  B 

Equation  (27)  is  designed  to  force  the  group  velocity  of  only 
the  acoustic  modes  (m  =  1  for  n  =  0  and  1)  to  zero  at 
maximum  wavenumber.  While  not  designed  to  do  so,  the 
group  velocity  of  modes  m  =  2  and  3  for  n  =  0  and  1 
are  near  zero.  For  m  >  3  for  n  =  0  and  1,  and  for  all 
m  modes  belonging  to  n  >  1,  modification  of  the  elastic 
dispersion  relation  using  equation  (27)  causes  slightly  nega¬ 
tive  group  velocities  at  maximum  wavenumber,  especially  for 
higher  frequency  modes.  This  has  been  found  to  be  a  math¬ 
ematical  artifact  having  no  physical  meaning.  Using  equa¬ 
tion  (31)  and  considering  the  temperature  domain  of  interest, 
0.1-600  K,  the  highest  thermally  contributing  frequencies  of 
the  dispersion  curve  may  be  estimated  as  those  frequencies  be¬ 
low  71  THz.  Therefore,  using  the  modification  term  creates 
an  effective  ‘Born-von-Karman’  boundary  condition  at  kmax 
in  the  modes  that  contribute  most  to  specific  heat  and  thermal 
conductance  and  is  therefore  a  better  approximation  for  phonon 
dispersion  than  the  elastic  dispersion  curves.  While  the  calcu¬ 
lation  of  thermal  conductivity  would  require  further  estimation 
of  phonon  scattering  effects  and  is  therefore  beyond  the  scope 
of  this  investigation,  it  is  unarguably  a  natural  extension  of  the 
current  findings.  It  has  been  suggested  in  the  literature  that 
using  the  elastic  dispersion  relations  without  further  modifica¬ 
tion  would  only  affect  thermal  conductivity  calculations  by  less 
than  10%  [29]. 
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ABSTRACT 

A  core-shell  elasticity  model  is  employed  to  investigate  the  effect  of  a  nanowire  and 
nanotube’s  increased  surface  moduli  on  its  specific  heat,  ballistic  thermal  conductance,  and 
thermal  conductivity  as  a  function  of  temperature.  Phonon  confinement  is  analyzed  using 
approximated  phonon  dispersion  relations  that  result  from  the  solutions  to  the  frequency  equation 
of  a  vibrating  tube.  The  results  of  the  investigation  indicate  an  approximate  10  %  decrease  in 
lattice  thermal  conductivity  and  ballistic  thermal  conductance  for  a  10  nm  outer  diameter 
crystalline  nanotube  with  an  inner  diameter  of  5  nm  when  the  average  Young’s  Modulus  of  the 
first  three  atomic  layers  on  both  the  inner  and  outer  free  surfaces  are  increased  by  a  factor  of 
1.53.  The  specific  heat  of  the  tube  structure  is  also  found  to  decrease  by  up  to  20  %  at  150  K. 
In  contrast,  a  10  nm  outer  diameter  nanowire  composed  of  similar  material  exhibits  up  to  a  10  % 
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increase  in  thermal  conductivity,  a  25  %  increase  in  ballistic  thermal  conductance,  and  a  48  % 
increase  in  specific  heat  when  its  outer  free  surface  is  acoustically  stiffened  to  the  same  degree. 
Our  simplified  model  may  be  extended  to  investigate  the  acoustic  tuning  of  nanowires  and 
nanotubes  by  inducing  surface  stiffening  or  softening  via  appropriate  surface  chemical 
functionalization  or  coatings. 

Keywords:  Nanowire,  Nanotube,  Thermal  Conductivity,  Surface  Effects,  Surface  Stress, 
Phonon  Confinement 

1.  Introduction 

Nano  wire  and  nanotubes  with  diameters  below  tens  of  nanometers  experience  a 
significant  increase  in  the  surface  area  to  volume  ratio.  As  a  result,  a  large  number  of  surface 
atoms  become  under-coordinated  due  to  bond  order  differences.  With  decreasing  diameter  the 
rising  number  of  under-coordinated  surface  atoms  begin  to  influence  the  overall  material 
properties  of  the  nanostructure1.  Both  the  application  of  stress  as  well  as  the  presence  of 
acoustically  stiffened  lateral  surfaces  are  expected  to  alter  the  mechanical1,  electrical2,  and 
thermal  properties3,  4  of  nanostructures.  While  residual  surface  stresses  are  present  to  some 
degree  in  all  nanostructures,  the  vast  majority  of  theoretical  nano  wire  and  nano  tube  thermal 
conductivity  predictions  have  not  accounted  for  the  actual  surface  moduli5'9,  potentially  leading 
to  significant  errors  in  predicted  values. 

In  general,  the  thermal  conductivity  of  a  nanowire  or  nanotube  is  dependent  on  the 
geometric  confinement  of  phonon  modes  (dictated  by  phonon  dispersion)  as  well  as  scattering 
mechanisms  due  to  phonon-phonon,  phonon-electron,  phonon-impurity,  and  phonon-boundary 
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interactions.  Since  all  these  scattering  mechanisms  are  frequency  and/or  group  velocity 
dependent10,  they  are  expected  to  be  altered  by  phonon  confinement  effects.  Even  in  the 
absence  of  boundary  scattering  effects,  phonon  confinement  is  believed  to  play  a  role  in  the 

O 

modification  of  thermal  conductivity.  In  this  regards,  Zou  and  Balandin  have  shown  that  in 
nanostructures,  geometric  phonon  confinement  effects  lead  to  9  %  reduction  in  bulk  thermal 
conductivity  in  the  absence  of  boundary  scattering,  i.e.,  purely  specular  boundaries,  for 
nanowires  whose  diameter  is  comparable  to  the  phonon  mean  free  path.  In  the  present  study,  we 
show  that  besides  geometrical  effects,  phonon  confinement  is  also  altered  by  the  presence 
acoustic  stiffening  at  the  walls  of  the  nanowires  and  nanotubes. 

In  recent  studies,  continuum  elastic  dispersion  relations  have  been  used  to  show  the 
importance  of  geometric  phonon  confinement  on  specific  heat  and  ballistic  thermal  conductance 
of  single  crystalline  nanowires  with  temperature11"14.  Lu,  et  al.n  have  shown  that  phonon 
dispersion  combined  with  phonon  radiative  transport,  accounting  for  non-equilibrium  phonon 
distributions,  is  capable  of  accurately  predicting  the  thermal  conductivity  of  10  x  20  nm 
rectangular  cadmium  telluride  nanowires.  Chen  et  aln  have  used  Monte  Carlo  simulations  to 
show  that  phonon  dispersion  relations  derived  from  continuum  elasticity  could  be  effectively 
used  to  evaluate  the  thermal  conductivity  of  silicon  nanowires  on  the  sub  100  nm  scale.  Based 
on  these  studies,  geometric  phonon  confinement  effects  are  necessary  when  predicting  the 
thermal  properties  of  nanowires  and  nanotubes  whose  characteristic  lengths  are  of  the  order  of 
the  dominant  phonon  wavelengths  of  the  system. 

Phonon  dispersion  is  understood  to  be  not  only  affected  by  geometry  of  the 
nanostructures,  but  also  on  the  boundary  conditions15"20.  Mukdadi  et  al ,16  have  used  a  semi- 
analytical  continuum-based  analysis  to  show  how  phonon  dispersion  is  altered  in  rectangular 
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nanowires  comprised  of  bilayers  of  two  perfectly  bonded  materials.  Pokatilov  et  al.  have  used 
the  elastic  continuum  approach  to  demonstrate  how  phonon  dispersion  is  altered  when 
rectangular  nanowires  ’  and  thin  films  are  embedded  in  both  acoustically  soft  and 
acoustically  hard  materials.  Their  results  suggest  that  the  sound  velocity  mismatch  at  the  core- 
clad  interface  is  a  potentially  significant  parameter  in  acoustic  tuning  of  composite 
nanostructures.  Using  an  elastic  continuum  analysis,  Zou  et  al.  found  that  in  GaN  thin  films 
coated  with  2  nm  AIN,  partial  phonon  confinement  can  lead  to  an  enhancement  in  thermal 
conductivity,  and  that  the  thermal  conductivity  could  be  tuned  to  higher  or  lower  values  by 
adjusting  the  clad-layer  thickness  while  maintaining  the  thickness  of  the  heterostructure.  Phonon 
confinement  effects  due  to  the  presence  of  stiffer  boundary  coatings  may  be  analogous  to  the 
acoustic  stiffening  of  surfaces.  If  so,  the  thermal  conductivity  would  be  dependent  on  the  surface 
properties  of  the  nanostructure,  i.e.,  surface  moduli,  especially  as  characteristic  dimensions  of  the 
structure  are  reduced. 

The  Bond-order-length-strength  (BOLS)  model  correlates  the  coordination  number  of 
atoms  to  bond-length  and  stiffness  as  well  as  other  properties,  such  as,  localization  of  charge  and 
electro-affinity4.  The  stiffness  of  under-coordinated  surface  atoms  is  in  particular  related  to  the 
elastic  properties,  e.g.,  the  Young’s  modulus  and  density.  Yueng4  utilized  a  continuum  elastic 
dispersion  model  along  with  the  BOLS  model  to  analyze  the  effect  of  surface  stiffening  on  the 
thermal  conductivity  of  Si  nano  wires.  The  study  determined  how  a  stiffened  surface  alters  the 
separate  contributions  from  the  torsional,  longitudinal  and  flexural  phonon  branch’s  thermal 
conductivity  as  a  nanowire’s  outer  diameter  decreases.  The  study  did  not  calculate  a  total 
thermal  conductivity  by  summing  the  contributions  from  each  phonon  branch. 
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In  the  present  work,  we  investigate  the  specific  heat,  ballistic  thermal  conductance,  and 
thermal  conductivity  of  a  single  crystalline  carbon  lattice  nanowire  with  an  acoustically  stiffened 
outer  shell.  These  results  are  compared  with  a  nanotube  having  the  same  outer  diameter  and 
composed  of  the  same  lattice  material  but  with  both  the  inner  and  outer  surfaces  perfectly 
bonded  to  a  thin  acoustically  stiffened  skin.  The  inclusion  of  the  stiffened  skin  on  the  inner  free 
surface  of  a  nanotube  is  found  to  have  interesting  consequences  on  its  predicted  thermal 
properties.  Thermal  conductivity  is  found  to  increase  by  up  to  10  %  when  the  acoustically 
stiffened  skin  is  placed  on  the  outside  of  the  nanowire,  but  decrease  when  the  acoustically 
stiffened  skin  is  placed  on  both  the  inner  and  outer  wall  surfaces  of  the  nanotube.  Such 
capability  to  tune  (decrease)  thermal  conductivity  via  confinement  of  phonon  modes  by  acoustic 
stiffening  in  nanotubes  may  improve  the  efficiency  of  next  generation  thermoelectric  devices. 
On  the  contrary,  a  more  complete  understanding  of  thermal  transport  in  the  presence  of 
acoustically  stiffened  surfaces  in  nanowires  may  lead  to  the  improvement  of  thermal 
management  materials  with  increased  thermal  conductivity. 


2.  The  Core-Shell  BOLS  Model 

In  the  present  study,  we  use  an  isotropic  continuum  wave  model  to  obtain  dispersion 
relations,  i.e.,  relationship  between  mode  frequency  O),  and  wavenumber  k,  in  a  10  nm  outer 
diameter  nanowire  and  nanotube  both  in  the  presence  and  absence  of  acoustically  stiffened 
surface(s).  Phonon  confinement  is  analyzed  by  using  the  co  =  f(k)  dispersion  relationship  as  an 
approximation  for  the  actual  phonon  dispersion  relation  in  the  nanostructure.  The  continuum 
approximation  is  an  accurate  description  of  the  actual  phonon  dispersion  relation  since  the 
average  phonon  mean  free  path  is  of  the  order  of  the  nanowire  or  nanotube  diameter  .  First, 
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calculations  of  specific  heat  and  ballistic  thermal  conductance  are  performed.  Subsequently, 
estimates  of  thermal  conductivity  are  obtained  by  use  of  a  combined  relaxation  time  accounting 
for  phonon  boundary,  impurities,  and  Umklapp  scattering.  In  the  present  study  the  absolute 
magnitudes  of  specific  heat,  ballistic  thermal  conductance,  and  thermal  conductivity  are  of  less 
concern,  instead,  the  focus  is  to  determine  the  relative  changes  in  the  thermal  properties  of 
nanowire  and  nanotubes  with  the  inclusion  of  acoustically  stiffened  surface(s).  The  material 
investigated  in  the  present  study  is  a  single  crystalline  carbon  structure  with  a  diamond  lattice, 
but  the  core-skin  elasticity  model  can  easily  be  extended  to  materials  with  a  variety  of  lattice 
constants  and  material  properties,  including  various  nanolayer  coatings  on  nanotubes  and 
nanowires. 

The  acoustically  stiffened  surfaces  of  the  nanowire  and  nanotube  are  modeled  using  a 
BOLS  model  for  the  elastic  moduli.  In  this  model,  the  Young’s  modulus  of  each  atomic  layer 
increases  as 


a =i,2, 3),  at 

E i 

where  i  is  the  layer  index  (/  =  3  corresponding  to  the  outermost  layer);  Ej  represents  the  Young’s 
modulus  of  the  bulk  material;  w  is  an  adjusting  parameter  related  to  the  material,  for  carbon  w  = 
2.56  *;  at  is  the  bond  contraction  coefficient  and  is  equal  to  the  bond  length  in  the  ith  atomic  layer 
divided  by  the  bond  length  in  the  bulk  material.  The  bond  contraction  coefficient  is  determined 
for  each  of  the  outer  three  surfaces  layers  using  coordination  numbers  (CN)  of  3.91,  6,  and  12 
corresponding  to  the  outer,  middle,  and  inner  three  surface  layers  of  a  face-center  cubic  lattice 
(Fig.  2).  Note,  for  the  outer  surface  layer,  the  CN  has  been  reduced  from  4  to  3.91,  so  as  to 
account  for  surface  curvature1. 
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Eq.  (1)  is  averaged  over  the  last  three  atomic  surface  layers  forming  a  ratio  for  the 
Young’s  moduli  of  two  separate  materials:  a  core  material  having  a  Young’s  Modulus  Ej,  and 
density  pi,  and  likewise  a  shell  structure  on  both  inner  and  outer  surfaces  having  a  Young’s 
Modulus  E2.  and  density  pi  (Fig.  1).  Note  that  the  new  subscript,  2,  refers  to  an  acoustically 
stiffened  material  having  the  elastic  properties  of  the  average  of  the  three  atomic  surface  layers. 

Assuming  the  Poisson’s  ration,  v,  to  be  the  same  for  the  core  and  the  shell  regions,  the 
Lame’  constants  for  the  core  and  shell  regions,  i.e.,  Xu  and  p/.i,  are  proportional  to  E/.i,  where 
X i.i  =  vE ij/(/l  +v)(l -2v) )  and  pi.i  =  Eu/2(l+v).  In  the  present  study,  the  Young’s  modulus  of 
the  core  material  is  chosen  to  be  Ej  =  1.3  TPa,  v  =  0.27,  and  pi  =  2260  Kg/m3.  Subsequently, 
according  to  the  BOLS  model1,  4,  the  acoustically  stiffened  surface  layer(s)  have  moduli  and 
wave  speeds  that  are  1.53  and  1.18  times  that  of  the  core  material,  respectively. 

The  discontinuity  in  the  material  properties  at  the  interface  between  the  skin  and  the  core 
structure  is  handled  by  enforcing  the  displacement  and  stress  fields  to  be  continuous  across  the 
interface.  The  cylindrical  nanostructures  maintain  an  outer  radius,  R0,  of  5  nm.  The  dispersion 
relations  are  solved  for  ratios  of  inner  to  outer  tube  radius,  y,  where  y  =  R/R0.  y  =  0  and  0.5  are 
investigated,  where  y  =  0  represents  a  solid  nanowire,  and  y  =  0.5  represents  a  thick-walled  tube. 
The  skin  thickness,  3,  of  the  BOLS  affected  layer  is  maintained  at  0.534  nm  and  is  assumed  to  be 
equal  to  3d,  where  d=  1.78  A,  the  lattice  constant  of  the  crystal.  The  solid  nanowire  (y  =  0)  has 
an  acoustically  stiffened  layer  on  the  outer  surface,  while  the  nanotube  (y  =  0.5)  is  modeled  with 
acoustically  stiffer  layers  having  similar  elastic  properties  on  both  the  inner  and  outer  surfaces 
(Fig.  1).  In  the  y  =  0.5  case,  the  distances  from  the  tube  axis  to  the  first  and  second  interfaces  are 
denoted  by  a  and  b,  respectively  (Fig.  1).  The  nanowire  case,  y  =  0,  is  modeled  using  Rt  =  a  =  0. 
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For  comparison,  a  similar  wave  model  is  established  for  a  nanowire  and  nanotube  (y  =  0  and  0.5) 
lacking  acoustically  stiffened  BOLS  layers. 

3.  Calculating  Elastic  Dispersion  Relations 

The  frequency  equation  for  the  nanostructures  is  derived  using  a  decomposition  of  the 
displacement  field  into  a  scalar  potential,  <pl2,  and  a  vector  potential,  Hi, 2,  in  both  the  core  and 

shell  materials  denoted  by  subscripts  2  and  1 ,  respectively: 

Ul,2  =V$,2+VXH1,2  (2) 


The  potentials  0  2  and  H12  satisfy  their  respective  wave  equations 


,  _  1  ^  01,2 

Kl~c,  2  dt2 

*1,2 


(3) 


n2  1  a2H12 

V  H  = - — 

1 2  c  2  at2 


(4) 


where  the  transverse  and  longitudinal  wave  speeds  are  defined  as  ct^  —  2  and 

c,12  =  V4,2+2/WA,2  >  respectively. 

According  to  Eq.  (2),  the  scalar  components  with  respect  to  a  cylindrical  coordinate 
system  of  the  displacement  field  in  the  core  and  shell  materials  can  be  described  as  (Fig.  1): 


001 2  10^12 

W„  =^r^+  '  ' 


dr  r  d0  dz 


(5) 
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Using  the  potentials  we  look  for  solutions  of  the  form  : 


$  2  =  / (r)12  cos  nO  cos (cot  +  kz) 

9 

(8) 

Hr  =hr(r)l2smnOsm(cot  +  kz) 

9 

(9) 

He  =  h0(r)l2  cos  nOsin(cot  +  kz) 

9 

(10) 

=h(r)l  2  cosnOcos(cot  +  kz) 

(11) 

The  radially  dependent  functions  in  the  above  product  solutions, /(r),  h,(r),  he(r),  and  h(r) 
are  linear  combinations  of  Bessel  functions.  The  radial  functions  are  decomposed  into  the 
following 


f, 2  "4,2 (1^1,2  I  (  ^1,2  I  ^*) 


(12) 


ha=CaZMr)  +  CuW(^r) 


(13) 


hrl2=-k 


(91,2 


—  ^1,2-^n+l 


(K 


■)+v«i 


(KT) 


(14) 


where  ^412,  ^412, 51>2, 512,  C12,  C12  are  12  unique  unknown  constants.  The  bar  above  the 
unknown  constants  represents  coefficients  to  the  modified  Bessel  functions  (W„)  of  order  n. 
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Because  Wn  is  unbounded  at  r  =  0,  the  barred  constants  are  zero  in  the  core  material  of  the  solid 
nanowire  ( y  =  0).  Since  the  shell  material  is  always  a  cylinder,  even  in  the  y  =  0  BOLS  case, 
unbounded  solutions  are  permitted  at  the  origin  and  therefore  the  hardened  shell  material 
contains  both  unbarred  and  barred  constants. 


The  dispersion  relation,  i.e.,  the  relationship  between  frequency,  co,  and  wavenumber  k,  is 
embedded  in  the  terms  a  and  /?  as  al  2  -  ^jorjc^^2  -k2  and  /3X2  -  -&2  ,  respectively. 

Stress  fields  are  written  via  the  isotropic  elastic  constitutive  relationship  (Hooke’s  Law): 


^7, ,2  A2  ekk J  2  ^ij  +  ^ 1,2  ei 


V\2 


(15) 


The  components  of  the  strain  tensor,  eijX  2 ,  are  given  by 


du„ 


dr 


,  _  1 

'ee'1  r  dO  r 


du. 


dz 
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p  — 
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(16) 


(17) 


(18) 


(19) 


(20) 
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(21) 


and  are  determined  by  differentiation  of  Eqs.  (5)  to  (7). 

The  displacement  field  is  assumed  to  be  a  continuous  function  of  the  radius,  r.  As  shown 
in  Fig.  1,  u(Ri<r<a)  corresponds  to  the  respective  inner  shell  equations,  u(a<r<b )  corresponds 
to  core  material  equations,  and  u(b<r<R0)  corresponds  to  outer  shell  equations.  The  arduous 
task  that  remains  is  to  establish  a  set  of  homogeneous  equations  that  include  the  imposed 
boundary  and  interface  conditions.  This  set  of  homogeneous  equations  can  be  written  as  Ac  =  0, 
where  c  is  the  vector  containing  the  unknown  constants  (up  to  18  elements  from 

[ai2,  Al2,  Bl2,  B12,  Cl2,  C12  ).  The  matrix  A  contains  the  coefficients  to  these  unknown 

constants  and  is  unique  given  the  necessary  boundary  conditions  for  the  cases  when  y  =  0 
(nano wire)  with  acoustically  stiffened  outer  surface,  y  =  0.5  (nanotube)  with  stiffened  inner  and 
outer  surfaces,  and  y  =  0  and  0.5  without  stiffened  surface(s).  Non-trival  co  vs.  k  solutions  are 
found  by  numerically  solving  the  appropriate  frequency  equation,  i.e.,  the  determinant  |A|  =  0. 


3.1.  Boundary  Conditions 

The  appropriate  homogenous  equations  written  from  the  applicable  boundary  and 
interface  conditions  are  presented  below  in  the  following  order:  1)  y  =  0  without  a  stiffened 
surface;  2)  y  =  0.5  without  stiffened  surfaces;  3)  y  =  0  with  a  stiffened  outer  surface;  4)  y  =  0.5 
with  stiffened  inner  and  outer  surfaces.  The  resulting  frequency  equation  is  solved  using  a 
custom  bisection  algorithm  written  in  MATLAB. 
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3.1.1.  Nanowire  without  a  Stiffened  Surface 


The  case  of  a  nanowire  without  a  stiffened  surface  involves  only  the  core  material.  The 
inner  radius  is  Rt  =  a  =  0.  The  outer  radius  is  b  =  R0  =  5  nm.  Moreover,  for  the  unstiffened 
nanowire 

<T,  (A, )  =  ^  (R, )  =  <Tz,  (A )  =  0,  (22) 

and 

m.,( 0)  <°°  i  —  r,0,z.  (23) 

Eq.  (23)  is  automatically  satisfied  by  setting  the  barred  constants  of  Eqs.  (12)  to  (14)  to  zero. 
The  traction  free  surface  condition  at  r  =  R0  is  satisfied  by  constructing  a  3x3  A  matrix  having 
each  row  corresponding  to  an  equation  from  Eq.  (22),  and  each  column  corresponds  to  a 
coefficient  from  one  of  the  3  unknown  constants,  A  /,  B/,  and  C/. 

3.1.2.  Nanotube  without  Stiffened  Surfaces 

The  inner  radius  of  the  nanotube  a  =  R{  =  2.5  nm  while  the  outer  radius  b  =  R0  =  5  nm. 
As  in  the  case  of  the  nanowire,  equations  which  correspond  to  only  the  core  material  remain. 
The  boundary  conditions  are  written  as: 

(R0 )  =  a ^  (R0 )  =  <7rz  (Ra)  =  0  (24) 

and 
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The  inner  and  outer  traction  free  surfaces  at  r  =  R,  and  r  =  R0  are  satisfied  by  constructing 
a  6x6  A  matrix  having  each  row  corresponding  to  an  equation  from  Eq.  (24)  and  Eq.  (25),  and 
each  column  corresponding  to  a  coefficient  from  one  of  the  6  unknown  constants, 

Ax,  Ax,  Bx,  Bl,Cl,C j .  Barred  constants  are  now  permitted  for  the  nanotube  since  unbounded 
solutions  are  of  no  concern  at  r  =  0. 

3.1.3.  Nanowire  with  a  Stiffened  Outer  Surface  (Outer  BOLS  Skin  Layer) 

The  outer  most  surface  of  the  stiffened  shell  material  is  a  free  surface,  and  therefore 
traction  free  boundary  conditions  apply4: 

<Tr2  (K)  =<rra2  (■ Ro )  =  0„2  (*0)  =  0  (26) 

The  above  equations  represent  3  rows  of  the  A  matrix  having  6  unknown 
constants  A2,A2,B2,B2,C2,C2.  The  coefficients  to  these  constants  form  the  first  6  columns  of 
the  9  x  9  A  matrix.  At  r  =  b,  the  tractions  at  the  interface  are  continuous  and  written  as, 


<Jrr2(b)-crrri(b)  =  0,  (27) 

°re2(b)-(Tret(b)  =  0,  (28) 

°rz2(b)~(Jrzfb)  =  0.  (29) 

The  first  terms  in  Eqs.  (27)  through  (29)  share  the  same  6  constants  as  Eqs.  (26)  and  therefore 
share  the  first  6  columns.  The  second  set  of  terms  in  Eqs.  (27)  through  (29)  belong  to  the  solid 
core  material  and  contain  the  coefficients  to  the  unbarred  constants  AjBi.Ci,  forming  the  4th-  6th 
rows  and  7th  -  9th  columns  of  the  9  x  9  A  matrix.  The  final  3  rows  of  A  are  formed  from  the 
homogeneous  displacement  equations, 
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uri(b)-urfb)  =  0 

(30) 

Ug2(b)-uei(b)  =  0 

(31) 

uZ2(b)-uzfb)  =  0 

(32) 

The  first  terms  of  Eqs.  (30)  to  (32)  are  associated  with  the  stiffened  outer  surface  and 
compose  of  rows  7-9  and  columns  1-6.  These  A  matrix  terms  contain  the  coefficients  to  the 

constants  A1,  A2,B2,  B2,  C2,  C2 .  The  second  set  of  terms  of  Eqs.  (30)  to  (32)  are  associated 
with  the  core  material  and  make  up  rows  7-9  and  columns  7-9,  corresponding  to  coefficients 
to  the  constants  AX,BX,CX.  The  bounded  condition  at  r  =  0  is  automatically  satisfied  by  the 

condition  that  A\  -  B\  -  Ci  =  0  .  Note  that  since  the  stiffened  surface  is  a  cylinder,  the  barred 
constants  A2,B2,Ci  ^0. 

3.1.4.  Nanotube  with  Stiffened  Surfaces  (BOLS  Skins) 

The  nanotube  with  stiffened  inner  and  outer  surfaces  (BOLS  skins)  lead  to  the  most 
complicated  A  matrix  (18  x  18  matrix)  satisfying  3  traction  free  equations  at  r  =  Ra,  6  interface 
equations  at  r  =  b,  6  interface  equations  at  r  =  a,  and  3  traction  free  equations  at  r  =  Ri.  Since, 
the  three  layers  are  comprised  of  hollow  concentric  cylinders,  each  term  in  the  A  matrix  is  a 

coefficient  belonging  to  one  of  the  6  constants,  Al2,  A l2,  Bx2,  Bl2,  CX2,  CX2.  The  constants  for 
the  stiffened  skins  are  considered  twice,  forming  18  columns. 

The  3  homogenous  traction  free  boundary  conditions  of  Eq.  (26)  apply  to  the  outer 
stiffened  surface.  The  interface  at  r  =  b  is  described  by  Eqs.  (27)  to  (32).  The  additional 
interface  between  the  core  material  and  inner  stiffened  surface  is  defined  using  Eqs.  (27)  to  (32), 
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but  with  b  replaced  by  a.  The  innermost  surface  is  defined  by  traction  free  equations  similar  to 
Eq.  (26),  but  with  R0  replaced  with  R,. 

3.2.  Frequency  Solutions  of  |A|  =  0 

An  infinite  number  of  allowable  modes  of  vibration  of  the  structure,  each  defined  by 
co(n,m,k),  can  be  obtained  from  solutions  of  the  appropriate  A  determinant.  The  azimuthal  mode 
number,  n,  is  equal  to  the  number  of  radial  nodal  lines  between  <9  =  0  and  n.  Therefore,  n  =  0 
describes  azimuthally  independent  modes  such  as  the  class  of  longitudinal  and  torsional  mode 
types,  n  =  1  has  one  nodal  line  between  6  =  0  and  n,  and  classifies  the  first  flexural  mode.  n>  1 
describes  higher  order  flexural  modes,  m  is  the  radial  mode  number  equal  to  the  number  of 
nodal  lines  perpendicular  to  an  arbitrary  radius  r,  between  r  =  Rt  and  r  =  Ra.  Because  two 

solutions  are  present  for  n  =  0,  past  analysis  of  wave  propagation  in  cylindrical  structures  have 

22 

separated  the  longitudinal  and  torsional  mode  types  into  separate  solvable  determinants 
Given  the  objective  of  gathering  and  subsequently  integrating  the  appropriate  number  of  phonon 
branches,  this  preprocessing  step  is  not  performed. 

In  the  present  study,  51  azimuthal  dispersion  plots  of  co  vs.  k  are  found  each  having  175 
m  branches;  n  =  0  consists  of  both  the  torsional  and  longitudinal  mode  types,  while  n  =  1  to  50 
represents  the  flexural  modes.  Flexural  modes  in  the  x-  and  y-  directions  are  identical  in  the 
wave  analysis  of  vibrating  isotropic  rods  and  cylinders.  Given  this  symmetry,  flexural  modes  are 
counted  twice. 

4.  Approximating  the  Phonon  Dispersion  Relation 

In  the  present  study,  we  are  interested  in  the  frequency  versus  wavenumber  relationship 
of  waves  whose  wave  fronts  propagate  along  the  longitudinal  z-axis  of  the  rod  or  tube.  The 
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wavenumber  is  the  magnitude  of  the  wavevector,  k  =  kz  with  the  subscript  z  dropped  for 
convenience.  Unlike  dispersion  relations  derived  from  lattice  dynamics,  the  elastic  dispersion 
relations  are  derived  assuming  continuum  theory  and  lack  upper  bounds  on  the  wavenumber. 
For  use  as  the  phonon  dispersion  relation,  the  continuum  based  dispersion  relations  are  modified 
by  multiplying  every  mode  by23 


sin(A:fi?/2) / 

/( kd/2 )' 


(33) 


For  simplicity,  the  conventional  unit  cell  length  is  chosen  to  be  half  the  lattice  constant  of 
diamond,  to  approximate  a  carbon-based  material  with  one  atom  per  unit  cell  13.  While  it  is 
impossible  to  accurately  choose  a  unit  cell  of  a  diamond  lattice  that  contains  one  atom  per  basis, 
the  chosen  lattice  constant  is  deemed  suitable  given  the  perspective  of  this  investigation.  The 
modification  given  by  Eq.  (33),  forces  the  slope  of  the  lowest  dispersion  branch  to  be  zero  at  a 
wavenumber  equal  to  kmax  =  n/d.  If  the  approximated  dispersion  relation  is  reflected  about  the 
frequency  axis,  i.e.,  co(k)  =  co(-k),  the  modified  elastic  dispersion  relation  better  approximates  a 
phonon  dispersion  relation  by  enforcing  a  periodic  boundary  condition  every  2n/d  wavenumbers. 


The  frequency  equations  for  y  =  0  and  y  =  0.5,  have  an  infinite  number  of  (n,m)  solutions. 
The  calculation  of  thermal  properties  is  performed  by  selecting  n  and  m  such  that  the  total 
number  of  phonon  branches  equals  three  times  the  number  of  atoms  in  the  nanotube  or  nanowire 
cross-section  to  account  for  the  three  translational  degrees  of  freedom  of  each  atomic  oscillator. 
All  phonon  branches  are  sorted  from  lowest  to  highest  cutoff  frequency  and  subsequently 
renamed  using  subscript  a  (a  =  1,  2,  3...).  The  maximum  value  of  a  is  chosen  as  7530  for  the 
nanowire  (y  =  0)  and  3897  for  the  nanotube  (y  =  0.5). 
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5.  Calculation  of  Specific  Heat,  Ballistic  Conductance,  and  Thermal  Conductivity 

The  contribution  toward  the  internal  energy  per  unit  volume,  U(T),  from  each  a  phonon 
branch  is  summed24. 


U(T)  =  —  Y  f - fanWc) - dk 

7tA  a  J  exp (j3hco(a,k))-l 


(34) 


where  A  is  the  cross-sectional  area  of  the  nanowire  or  nanotube,  and  (3  =  1 /kBT  .  The  limits  of 

integration  in  Eq.  (34)  are  normally  from  -n/d  to  n/d.  However,  due  to  the  symmetry  of  the  first 
Brillouin  zone,  the  integral  is  doubled  and  taken  from  0  to  n la. 

According  to  the  Bose-Einstein  distribution  function  (denominator  of  integrand  in  Eq. 
(34)),  a  greater  population  of  phonons  occupy  lower  energy  states  at  a  particular  temperature. 
More  specifically,  the  distribution  function,  gives  the  average  phonon  population  of  energy  states 
fid)  at  a  given  temperature  (khT).  With  increasing  temperature  the  distribution  function  permits 
lower  energy  states  (mode  frequencies)  to  be  populated  prior  to  higher  energy  states. 

Phonon  confinement  that  results  in  the  reduction  of  phonon  mode  frequencies,  and  thus 
phonon  branches  (mode  softening),  result  in  a  greater  number  of  populated  lower  phonon  modes 
at  a  particular  temperature25.  The  lower  frequency  phonon  modes,  because  of  their  greater 
phonon  population,  will  increase  internal  energy  (and  thus  specific  heat)  as  compared  to  the  less 
populated  higher  frequency  modes.  Moreover,  the  effects  of  phonon  confinement  on  thermal 
properties  become  more  pronounced  when  the  shifted  phonon  branches  are  those  of  an  already 
lower  frequency  since  the  lower  frequency  phonon  modes  are  occupied  by  a  greater  average 
number  of  phonons.  Therefore  it  is  of  interest  to  identify  which  phonon  branches  are  more 
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susceptible  to  the  inclusion  of  acoustically  stiffened  surfaces  and  or  acoustically  mismatched 
coatings. 

Specific  heat  (per  unit  volume)  is  determined  by  differentiating  Eq.  (34)  with  respect  to 
temperature: 


C(T)  = 


1 

nAkBT 2 


Tt/d 

XJ 


a  o 


Pfofadff  exp(/?#ft?(cr,A:))^ 
(exp  (J3hw(a,  k )))  - 1)2 


(35) 


The  ballistic  conductance  in  the  z-direction  corresponds  to  phonon  transport  occurring 
with  a  mean  free  path  (2)  equal  to  the  length  of  the  structure: 

K,d 

g(T)  =— Yj  {  C(a,k,T)  v(a,k)dk,  (36) 

^  a  o 

where  v(a,k)  is  the  phonon  group  velocity.  The  group  velocity  is  defined  as  the  slope  of  the  a 
phonon  branch  (mode)  at  wavenumber  k.  Ballistic  thermal  conductance  is  dependent  on  phonon 
confinement  effects  as  well  as  the  phonon  group  velocity.  Specific  heat  and  ballistic  thermal 
conductance  reveal  differences  in  geometric  phonon  confinement  (or  simply  geometric  effects) 
via  the  phonon  dispersion  relations  since  these  quantities  do  not  require  an  additional  assumption 
regarding  the  scattering  nature  of  phonons  through  either  a  phonon  relaxation  time  r,  or  a  phonon 
mean  free  path  X,  where  for  an  individual  phonon,  X  =  vr. 

Often  in  many  phonon  dynamic  simulations  an  effective  relaxation  time  of  all  phonons, 
t eff,  can  be  linked  to  various  individual  scattering  mechanisms  using  Mathissen’s  rule 

t  -1  =  t~1  +r-1  +r_1  I37t 

eff  boundary  umklapp  impurity  5  k  ' 

In  the  present  study,  the  boundary  scatter  term,  re//  is  described  by, 
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<ounda,y={l-p)^p-’  (38) 

where  p  is  the  specularity  parameter  ranging  from  0  to  1 .  p  =  1  is  complete  specular  (non- 
resistive)  and  p  =  0  is  complete  diffuse  (resistive)  phonon  scattering.  Specular  scattering  is 
assumed  not  to  resist  phonon  propagation  since  the  overall  phonon  momentum  is  conserved24. 
The  boundary  relaxation  time  is  largely  governed  by  the  length  dimension  term,  D,  taken  as  the 
nanowire  diameter,  or  the  nanotube  wall  thickness.  The  phonon  relaxation  times  corresponding 
to  Umklapp  and  impurity  scattering  are  taken  to  be  of  the  form 

Ck,aPP  =  C\Tco{cc,k)2  e~Cl'T  (39) 


purity  =CM^k)4  (40) 

The  constants  in  the  Umklapp  and  impurity  relaxation  times,  Ci  =  1.73xl0'19,  C2  =  137.3,  and 
C3  =  1.32xl0'45  are  chosen  as  representative  values.  These  values  have  been  previously  used  to 
obtain  the  thermal  conductivity  of  bulk  monocrystalline  silicon5. 


The  thermal  conductivity  of  the  nanowire  or  nanotube  structure,  in  terms  of  the  assumed 
effective  relaxation  time  is  given  by 


k(T)  = 


— y7 '  exp(ffi<»(g,*))  v{a  k)i  dk 
2 npT  a  l  [Qxp(j3hco(o(,k)))-l) 


(41) 


which  is  of  the  form  k  =  CvA  21  ’ 24  but  computed  using  the  summation  of  contributions  from  each 
phonon  branch  integrated  over  the  first  Brillouin  zone. 


6.  Results  and  Discussion 
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Figures  3  and  4  depict  a  portion  of  the  dispersion  branches,  prior  to  the  application  of  Eq. 
(33),  for  the  nanowire  (y  =  0)  and  nanotube  (7  =  0.5)  both  in  the  presence  ( dashed  black  lines ) 
and  absence  {solid  red  lines )  of  the  acoustically  stiffened  BOLS  layer(s).  The  non-dimensional 
frequency  ranges  are  chosen  for  two  reasons.  First,  convenience  in  visualizing  changes  in 
phonon  confinement  due  to  the  presence  of  the  BOLS  layer(s)  given  the  large  numbers  of 
branches  computed,  and  secondly,  to  highlight  the  more  accurate  low  frequency  region  of  the 
dispersion  relation.  The  continuum  approximation  is  most  valid  at  lower  frequencies  where 
phonon  confinement  occurs  due  to  the  dominant  phonon  wavelength  being  comparable  to  the 
diameter  of  the  nano  wire  or  nanotube  21 .  At  temperatures  approaching  the  Debye  temperature, 
full  lattice  dynamic  calculations  are  more  appropriate  for  calculating  high  frequency  phonon 
branches  since  the  dominant  phonon  wavelength  is  of  the  order  of  the  lattice  constant  of  the 
crystal.  Therefore  at  increased  temperatures,  a  continuum  analysis  of  the  crystal  becomes 
decreasingly  valid.  Moreover  when  the  dominant  phonon  wavelength  of  the  system  becomes 
smaller  than  the  diameter  of  the  nanowire,  or  wall  thickness  of  the  nanotube,  phonon 
confinement  effects  vanish  25 .  In  the  present  study,  the  prediction  of  all  thermal  properties  is 
performed  to  a  maximum  temperature  of  600  K,  which  is  well  below  the  Debye  temperature  of 
carbon.  A  comparison  to  experimental  data  would  be  necessary  to  ascertain  the  exact 
temperature  limit  at  which  a  continuum  approximation  of  the  crystal  lattice  is  most  valid,  but  that 
analysis  is  beyond  the  scope  of  the  present  investigation. 

Figure  5  displays  the  total  number  of  phonon  branches  existing  within  a  frequency 
window  having  a  width  of  0.02  0)RJcoi  and  spanning  from  0)RJcoi  =  0  to  80.  The  window  size 
is  arbitrarily  chosen  but  remains  constant  for  each  set  of  dispersion  curves  investigated.  A 
formal  density  of  states  (DOS)  calculation  is  not  performed  since  the  modeled  nanotube  is 
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infinitely  long,  i.e.,  significantly  longer  than  the  phonon  mean  free  path.  A  DOS  calculation 
would  require  a  count  of  the  total  number  of  phonon  modes  within  a  frequency  range  co  and  co  + 
A  co.  But  since  a  crystal  length  is  not  assumed,  the  wavevector  spacing  of  the  phonon  modes  are 
infinitesimally  close  together  instead  of  having  a  2n/L  spacing  .  Therefore,  instead  of  a  DOS 
calculation,  Fig.  5  counts  the  number  of  phonon  branches  versus  individual  modes  to  avoid 
tallying  the  number  of  phonon  modes  within  each  branch.  Figure  5  is  a  useful  tool  for 
identifying  which  geometry  creates  a  greater  number  of  phonon  branches  at  a  particular 
frequency. 

6.1.  Nano  wire  Results 

The  presence  of  a  larger  numbers  of  phonon  modes  and  therefore  branches,  especially  at 
lower  frequencies,  is  expected  to  increase  specific  heat,  ballistic  thermal  conductance,  and 
thermal  conductivity.  In  the  case  of  the  nanowire,  the  outer  stiffened  BOLS  affected  surface 
layer  (the  S/Ra  =  0.107  case)  is  found  to  decrease  the  frequency  of  the  long  wavelength  (low 
wavenumber)  branches  for  normalized  frequencies  a)RJc0i  <  5  (Fig.  3  inset).  Figure  3  also 
indicates  that  higher  order  (larger  a)  small  wavenumber  branches  show  a  transition  to  higher 
frequency  in  the  presence  of  the  acoustically  stiffened  outer  layer.  The  transition  to  higher 
frequency  is  most  noticeable  in  short  wavelength  (large  wavenumber)  modes  and  is  also 
reflected  by  a  reduction  of  phonon  branches  in  10  <  coR0/c0i  <  23  frequency  range  (Fig.  5  inset). 
This  shift  to  high  frequencies  is  likely  due  to  the  general  increase  in  wave  velocities  as  shown  by 
Fig.  3.  The  slight  reduction  in  the  number  of  phonon  branches  for  10  <  a>RJc0i  <  23  likely  result 
in  the  reduction  in  specific  heat,  and  thus  ballistic  thermal  conductance,  and  thermal  conductivity 
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in  the  presence  of  the  acoustically  stiffened  outer  surface  of  the  nanowire  for  temperatures  below 
approximately  140  K  (Fig.  6a). 

Figure  6a  also  suggests  that  the  increase  in  phonon  branches  for  normalized  frequencies 
greater  than  24  t 0RJcoi ,  as  shown  by  Fig.  5  in  the  case  of  the  nanowire  with  the  acoustically 
stiffened  outer  surface,  leads  to  the  significant  increase  in  specific  heat,  ballistic  thermal 
conductance,  and  thermal  conductivity  for  temperatures  greater  than  140  K.  Because  the  total 
number  of  phonon  branches  is  the  same  for  the  nanowire  with  the  acoustically  stiffened  skin 
(BOLS  skin)  and  without  the  acoustically  stiffened  skin,  the  increase  in  the  number  of  phonon 
branches  for  normalized  frequencies  greater  than  24  coR0/c0i  is  due  to  a  shift  of  phonon  branches 
from  higher  to  lower  frequencies.  The  shift  in  frequency  explains  the  increase  in  specific  heat, 
which  is  found  to  increase  by  up  to  48  %  at  460  K.  The  lesser  percent  increase  in  ballistic 
thermal  conductance  and  thermal  conductivity  of  the  nanowire  is  likely  due  to  the  greater 
numbers  of  contributing  high  frequency  phonon  branches  having  near  zero  to  negative  group 
velocities,  similar  to  the  optical  modes  present  in  lattice  dynamic  calculations. 

The  relative  change  in  specific  heat  and  ballistic  thermal  conductance  is  reduced 
following  460  K.  Fig.  5  indicates  a  sharp  decrease  in  the  number  of  phonon  branches  present 
for  coRo/c 0i  >  70  in  the  case  of  the  nanowire  with  an  acoustically  stiff  surface  layer  and  a  gradual 
increase  of  phonon  branches  for  the  nanowire  lacking  the  surface  layer.  The  decrease  in  the 
number  of  BOLS  affected  nanowire  braches  and  the  increase  in  the  number  of  phonon  branches 
of  the  non-BOLS  nanowire  explain  the  decline  in  the  relative  change  of  specific  heat  for 
temperatures  greater  than  460  K  (Fig.  6a).  However,  this  result  is  likely  a  manifestation  of  the 
continuum  approximation,  which  is  less  valid  at  higher  temperatures  and  therefore  has  limited 
credibility.  However,  assuming  the  complete  validity  of  the  continuum  assumption  at  the 
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temperature  range  presented,  the  thermal  conductivity  is  slower  to  experience  then  reduction 
following  460  K  because  thermal  conductivity  is  a  function  of  the  square  of  the  phonon  group 
velocities  (Eq.  41),  which  is  found  to  be  elevated  for  nanowires  with  stiffened  surfaces  (Fig.  3). 

While  it  is  admittedly  difficult  to  analyze  the  physical  meaning  of  an  individual  phonon 
branch,  it  is  worth  noting  that  the  acoustic  phonon  branches  (m  =  0  for  both  n  =  0  and  n  =  1) 
(Fig.  3a  and  3b,  respectively)  experience  a  noticeable  increase  in  frequency  corresponding  to  the 
phonon  modes,  and  thus  in  wave-speed  for  wavenumbers  kR0>  10.  This  equates  to  wavelengths 
smaller  than  1/10  of  the  outer  circumference  of  the  structure  or  smaller  than  2n  times  the  width 
of  the  stiffened  surface  (considering  the  width  of  the  stiffened  surface  layer  is  approximately 
i?o/10).  This  shift  is  plausibly  due  to  the  modes  becoming  confined  to  the  stiffer  outer  surface 
layer  that  has  a  greater  longitudinal  and  transverse  sound  speed  compared  to  the  core  material. 
For  acoustic  modes  ( coR0/c0i  —*■  0  for  kRa  — >  0)  whose  wavelengths  are  comparable  to  or  larger 
than  the  outer  circumference  of  the  structure  ( kR0  <  1 ),  the  phonon  dispersion  relations  are  nearly 
equivalent  (Fig.  3).  This  behavior  is  not  surprising  since  these  waves  are  expected  to  be 
geometrically  confined,  i.e.,  too  ‘large’  to  occupy  the  acoustically  stiffened  BOLS  affected 
region. 

6.2.  Nanotube  Results 

Figure  5  indicates  that  when  the  inner  and  outer  surfaces  of  the  nanotube  ( y  =  0.5)  are 
acoustically  stiffened,  the  phonon  branches  experience  a  shift  to  higher  frequencies.  The  shift  is 
plausibly  due  to  an  average  increase  in  each  modes  group  velocity.  Figure  5  reveals  that  a  more 
subtle  difference  occurs  in  the  number  of  phonon  branches  versus  normalized  frequency  for  the 
nanotube  compared  to  the  nanowire  when  the  surfaces  of  these  structures  become  acoustically 
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stiffened.  However,  inspection  of  Fig.  4  reveals  how  drastically  different  phonon  branch  shapes 
are  created  when  the  stiffened  surfaces  are  attached  to  both  the  inner  and  outer  nanotube.  These 
newly  created  branches  are  seemingly  non-analogous  to  the  branches  present  in  the  non-BOLS 
affected  nanotube  and  readily  appear  for  m  >  3  for  n  =  0,  and  m  >  2  for  n  >  1.  The  only 
similarly  existing  branches  are  the  long  wavelength  ( kR0  >1)  low  frequency  acoustic  branches 
(i coR0/c0i  — ►  0  for  kR0  — ►  0).  These  acoustic  modes  are  geometrically  confined  to  the  length  of 
the  nanotube/nanowire  structures  due  to  their  long  wavelengths  being  larger  than  the  structures 
circumference  or  diameter/wall-thickness.  Due  to  the  confinement  of  acoustic  phonon  modes, 
these  phonons  do  not  recognize  the  presence  of  the  stiffened  surface  layers. 

The  relative  change  in  specific  heat,  ballistic  thermal  conductance,  and  thermal 
conductivity  of  the  nanotube  is  negative  throughout  the  entire  range  of  temperatures  investigated 
(Fig.  6b).  The  higher  frequency  shift  in  the  DOS  of  the  BOLS  affected  nanotube  reduces 
specific  heat  since  higher  frequency  phonon  branches  have  a  smaller  occupation  of  phonons  as 
described  by  the  Bose-Einstein  distribution  function.  A  16  %  reduction  in  specific  heat  at  300  K 
is  predicted  for  a  nanotube  with  acoustically  stiffened  surfaces.  Ballistic  thermal  conductance 
and  thermal  conductivity  do  not  experience  the  same  degree  of  reduction  due  to  their  dependency 
on  the  increased  phonon  group  velocities  which  are  readily  apparent  for  nondimensional 
wavenumber  kR0  >  10  (Fig.  4).  As  discussed,  modes  having  wavenumbers  in  this  regime  have 
wavelengths  capable  of  populating  the  width  of  the  BOLS  affected  stiffened  surface  layer.  With 
increasing  temperature,  the  increased  phonon  group  velocities  attempt  to  recover  the  loss  of 
ballistic  thermal  conductance  and  thermal  conductivity  caused  by  the  reduction  in  available 
phonon  branches.  Near  room  temperature,  the  results  indicate  an  approximate  6  %  decrease  in 
thermal  conductivity  for  the  nanotube  with  acoustically  stiffened  surfaces.  Comparison  of  the 
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trend  between  ballistic  conductance  and  thermal  conductivity  indicates  that  these  two  quantities 
remain  within  2  %  of  one  another  throughout  the  temperature  domain  studied  (Fig  6b).  While 
the  use  of  Eq.  (33)  is  believed  to  increase  the  accuracy  of  the  model,  especially  for  higher 
wavenumber/frequency  modes,  it  is  difficult  to  conjecture  a  valid  hypothesis  based  on  a  2  % 
variation. 

7.  Conclusion 

A  model  to  capture  the  effects  of  increased  surface  moduli  on  specific  heat,  ballistic 
thermal  conductance,  and  thermal  conductivity  of  nanowires  and  nanotubes  is  presented.  Of  the 
three  thermal  properties  investigated,  our  results  indicate  that  the  specific  heat  of  a  nanowire  and 
nanotube  is  the  most  sensitive  thermal  property  to  surface  modification,  i.e.,  surface  stiffening. 
The  thermal  conductivity  of  the  nanowire  and  nanotube  having  acoustically  stiffened  surface(s) 
remained  within  10  %  of  the  thermal  conductivity  of  the  nanostructure  lacking  surface 
modification.  While  a  10  %  variation  may  seem  insignificant,  our  simplified  model  suggests  that 
a  more  advanced  study  of  the  topic  is  warranted  that  would  contain  variations  of  the  elastic 
properties  of  the  core  and  or  coating  materials.  Such  models  might  include  the  effect  of  lattice 
anisotropies,  which  would  permit  the  analysis  of  more  specific  materials,  as  well  as  more 
complex  geometries. 

Unfortunately,  there  are  is  little  to  no  experimental  data  available  of  specific  heat  and 
thermal  conductivity  for  nanowires  and  nanotubes  having  diameters  less  than  5  nm.  Moreover, 
the  effect  of  varying  boundary  conditions  on  individual  nanowires  and  nanotubes  of  these 
dimensions  is  best  left  to  theoretical  comparisons  until  more  advanced  measurement  and 
synthesis  techniques  are  developed. 
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The  most  accurate  model  for  the  acoustic  stiffening  of  surfaces  may  be  molecular 
dynamic  simulations,  where  surface  stiffening  effects  are  implicitly  accounted  for  in  the 
simulations.  However,  the  consequence  of  omitting  surface  modification  may  be  best  captured 
through  the  comparison  of  such  atomistic  models  to  elasticity  theory  as  presented  in  this 
investigation. 

The  present  study  also  supports  the  small  body  of  work  that  investigates  the  use  of 
acoustically  stiffened  and/or  softened  surface  coatings  using  chemical  functionalization  to  tune 
thermal  conductivity  in  nanostructures  with  varying  geometry.  The  results  of  this  study  suggest 
that  surface  modulus  modification  (stiffening)  of  the  inner  surface  of  a  nanostructure  will  shift 
the  available  phonon  modes  to  higher  frequency  thus  reducing  thermal  conductivity.  Similarly, 
based  on  the  acoustic  stiffening  results  of  the  nanowire’s  outer  surface,  it  is  plausible  that  the 
opposite  effect  can  be  achieved,  i.e.,  a  reduction  in  thermal  conductivity,  if  an  acoustically  soft 
material  is  bonded  to  the  outer  nanowire  surface.  Moreover,  considering  these  two  conditions, 
one  might  expect  that  the  lowest  thermal  conductivity  reduction  may  be  possible  using  nanotubes 
with  acoustically  stiffer  inner  coatings  and  acoustically  softer  outer  surface  coatings.  Thermal 
conductivity  reductions  created  by  phonon  confinement  effects  due  to  surface  coatings  will 
undoubtedly  be  compounded  by  increased  phonon  boundary  scattering.  Validation  of  these 
hypotheses  are  currently  underway. 

The  relationship  between  phonon  confinement  effects  and  thermal  properties  is  fairly 
well  understood.  Understanding  how  to  alter  phonon  confinement  seems  to  be  the  next  logical 
step  in  engineering  more  efficient  thermoelectric  and  thermal  management  materials. 
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Figure  1:  Schematic  diagram  of  core  skin  structure  nanotube.  For  all  structures  investigated  Ra  =  5nm  is  maintained.  The 
acoustically  stiff  surfaces  (BOLS  affected  layers)  are  applied  to  the  inner  and  outer  tube  walls  and  has  a  thickness  of  S  =  0.534 
nm  (3  atomic  layers).  The  nanowire  is  similarly  modeled  except  with  Rt  =  a  =  0.  The  models  are  compared  to  nanowires  and 
nanotubes  lacking  stiffened  surfaces  ( a  =  Ru  b  =  R0  for  the  nanotube  and  a  =  Rt  =  0,  b  =  R0  for  the  nanowire). 
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Figure  2.  BOLS  Model  for  bond  contraction  as  a  function  of  atomic  coordination  number.  The  coordination  numbers  are 
assumed  to  decrease  in  the  first  3  atomic  layers  exposed  to  a  free  surface.  Beginning  at  the  surface,  the  coordination  numbers 
used  in  the  present  study  are  approximately  3.91,  6,  and  12  ( blue  circles).  The  3rd  most  atomic  layer  from  the  surface  has  a 
coordination  number  of  12  and  is  therefore  assumed  to  belong  to  the  bulk  medium  of  an  FCC  structure.  C.Q.  Sun,  S.  Li,  and  C.M. 
Li,  The  Journal  of  Physical  Chemistry  B  109,  415  (2005). 
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(A) 


Figure  3.  Elastic  dispersion  relations  for  the  solid  rod  ( y  =  0)  for  the  non-BOLS  affected  case  ( S/R0  =  0)  ( solid  red  lines)  and  the 
BOLS-affected  case  (d  /Ra  =  0.107)  {dashed),  (a)  n  =  0  combination  of  both  torsional  and  longitudinal  modes,  (b)  n  =  1  flexural 
modes.  In  the  BOLS-affected  case,  the  acoustically  stiffened  skin  thickness  is  8  =  0.534  nm. 
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(A) 


(B) 


Figure  4.  Elastic  dispersion  relations  for  the  tube  (y  =  0.5)  for  the  non-BOLS  affected  case  ( S/Ro  =  0)  (solid  red  lines)  and  the 
BOLS-affected  case  ((5  /Ro  =  0.107)  (dashed),  (a)  n  =  0  combination  of  both  torsional  and  longitudinal  modes,  (b)  n  =  1  flexural 
modes.  In  the  BOLS-affected  case,  the  acoustically  stiffened  skin  thickness  is  5  =  0.534  nm. 
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Figure  5.  Number  of  phonon  branches  as  function  of  non-dimensional  frequency  for  a  nanowire  ( red  squares)  and  nanotube 
( blue  circles)  having  a  5  nm  outer  diameter.  The  nanotube  wall-thickness  is  2.5  nm.  The  BOLS  skin  layer  refers  to  an 
acoustically  stiffened  material  bonded  to  the  outer  wall  of  the  nanowire  and  both  the  outer  and  inner  walls  of  the  nanotube.  The 
inset  is  shown  with  the  same  units  and  depicts  the  nanowire  and  nanotube  with  and  without  the  BOLS  affected  skin  layer  (s). 
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Temperature  (K) 


Temperature  (K) 


Figure  6.  Relative  change  in  specific  heat  (A C/C),  ballistic  thermal  conductance  (A G/G),  and  thermal  conductivity  (A K/K)  as 
caused  by  the  presence  of  BOLS  hardened  surface(s)  for  (a)  The  nanowire  (7  =  0);  (b)  The  nanotube  (7  =  0.5). 
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